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EXECUTIVE  SUMMARY 


Investigations  were  conducted  to  Identify  those  processes  that 
control  the  loss  and  movement  of  WP/F,  RP/BR,  and  their  combustion 
products  in  air,  soil,  and  aquatic  environments. 

The  fate  of  RP/BR  will  be  controlled  by  oxidative  transformation. 

In  air,  half-lives  of  1.8  years  were  found  while  in  aerated  water,  the 
half-life  approached  3  years  and  was  found  to  be  dependent  on  both 
particle  size  and  oxygen  concentration.  In  soil,  the  transformation  wac 
limited  by  oxygen  diffusion  and  half-lives  into  the  thousands  of  years 
were  projected. 

WP/F  was  found  to  react  extremely  rapidly  with  air  (t  ^  ~  5  min) 
while  it  was  found  to  be  quite  persistent  as  a  suspension  in  water 
(t  1/2  “  2. A  years).  Dissolution  rates  for  WF/F  in  water  were  meaoured 
to  be  1  x  10"  7  g  hr”1  cm"2.  Soluble  WP  will  be  lost  to  the  atmosphere 
through  volatilization  (t  jy2  "  48  min)  ln  shallow  and  turbulent  streams 
while  both  oxidative  and  hydrolytic  processes  will  equally  control  the 
loss  of  WF  in  deep  water  bodies  (t  jy2  *  A2  hrs).  Humic  materials  in 
natural  waters  were  found  to  have  a  stability  effect  on  soluble  WP.  In 
soil,  WP  will  be  highly  persistent  and  its  transformation  will  be 
limited  by  oxygen  diffusion.  A  minimum  life-time  of  10  years  was 
projected  for  a  1  cnr'  chunk  of  WP  burled  12  cm  beneath  the  surface. 

Longer  life-times  are  projected  at  deeper  depths  and  the  buildup  of 
surface  oxide  layers  will  add  to  the  persistence. 

The  oxidative  transformations  and  combustion  of  WP/F  and  RP/BR  lead 
to  the  formation  of  linear  condensed  polyphosphates,  (pi"P22)>  phosphoric 
acid,  phosphorous  acid,  hypophoaphorous  acid,  and  phosphine.  The  linear 
polyphosphates  hydrolyze  rapidly  at  low  pH  (t  2  “  2  days)  which  is 
expected  in  the  aerosol  phase.  In  a  neutralized  state,  the  polyphosphates 
are  subject  to  slow  chemical  hydrolysis  (t  “  100  days).  Microorganisms 
will  hydrolyze  linear  polyphosphates  in  soils  and  water  with  half  lives  of 


15  days  in  systems  containing  1  x  10^  organisms/ml  (gm) .  Cyclic  meta¬ 
phosphates  (tri-,  tetra-,  and  hex*-)  ate  formed  from  the  linear 
polyphosphates  in  aqueous  solution  and  t. re  stable  to  chemical  hydrolysis 
(t  1/2  "  3  years).  However,  microorganisms  can  utilize  these  compounds  as 
a  phosphate  source  with  a  half-life  of  20  days. 

Phosphorous  acid  is  stable  to  chemical  transformations  under 
environmental  conditions,  however  microorganisms  will  biotransform  this 
chemical  to  orthophosphoric  acid  with  a  half-life  of  25-35  days.  Phosphine 
is  rapidly  transported  to  the  vapor  phase  where  reactions  with  ozone  and 
hydroxyl  radical  will  limit  its  persistence  to  half-lives  of  2-3  hours. 

Using  laboratory  generated  transport  and  tranformatlon  rates 
constants,  computer  models  for  aquatic  (SRI  aquatic  model),  terrestrial 
(SRI  terrestrial  model),  and  atmospheric  (MAER0S)  environments  were 
constructed  to  simulate  the  loss  and  movamant  of  WP/F,  RP/BR  and  their 
combustion  products  under  selected  environmental  senarlos. 
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I  INTRODUCTION 


The  increased  use  of  smokes  and  obscurants  in  military  training 
exercises  has  brought  about  a  concern  for  the  effects  of  both  combusted  and 
uncombusted  smoke  products  on  the  environment  and  human  health.  In 
training  exerciser,  smoke  is  generated  to  simulate  actual  battlefield 
conditions.  This  results  in  potentially  significant  contamination  of 
terrestrial,  aquatic,  and  atmospheric  environments  in  training  areas.  The 
persistence  of  such  smoke  in  these  environments  is  not  well  understood. 
Also,  the  deposition  of  smoke  products  into  the  environment  from  the 
incomplete  combustion  of  smoke  devices  and  because  of  past  and  present 
smoke  materials  manufacturing  practices  of  smoke  has  led  to  further  concern 
as  to  the  fate  of  these  agents  in  the  environment. 

Two  such  agents  used  in  training  exercises  are  white  phosphorus/ felt 
(WP/F)  and  red  phosphorus/butyl  rubber  (RP/BR).  WP/F  is  currently 
manufactured  by  Pine  Bluff  Arsenal  (PBA)  in  Fine  Bluff,  Arkansas.  RP/BR  is 
currently  produced  in  Great  Brltlan  and  shipped  to  PBA  for  incorporation 
into  smoke  devices.  Because  of  a  paucity  of  information  available  to 
estimate  the  persistence  of  these  agents  in  the  environment,  the  magnitude 
of  potential  environmental  effects  is  difficult  to  estimate. 


To  help  Increase  the  knowledge  and  understanding  within  this  area, 
this  report  describes  the  efforts  of  SRI  International  to  assess  the 
dominant  processes  that  affect  the  loss  and  movement  of  WP/F  and  RP/BR  when 
dispersed  into  the  environment  in  both  their  combusted  and  uncombustad 
forms.  After  initial  identification  and  compositional  studies  of  smoke 
components.  Investigational  studies  were  conducted  into  the  transport  and 
transformation  processes  that  affect  the  persistence  of  each  component. 
These  studies  included  solubility  measurements,  dissolution  studies, 

Henry's  Constant  determinations,  sorption  and  octanol-water  partition 
coefficients,  and  volatilization  rate  constant  measurements  to  describe 
physical  transport  processes.  Transformation  processes  studied  were 


chemical  oxidation  In  soil  and  water,  hydrolysis,  biotransformation,  and 
transformation  by  the  tropospheric  agents,  hydroxyl  radical  and  ozone. 

We  conducted  kinetic  Investigations  of  dominant  processes,  from  which 
rate  constants  were  established.  From  these  data,  we  projected  half~lives 
for  specific  components  in  terrestrial,  aquatic,  and  atmospheric 
environments.  These  data  were  then  incorporated  into  computer  simulation 
models  to  project  the  overall  behavior  and  persistence  of  phosphorus  smoke 
products  in  the  environment. 
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II  BACKGROUND 


The  use  of  screening  smoke  in  military  operations  is  effective  in 
visually  concealing  personnel  and  materiel,  screening  targets,  and  creating 
confusion  and  deception  among  enemy  troops*  Two  such  smokes  used  in 
military  tactics  are  WP/F  and  RP/BR.  They  are  usually  deployed  from  field 
artillery,  grenades,  mortars,  and  smoke  pots*  WP/F  consists  of  75-80% 
elemental  phosphorus  solidified  into  a  cellulose  (felt)  matrix*  RP/BR 
consists  of  a  mixture  of  red  phosphorus  and  butyl  rubber  (isobutylenes,  BR) 
mixed  in  the  ratio  of  95  RP  to  5  BR*  RP  is  an  amorphous  material  that  is 
not  readily  characterised.  Although  there  are  many  forms  of  RP,  it  is 
generally  represented  as  polymeric  units. 

The  forma  of  elemental  phosphorus  shown  in  Figure  1  can  enter  the 
environment  through  manufacturing  processes  and  practices,  spillage  of 
material  from  smoke  devices,  and  the  Incomplete  combustion  of  material  from 
smoke  devices.  Also,  the  combustion  of  such  smoke  generating  materials 
produces  a  complex  mixture  of  phosphorus  oxides  and  condensed  linear  and 
cyclic  polyphosphates  that  deposits  on  land  and  water  from  an  aerosolized 
state.  Berkowitz  et  al  (1981)  have  reported  the  compounds  shown  below  as 
resulting  from  UP  and  RP  combustion. 


‘4 

or 

RP 


J2L 


H3P04  +  P4°10  +  P4°6  +  H3P03  +  Ptt3  + 


S  ? 

-v0P~O~POrj~ 
HO  t>H 


Thus,  from  this  partial  list  of  tranformation  products,  it  is  obvious  that 
many  forms  of  phosphorus,  ranging  in  oxidation  state  from  -3  (PHj) 
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FIGURE  1  COMPOSITION  OF  PHOSPHORUS  SMOKES 
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to  0  (P^)  to  +5  <P04"3),  will  be  present  In  the  environment.  For  this 
study,  we  consider  the  transformation  of  phosphorus  species  to  ortho- 
*  phosphoric  add  (HjP04)  to  be  the  final  end  product.  Orthophosphorlc  acid 

or  its  salts  are  ubiquitous  in  the  environment  and  any  contribution  to  the 
,  phosphoric  acid  pool  from  phosphorus  smoke  devices  would  be  infinitesimally 

small  when  dispersed  over  a  wide  area  (aerosol  deposition).  Elevated 
levels  of  phosphates  are  expected  at  specific  points  and  will  undergo 
neutralization  reactions  at  ratas  governed  by  the  buffering  power  of  the 
environment  at  the  point  of  contact. 

In  general,  the  processes  that  contribute  to  the  loss  and  movement  of 
phosphorus  species  in  the  aquatic  environment  are  biotransformation, 
oxidation,  reduction,  hydrolysis,  photoxidation,  volatilization,  and 
sediment  sorption.  A  number  of  these  processes  also  occur  in  the 
terrestrial  environment  and  in  the  atmosphere.  By  Identifying  which 
processes  dominate  in  a  particular  environment  and  establishing  the  pseudo- 
first-order  rate  constant  for  the  proceas(es),  we  can  determine 
environmental  half-lives  for  specific  chemicals  and  estimate  their  . 
persistence.  Specific  uses  of  phosphorus  smokes  or  smoke  products  can  then 
be  modeled  using  computer  simulations  of  selected  environments.  These 
simulations  provide  the  basis  for  environmental  assessments  that  permit  the 
prediction  of  the  persistence  of  chemicals  and  their  stable  end  products. 

We  estimated  the  persistence  of  WP/F  and  RP/BR  and  their  products  in 
aquatic,  soil,  and  air  environments.  We  determined  deposition  velocities 
for  combustion  aerosols  and  modeled  the  aerosol  dynamics  using  computer 
simulations.  From  these  results,  we  predicted  the  time  required  to  convert 
WP/F,  RP/BR,  and  their  combustion  products  in  selected  environments  to 
orthophosphorlc  acid. 
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Ill  CHARACTERIZATION  OF  WP/F  AND  RP/BR  AND  AEROSOL  PRODUCTS 


A.  WP/F 

An  actual  WP/F  155  mm  round,  minus  the  burster  charge  was  obtained 
from  Aberdeen  Proving  Ground,  MD.  This  Item  Is  a  sealed  metal  cannlster 
155  mm  In  diameter  and  200  cm  In  length 4  The  cannlster  was  placed  In  a 
wooden  cradle  submerged  In  water  that  was  contained  In  a  55-gallon  drum 
that  had  been  cut  lengthwise#  The  cannlster  was  opened  under  water  using 
a  hacksaw,  the  felt  wedges  were  pried  from  the  cannlster,  and  both  were 
stored  under  water. 

The  phosphorus  composition  of  the  wedge  was  determined  by  drying 
4-10  go  samples  in  an  inert  atmosphere  for  several  days,  weighing  a  sample 
into  a  glass  apparatus,  and  sealing  the  system  under  vacuum.  Then  the 
apparatus  was  filled  with  100X  oxygen,  sealed,  and  allowed  to  burn  to 
completion.  The  residue  was  dissolved  in  water,  boiled  for  1  hour,  and 
titrated  potentlometrlcally  with  standardized  sodium  hydroxide.  The 
results  of  these  determinations  showed  the  WP/F  to  contain  75.3  ±  3.51 
elemental  phosphorus.  The  high  standard  deviation  is  probably  the  result 
of  a  nonuniform  distribution  of  P4  in  the  felt  matrix.  The  remaining  25X 
of  the  WP/F  is  believed  to  be  cellulose,  although  combustion  studies  (see 
Section  III-C)  could  not  account  for  all  material  being  converted  to  CC^. 

B.  RP/BR 

RP/BR  was  also  obtained  from  Aberdeen  Proving  Grounds.  This  material 
is  produced  as  pellets  and  is  very  spongy  and  hydroscopic.  Analysis  of 
eight  samples  for  elemental  phosphorus  by  combustion  followed  by  boiling  of 
the  residue  in  water  and  potentiometric  titration  yielded  a  phosphorus 
content  of  87.9  A  3. 0Z.  An  additional  3.2%  phosphorus  is  present  as 
surface  oxides,  titratable  by  boiling  the  uncombusted  RP/BR  in  water. 

Thus,  surface  oxides  coat  most  RP/BR  pellets.  Combustion  studies 
showed  the  butyl  rubber  to  amount  to  5.2X  of  the  formulation  (see  Section 
III-C) . 


Samples  of  RP/BR  pellets  (~  0.6  cm  x  0.6  cm)  and  powder  (~  0.09  cm) 
were  submitted  to  Ooicron  Technology  Corporation  (Berkeley  Heights,  Hew 
Jersey)  for  surface  area  measurements.  As  measured  by  a  continuous 
volumetric  method,  the  surface  areas  for  powder  and  pellets  were  0.2 
m^/g  and  0.15  mfyg  ±  20%,  respectively.  A  repeat  of  analysis  using  a 
one-point  BET  method  (Brunaller,  Emmett,  Teller  surface  area 
determination)  resulted  in  a  value  of  0.1  m^/g  for  both  samples,  which 
is  slightly  lower  that  the  volumetric  result  but  within  the  uncertainty 
of  the  measurement. 

Samples  of  RP/BR  pellets  and  powder  were  examined  by  scanning 
electron  microecopy  (SEM).  As  shown  in  Figure  2,  the  surfaces  of  the 
pellets  and  powder  are  similarly  rough,  as  is  suggested  by  the  surface 
area  measurements.  It  Is  not  possible  to  distinguish  between  individual 
BR  and  RP  crystallites  by  this  technique. 

X-ray  fluorescence  (XRF)  analysis  of  both  powder  and  pellets  showed 
the  presence  of  phosphorus  as  the  major  component  and  trace  amounts  of 
silicon.  Carbon  was  difficult  to  detect  by  this  technique. 

C.  RP  Combustion  Studies  to  Generate  C02 

Samples  of  RP/BR  and  WP/F  were  combusted  to  examine  the  effects  of 
oxygen  partial  pressure  on  the  production  of  carbon  monoxide  and  carbon 
dioxide.  The  gaseous  species  were  monitored  by  gas  chromatography  <GC) 
using  a  thermal  conductivity  detector  and  a  0.31  cm  x  1.8-m  stainless 
steel  column  packed  with  80/100  mesh  Carbosphere. 

RP/BR  was  Ignited  with  a  Nichrome  fuse  wire  in  a  92  cc  Monel 
reactor  sealed  with  a  Teflon  washer  at  different  oxygen  pressures.  The 
solid  residue  following  combustion  was  dissolved  in  water,  heated  to 
boiling  for  1  hour,  cooled,  and  titrated  with  standard  0.1  K  NaOH  to  the 
pKj  for  H3P04  (pH  4. A). 

Figure  3  shows  the  formation  of  C02  and  CO  and  the  residual  oxygen 
following  RP/BR  combustion  under  different  oxygen  pressures.  As  shown, 
CO  passes  through  a  maximum  but  la  always  an  order  of  magnitude  less 
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(a)  SEM  Photograph  of  RP/BR  Pellets  (230x) 


(b)  SEM  Photograph  of  RP/BR  Powder  (220x) 
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FIGURE  2  SCANNING  ELECTRON  MICROGRAPHS  OF  RP/BR  PELLETS  (TOP) 
AT  230x  AND  RP/BR  POWDER  (BOTTOM)  AT  220x 
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than  CO2  at  tha  same  initial  O2  pressure.  As  the  O2  pressure  Increases, 
the  oxidation  of  carbon  In  BR  (5.2%)  proceeds  nearly  quantitatively 
(97%)  to  C02. 

As  shown  in  Table  1,  we  calculated  the  percent  conversion  of  RP  to 
!>4^10  based  on  the  oxygen  balance  shown  in  Figure  3.  The  oxygen 
consumed  is  calculated  from  the  difference  between  02  (initial  02 
concentration,  abscissa)  and  the  02  curve.  If  the  oxygen  consumption  is 
corrected  for  C02  formation  (CO  can  be  neglected),  the  net  value  is  that 
which  has  combined  with  RP.  Because  we  can  calculate  the  02  that  is 
theoretically  required  to  oxidize  the  available  IP  (87.9%)  to  P4O1Q,  we 
can  estimate  the  extent  of  this  conversion. 


As  shown  in  Figure  4,  the  oxidation  of  RP  in  RP/BR  to  P^O^q 
proceeds  linearly  with  the  initial  oxygen  pressure.  The  incomplete 
conversion  at  a  stoichiometric  oxygen  pressure  (88.5%  conversion)  and 
below  la  probably  due  to  a  protective  shield  of  P^O^q  that  forms 
initially  on  the  RP  surface.  Also,  our  calculations  do  not  take  into 
account  the  presence  of  lower  oxides  such  as  P40g>  which  undoubtedly 
form  during  Incomplete  combustion. 


The  combustion  results  for  RP/BR  are  summarized  in  Table  2  in  terms 
of  RP  and  BR  conversion  and  initial  oxygen  pressure. 


RP  conversion  is  generally  lower  than  the  oxidation  of  BR. 
However,  the  starting  material  contains  an  appreciably  higher  fraction 
of  RP  (RP/BR  »  16.9).  In  11%  excess  02,  the  oxidation  of  both 
components  is  nearly  equal  and  close  to  100%. 


Using  GC/AFID  methods,  we  made  some  preliminary  investigations  of 
phosphorus-containing  apecies  being  generated  in  the  combustion 
process.  Both  P4  (from  the  depolymerization  of  red  phosphorus)  and 
phosphine  (PH3)  were  observed.  The  latter  product  probably  arises  from 
the  disproportionation  of  P40g  in  the  presence  of  moiBture,  as  shown  in 
Equation  1. 
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Table  2 


PERCENT  CONVERSION  OP  BR  AND  P/ 
AS  A  FUNCTION  OF  OXYGEN  PRESSURE 


On  Pressure 

Percent 

Conversion 

<*tm) 

BR 

RP 

0.5 

31 

24 

1.0 

78 

49 

1.5 

92 

74 

1.8  (stoichiometric) 

96 

89 

2.0 

97 

98 

D.  WP/F  Combustion  Studies 

Studieo  were  performed  to  evaluate  the  coabuetlon  of  WP/F  aa  a 
function  of  Initial  oxygen  presaure  and  to  evaluate  the  formation  of 
gaaeoua  products.  WP/F  aamplaa  (~  1  g  each)  were  weighed  into  the 
Monel  reactor  in  a  dry  box.  The  reactor  waa  sealed  and  tranferred  to  a 
vacuum  line.  After  removal  of  the  N2  Atmosphere  and  cooling  of  the 
reactor  to  liquid  nitrogen  temperature  (-195*0),  02  waa  added  from  a 
calibrated  volume.  Reaction  of  WP/7  with  02  generally  occurred  when  the 
reactor  warmed  to  near  room  temperature. 

Figure  5  shows  the  formation  of  CO  and  C02  from  the  combustion  of 
1  g  of  WP/F  as  a  function  of  the  initial  02  pressure.  CO,  present  at 
about  two  orders  of  magnitude  less  than  C02,  reaches  a  maximum 
concentration  near  2.5  atm  02,  whereas  C02  increases  with  pressure.  As 
shown  in  the  Figure,  conversion  of  felt  to  C02  is  far  less  than  the 
theoretical  value  of  39.5  x  10“ 3  moles. 

A  comparison  of  the  conversion  of  butyl  rubber  in  RP/BR  and  felt  in 
WP/F  to  C02  is  shown  in  Figure  6.  The  results  demonstrate  a  significant 
difference  in  the  burning  characteristics  of  the  two  munitions:  at  the 
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FIGURE  5  FORMATION  OF  CO  AND  COj  FROM  THE  COMBUSTION  OF  1-g 
OF  WP/F  AS  A  FUNCTION  OF  02  CONCENTRATION 
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MOLES  CO  x  105/ g  SAMPLE 


FIGURE  6  PERCENT  COMBUSTION  OF  BUTYL  RUBBER  (BR)  AND  FELT  IF) 


AS  A  FUNCTION  OF  02  CONCENTRATION  FOR  1-«  OF  RP/BR 
AND  WP/F 
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COMBUSTION  OF  FELT 


stoichiometric  Oj  pressure  for  RP/BR,  98%  of  the  butyl  rubber  is  burned 
but  only  11%  of  the  felt  is  converted  to  C(>2«  We  estimate  for  WP/F  that 
even  in  the  presence  of  50%  excess  02,  only  about  25%  of  the  felt  Hill 
burn  completely. 

The  degree  of  oxidation  of  phosphorus  to  P4O1Q  (assumed  product)  is 
similar  for  both  munitions.  As  shown  in  Figure  7,  about  90%  of  the  P4 
in  both  munitions  is  converted  to  P4O1Q  at  the  respective  stoichiometric 
O2  pressure.  At  higher  O2  pressures,  however,  a  major  difference  is 
observed.  The  results  suggest  that  nearly  complete  oxidation  of  red 
phosphorus  can  be  achieved  by  Increasing  the  oxygen  pressure,  whereas 
the  upper  limit  of  conversion  of  white  phosphorus  is  about  92%.  Thus, 
in  the  burning  of  WP/F  in  the  environment,  a  significant  amount  of 
unraactad  P4  is  expected  to  remain  In  the  burned  felt  matrix. 

E.  Aerosol  Components 

The  aerosol  components  of  RP/BR  and  WP/F  were  characterised  by  ion- 
exchange  chromatography  followed  by  post-column  hydrolysis  and 
derlvatization  and  detection  at  660  nm.  This  procedure  is  similar  to 
that  described  by  Yamaguchi  (1979).  The  fol1 owing  conditions  were 
employed : 

Instrument :  Spectra- Physics  3500B  Liquid  Chromatograph 

Column:  Bio-Rad  150  x  10  mm  Aainex  A-14  ion-exchange  column 

Mobile  Phase:  A)  0.22  M  NaCl  in  5  oM  EDTA  (disodium) 

B)  0.55  M  NaCl  in  5  mM  EDTA  (disodium) 

Solvent  Program:  0-100%  B  in  A  over  a  6-hr  convex  gradient 

Flow:  1 .0  ml/min 

Detection:  Molybdate  blue  at  660  nm. 

Sampled  of  aerosol  were  collected  by  freeze-trapping  on  silanized 
glass  wool  at  -20*C.  The  aerosols  were  washed  from  the  traps  (two  in 
series)  and  2.0  ml  was  injected  into  the  HPLC  using  a  Water's  WISP 
autoinjector . 
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MOLES  02  x  102 


FIGURE  7  PERCENT  CONVERSION  OF  RP/BR  AND  WP/F  TO  P4O,0 
AS  A  FUNCTION  OF  OXYGEN  CONCENTRATION 
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WP/F  and  RP/BR  gave  nearly  identical  profiles  of  phosphorus- 
containing  species*  identified  as  orthophosphorie  acid  (H3PO4), 
pyrophosphoric  acid  (H4P2O7),  triphosphoric  acid  (H5P301q),  and 
tetrapolyphosphoric  acid  (HgP^O^)  by  comparing  retention  times  to 
authentic  standards  (Sigma  Chemical  Co.,  St.  Louis,  Mo),  plus  a 
homologous  series  of  linear,  condensed  polyphosphates  consisting  of  Pg- 
P^g  polyphosphates  tentatively  identified  on  the  basis  of  their 
chromatographic  properties  and  their  formation  of  the  heteroblue- 
molybdate  complex  (Figure  8).  On  standing  in  water,  the  trl-,  tetra-, 
and  hexa-cyclic  metapolyphosphates  were  formed,  with  the  cyclic 
trlmetapolyphosphate  predominating. 

The  only  lower  oxidation  state  phosphorus  compounds  identified 
besides  phosphine  were  hypophosphorous  acid  (H3po2)  end  phosphorous  acid 
(H3PO3).  These  compounds  were  identified  by  indirect  photometric 
chromatography  (Small,  1982)  as  compared  to  authentic  standards 
(Figure  9).  The  structures  of  these  components  are  shown  in 
Figure  10.  The  products  shown  in  Figure  10  and  WP/F  and  RP/BR  were 
therefore  studied  with  raapect  to  their  persistence  in  selected 
environments. 


COUNTS 


(b) 


LA-4937-10 

FIGURE  9  HPLC  PROFILES  OF  PHOSPHOROUS  SMOKE 

RESIDUE  (BURN  A)  AND  SMOKE  CONDENSATE 
(BURN  0)  FOR  HYPOPHOSPHITE,  PHOSPHATE, 
AND  PHOSPHINE 


WP/F  10] 


FIGURE  1G  COMBUSTION  PRODUCTS  FROM  W/F  AND  RP/BR  BURNS 


IV  PHYSICAL  PROPERTIES  AND  TRANSPORT  STUDIES 


Research  on  physical  and  transport  properties  focused  on  white 
phosphorus  and  phosphine .  The  physical  and  transport  properties  of  red 
phosphorus  are  briefly  reviewed  following  discussion  of  this 
experimental  work. 

A.  Aqueous  Solubility  of  Red  Phosphorus  and  P^ 

The  aqueous  solubility  of  red  phosphorus  was  considered  negligible 
due  to  its  polymeric  nature  and  its  lack  of  discreet  unoxidlzed  units 
detectable  in  the  aqueous  phase. 

The  aqueous  solubility  of  P^  was  measured  at  15,  25,  and  35#C  by 
placing  excess  F4  and  daoxygenated  water  in  a  4-al  vial  equipped  with  a 
teflon  Mininert  valve  (no  headspace).  The  sample  was  warmed  to  45  "C, 
sonicated  for  30  min,  then  stirred  for  48  hr  in  a  temperature-controlled 
bath.  The  vial  containing  the  sample  was  than  centrifuged  at  15,  25,  or 
35*C  and  7000  rpm  for  60  min-.  A  0.5-ml  aliquot  was  withdrawn  with  a 
1-ml  gastlght  syringe  and  immediately  added  to  an  equal  volume  of 
hexane.  The  hexane  extract  was  analyzed  by  GC  using  an  alkali-flame 
ionization  detector.  The  GC  conditions  were: 

Column:  Porapak  P,  glass  (4  ft  x  2 tom  ID) 

Column  temp:  200 *C 

Carrier  gas:  He,  40  ml/mln 

Retention  time.  P4*.  5.0  min. 


Calibration  standards  were  prepared  in  benzene  with  vacuum-dried  P^. 

The  stock  benzene  solution  of  P4  was  prepared  in  a  ^"filled  dry  box  and 
diluted  for  calibration  standards.  The  calibration  was  linear  for 
concentrations  less  than  10  ppm.  Hexane  or  benzene  solutions  of  P4  ware 
stable  within  detector  sensitivity  for  several  days.  The  calibration 
was  checked  daily. 


Table  3  summarizes  che  solubility  data. 


Table  3 

SOLUBILITY  OF  P4  AS  A  FUNCTION  OF  TEMPERATURE 


Tamp 

s  * 

Relative 

Standard 

Deviation 

(*C) 

(ppm) 

<*> 

15 

2.4 

10 

25 

4.1 

4.8 

35 

5.0 

7.6 

•Mean  value  of  four  determinations. 


These  values  agree  reasonably  mall  with  the  solubility  of  3  ppm  at  13*C 
reported  by  Berkovits  at  al.  (1931).  Figure  11  gives  a  plot  of  the 
aqueous  solubility  of  P4  versus  the  reciprocal  of  the  absolute 
temperature.  The  enthalpy  of  solution  (AHg)  of  P4  In  water  can  be 
calculated  from  Equation  2. 

d(lnSJ/d(l/T)  -  -AH  /R,  (2) 

w  s 

where  Sw  Is  the  aqueous  solubility  and  R  is  the  ideal  gas  constant.  The 
calculated  value  of  AH^  is  6  Kcal/mol  from  the  best  fit  of  data  through 
the  three  data  points  of  Figure  11. 

B.  Soil  Sorption  and  Stability  of  P4  in  Soil  Water 
1.  Soil  Characterization 

A  sandy  loam  (supplied  by  Woods  Science,  Monterey,  CA),  PBA  soil,  and 
a  Menlo  Park  soil  collected  by  SRI  staff  were  used  as  test  soils  for  white 
phosphorus  sorption.  Characterization  Included  measurements  of  organic 
carbon  content  (TOC),  particle  size  distribution,  oxalate-extractable  iron, 
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1/T  x  103 

FIGURE  11  AQUEOUS  P4  SOLUBILITY  AS  A  FUNCTION  OF  TEMPERATURE 
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pH,  and  porosity.  The  characterization  methods  used  and  the  results 
obtained  are  summarized  in  Table  4. 


2.  Soil  Sorption  Partition  Coefficient  (K  ) 


Studies  were  conducted  to  measure  the  soil  partition  coefficients 
for  P^  using  soils  both  high  (4.2%)  and  low  (0.3%>  in  organic  content. 
These  studies  were  performed  Inside  a  glove  bag  under  nitrogen  to 
minimize  oxidation  of  the  P^.  Degassed  soils  were  weighed  into  tared 
25-ml  centrifuge  tubes  and  filled  with  aqueous  P^  solution  (no 
headspace).  Test  samples  ranging  from  0.03  to  1.0  g  for  the  high 
organic  content  soils  (Menlo  Park)  and  from  3  to  15  g  for  the  low 
organic  content  soils  (Sandy)  were  used  in  the  batch  experiments. 
Approximately  30  g  of  aqueous  P^  solution  at  1  to  2  ppm  was  used  in  both 
sets  of  experiments.  The  tubes  were  weighed,  wrapped  in  aluminum  foil, 
placed  in  a  rotating-arm  mixer  overnight,  and  then  centrifuged  at  10,000 
rpm  for  30  minutes.  The  supernatant  was  transferred  to  a  septum-sealed 
vial  in  the  glove  bag.  The  soil  was  washed  with  clean  deoxygenated 


water  and  the  residue  was  extracted  with  hexane. 


The  aqueous  supernatant  and  hexane  extract  were  analyzed  by  GC  with 
a  nitrogen-phosphorus  detector.  Our  first  experiments  indicated  a  K 


value  of  20  mi  g~L  and  a  K_„  (HR-VOC  where  OC  is  the  organic  carbon 

oc  p 

content)  value  of  480  mi  g-*1  for  the  Menlo  Park  soil  and  Kp  -  1,  Koc  - 


360  for  the  sandy  loam.  The  experiment  with  the  high  organic  carbon 
Menlo  Park  soil  was  repeated  because  of  poor  mass  balance  recoveries  by 
performing  analyses  of  hexane  extracts  of  the  aqueous  supernatant 
instead  of  direct  aqueous  injection  to  minimize  loss  of  P^  during  this 
step.  Because  of  unaccountable  losses  of  P^  (discussed  below),  the 


upper  K_  limit  that  could  be  calculated  for  the  Menlo  Park  soil  was  25 
“1  ^ 

ml  g  .  This  estimate  assumes  that  all  P^  not  recovered  from  the 


supernatant  is  sorbed.  The  true  value  is  thus  equal  to  or  lower  than 


this  value. 


Our  attempts  to  obtain  reliable  data  for  the  soil  partition 
coefficient  for  white  phosphorus  were  confounded  by  significant  losses 
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of  P4  in  the  "no  soil"  controls.  This  uncertainty  was  too  great  to 
permit  a  reasonable  attempt  to  determine  whether  the  sorption  isotherms 
were  linear  or  whether  sorption  of  P4  significantly  correlated  with  the 
organic  carbon  content  of  the  soil,  or  to  determine  the  leachabllity  of 
P4  in  soil.  Instead,  efforts  were  directed  toward  resolving  the  problem 
of  loss  in  the  "no  soil"  controls. 

A  stability  experiment  was  performed  with  white  phosphorus  in  water 
under  the  best  anaerobic  conditions  achievable  in  a  glove  bag.  The  bag 
Itself  was  flushed  with  nitrogen  (containing  0.1  ppm  O2  after  passage 
through  an  oxytrap). 

P4  solutions  were  prepared  by  allowing  water  to  stand  over  P4  lumps 
for  several  days  in  a  glove  bag.  The  P4  solution  was  then  poured  into 
four  5-ml  glass  vials  (no  headspace).  The  vials  were  sealed  tightly 
with  Mlninert-valve  cape,  and  the  cap-glass  interface  was  wrapped  with 
Teflon  tape.  The  vials  remained  in  the  glove  bag  at  all  times. 

A  fresh  vial  was  sampled  at  intervals  by  removing  P4  solution 
through  the  valve  and  extracting  this  aliquot  with  an  equal  volume  of 
hexane.  The  hexane  extract  was  analyzed  by  GC/NPD.  The  results  are 
presented  in  Table  5. 

Table  5 

OBSERVED  LOSS  OF  P4  FROM  "NO  SOIL"  CONTROL 


an  adsorption  isotherm  experiment .  Because  these  losses  occur 
simultaneously  with  adsorption,  the  results  of  the  sorption  experiments 
are  very  uncertain  and  thus  difficult  to  interpret. 

3.  Stability  of  Dissolved  P4  in  Water 

It  was  observed  in  the  sorption  measurements  that  the  P4 
concentration  in  the  Mno  soil"  controls  was  significantly  lower  than 
that  observed  in  the  supernatant  above  the  "low  soil"  controls.  In  one 
adsorption  experiment  performed  over  16  hr,  an  initial  concentration 
of  1.7  ppm  was  reduced  to  0.38  ppm  in  the  Mno  soil"  controls,  whereas 
concentrations  of  1.5  ppm  and  0.92  ppm  were  observed  in  “low  soil"  and 
“high  soil"  adsorption  tube  supernatants,  respectively,  under  identical 
conditions.  These  data  suggest  that  squaous-extrpctable  soil  matter 
stabilises  P4  in  solution.  Therefore,  an  experiment  was  performed  to 
determine  whether  the  increased  stability  of  white  phosphorus  observed 
in  the  supernatant  of  soil  adsorption  studies  is  due  to  the  aqueous- 
extractable  soil  matter. 

A  P4  solution  was  prepared  by  allowing  water  to  stand  over  P4  lumps 
for  several  days  in  a  glove  bag  filled  with  argon.  The  suspension  was 
warmed  and  then  sonicated  to  hasten  dissolution  of  the  P4.  The  solution 
was  filtered  and  an  aliquot  was  diluted  with  deoxygenated  water 
containing  aqueous-extractable  matter  from  Menlo  Parle  soil  (TOC: 

4.2Z).  This  soil  extract  was  prepared  by  centrifuging  and  filtering  a 
10-wtZ  suspension  of  the  soil.  Ar.  aliquot  of  P4  solution  diluted  with 
deoxygenated  Mllll-Q  water  was  used  as  a  control.  A  series  of  ,2-ml 
vials  were  filled  to  overflowing  with  the  diluted  P^  solutions.  All 
solution  manipulations  wore  carried  out  in  the  glove  bag  in  which  the 
vials  were  stored.  Aliquots  (0.5  ml)  were  periodically  removed  and 
extracted  with  0.5  ml  of  hexane  for  analysis  by  GC. 

The  experimental  results,  shown  in  Figure  12,  indicate  that 
aqueous-extractable  soil  matter  stabilizes  P4  dissolved  in  deoxygenated 
water.  After  24  hr,  70X  of  the  dissolved  P4  remained  in  the 
distilled-water  controls,  whereas  !00Z  remained  in  the  solutions  diluted 


PHOSPHORUS  REMAINING  (%) 


FIGURE  12  STABILITY  OF  P4  IN  AQUEOUS  SOIL  EXTRACT  AND  CONTROL 
WATER  SOLUTIONS 
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with  water  containing  soil  extract.  After  216  hr,  41%  remained  in  the 
controls,  but  86%  remained  in  the  solutions  containing  soil  extract. 
However,  after  100  hr,  both  solutions  appeared  to  lose  P^  at  the  same 
rate. 

C.  Octanol/Water  Partition  Coefficient  (Kow)  of  P^ 

The  Kow  for  P^  was  determined  by  the  method  of  Leo  et  al.  (1971). 
Under  anaerobic  conditions,  a  stock  solution  of  P^  in  water-saturated 
octanol  was  prepared.  Small  volumes  of  this  solution  were  equilibrated 
with  various  volumes  of  octanol-saturated  water.  After  equilibration  at 
ambient  temperature,  the  phases  ware  separated  by  centrifugation  at  1000 
rpm  for  10  min.  The  concentration  of  sample  in  the  octanol  or  aqueous 
phase  was  determined  by  GC.  The  GC  conditions  were: 

Detector:  Alkali-flame  ionization/HP5880A 

Column:  3%  SP2100/Supelcoport,  2'  x  2  mm  10  glass 

Temp  program: 

Temp  1  33*C  (1  min) 

Rate  3  deg  min'1 

Temp  2  60* C 

Carrier:  Helium,  50  ml  min”* 

P4  retention  time:  1.53  min 

The  P^  concentration  was  calculated  using  the  external  standard 
method.  The  results  are  summarized  in  Table  6. 
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Table  6 


Kov 

DETERMINATIONS 

FOR  WHITE 

PHOSPHORUS 

(fa) 

*4 

P4 

Volume  Ratio 

(water) 

(octanol) 

Bun  No. 

(o/w) 

(ng) 

(*g) 

1000 

1 

1/  7.5 

0.27 

281 

2 

1/10 

0.37 

484 

1300 

3 

1/10 

0.25 

318 

1300 

4 

1/10 

0.32 

394 

1200 

5 

1/10 

0.29 

398 

1400 

6 

1/20 

0.55 

669 

1200 

7 

1/20 

0.47 

503 

1100 

Av.  1200  ±  100 

It  la  interesting  to  note  that  several  researchers — for  example, 
Kenaga  and  Caring  (1978)  and  Yalkowsky  and  Valvanl  (1980) — have  reported 
good  correlations  between  Kw  and  Sy  for  a  variety  of  organic 
compounds.  The  values  of  Sy  predicted  by  substituting  KQW  -  1200  into 
these  equations  are  between  one  and  two  orders  of  magnitude  higher  than 
those  we  found  experimentally.  This  is  not  surprising  considering  that 
the  Koy-Sw  correlation  is  based  on  the  assumption  that  octanol /solute 
interactions  are  similar  in  magnitude  to  solute/solute  interactions. 

Our  data  indicate  that,  as  expected,  P4/P4  interactions  are  much 
stronger  than  octanol/P^  interactions. 

D.  Vapor  Presaure  of  P4 

The  vapor  pressure  of  P4  was  measured  at  25*C  by  the  gas  saturation 
method.  Briefly,  a  measured  volume  of  N2  gas  is  saturated  by  passing  it 
through  a  U-tube  containing  glass  beads  coated  with  elemental  phosphorus 
in  a  constant-temperature  box.  The  vaporized  phosphorus  is  trapped  in 
xylene  impingers  and  analyzed  by  GC.  Assuming  Dalton's  law  of  partial 
pressures  and  the  ideal-gas  law,  the  partial  pressure  of  the.  vapor  is 
calculated  by  Equation  3. 
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where  P  a  partial  praaaura  of  phosphorus 
g  •  grass  of  phosphorus 
R  -  gas  constant 
T  «■  taaparatura 
N  -  molacular  weight 

V  ■  total  volume  of  gas  passed  through  the  lmplnger. 
Figure  13  Is  a  diagram  of  the  apparatus  that  we  used. 


Constant  Temperature  Box 


t  ■  Temperature 
M  ■  Molecular  Weight 

V  ■  Total  Volume  of  Gat  Pined  Through  The  lmplnger 

FIGURE  13  GAS  SATURATION  APPARATUS 
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Th«  observed  vapor  pressures  for  duplicate  runs  at  25  ±  0.5°C  were 
0.024  torr  and  0.021  torr.  The  vapor  pressure  calculated  from  the 
Clausius-Clapeyron  constants  as  reported  In  the  CRC  Handbook  (1973)  is 
0.04  torr  at  25*0.  It  should  be  noted  that  the  white  phosphorus  coating 
on  the  glass  beads  in  the  generator  turned  red  during  these 
experiments.  This  indicates  that  the  was  being  transformed  and  thus 
our  measured  vapor  pressures  are  probably  low. 

E.  Henry's  Constant  and  Volatility  of 

Henry's  constant  was  measured  at  25*0  using  the  method  of  Mackey  et 
al.,  1979).  The  apparatus  (Figure  14)  consisted  of  a  glass  column 
fitted  with  a  gas-dispersion  tube.  Nitrogen  was  passed  through  an 
Oxysorb  trap  and  then  saturated  with  water  before  bubbling  into  the 
Henry's  apparatus.  Nitrogen  flow,  water  volume,  and  chemical  (P^) 
concentration  were  measured  periodically.  A  plot  of  In  concentration  vs 
purge  volume  produced  a  straight  line,  from  which  Henry's  constant  was 
calculated  to  1600  liter-torr/mole.  This  measured  value  is  about  twice 
as  high  as  the  value  of  700  liter-torr/mole  calculated  from  the  ratio  of 
measured  values  of  the  saturated  vapor  pressure  and  the  aquaous 
solubility  of  P^.  This  is  another  indication  that  our  measured  vapor 
pressure  le  too  low.  If  the  literature  value  of  the  vapor  pressure  is 
used  with  our  measured  value  of  the  aqueous  solubility,  then  the 
calculated  value  of  H  is  very  close  to  the  value  of  H  measured 
directly.  The  aaeumption  that  H  ■  1600  torr  M"1  suggests  that  P^  is  a 
moderately  volatile  chemical  and  that  the  volatilization  rate  will  be 
controlled  by  maaa  transport  resistance  in  both  the  gas  and  liquid 
phases.  The  first-order  rate  constant  for  volatilization  is 
conveniently  expressed  as  Equation  4 

\  lrJL  +  rg/Hf 1  ,  (4) 
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FIGURE  14  SAMPLING  SYSTEM  FOR  DETERMINATION 
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where  L  -  depth  of  the  water  body 
■  liquid  phase  resistance 
rg  -  gas  phase  resistance 

H  ■  Henry's  constant  in  dimansionless  units • 

Assuming  that  r^  ~  200  sec  cm"1,  rg  ~  5  sec  cm”1  (moderately  turbulent 
conditions  in  both  the  liquid  and  air  phases)  and  L  (water  depth)  •  1  m, 
we  can  estlmata  tho  volatilisation  lialf-life  of  in  a  water  body  such 
as  a  pond.  Assuming  that  H  ■  1600  torr  M”1  (0.088  in  dimensionless 
form),  Ky  can  be  calculated  from  Equation  4  and  the  volatilisation 
half-life  (tj./2)  *■*  given  by 

■  -g - -  5.b  hours. 

Henry's  constant,  together  with  an  appropriate  value  of  the  soil 
partitioning  coafficlant,  can  also  be  used  to  estimate  the  partitioning 
of  P^  in  a  soil  column. 

If  equilibrium  partitioning  between  the  soil  air  (A),  soil  water 
(w),  and  soil  solids  (s)  phase  is  assumed,  then  Bombarger  at  al.  (1982) 
have  derived  simple  expressions  for  the  fraction  of  chemical  X  in  these 
phases,  as  shown  in  Equations  5-7. 


X.  “  fraction  compound 
in  air  phase 


81  +  9w  +  (1  -  eA  -  ew)~^ocroc  ° 


(5) 
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X^  ■  fraction  compound 
in  water  phase 
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X  *  fraction  compound 
h  in  solid  phase 
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If  we  assume  for  simplicity  that  the  soil  volume  fractions  are 
equal  (0.  ■  9  ■  9  -  0.33),  the  soil  specific  gravity  d  ■  2.65  g  cm"'*, 

AW* 

the  fraction  organic  content  PQC  -  0.001,  and  K  ■  400  cm^  g"1  (see 
Section  IV. B),  then 

\A  ~  0.01 

\  ~  0.09 
w 

~  0.90 

This  estimate  indicates  that  dissolved  will  partition  primarily 
between  the  soil  solids  and  vatsr  phases  of  a  wet  soil.  The  fraction  of 
P4  that  partitions  into  tha  soil  water  (\^)  is  sufficiently  high  that 
appreciable  leaching  of  dissolved  P^  should  be  expected.  Moreover  X, 
is  sufficiently  high  that  modest  volatilisation  losses  of  may  also  be 
expacted. 

This  modal  of  course  assumes  that  is  dissolved  in  soil  water  and 
achieves  an  equilibrium  distribution  between  the  three  soil  phases.  The 
P4  in  each  phase  will  therefore  be  subject  to  oxidative  or  hydrolytic 
transformation,  dependant  on  available  oxygen  and  water  contents 
associated  with  each  phase. 

P.  Dissolution  Rate  of  P4 

K  flow  system  (Figure  15)  was  used  to  study  the  dissolution  rate  of 
solid  phosphorus  at  25*C.  Oxygen-saturated  water  or  oxygenated  water 
was  passed  through  a  cell  containing  a  piece  of  P4  with  an  exposed 
surface  area  of  1.25  cm^.  The  eluent  was  collected  in  volumetric 
flasks.  In  the  first  experiments,  aliquots  of  the  aqueous  eluent  were 
reacted  with  excess  molybdate  reagent  to  form  the  heteropoly  blue 
complex.  The  total  phosphorus  concentration  was  determined  by 
monitoring  the  UV  absorbance  at  830  nm.  In  aubaequent  axperimants,  the 
aqueous  sluent  was  discharged  into  hexane.  The  hexane  was  then  analyzed 
for  P4  by  GC  analysis.  Water-soluble  phosphorus  compounds  (presumably 
orthophosphate  was  the  principal  compound)  were  again  reacted  with 
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FIGURE  15  DISSOLUTION  APPARATUS 
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molybdate  reagent  and  total  phosphorus  in  the  aqueous  extract  was 
determined* 


The  initial  measurements  were  done  using  a  Van  Roy  pump  with  a  flow 
rate  of  about  3  ml  min-*.  The  steady-state  total  phosphorus 
concentration  was  0.03  ±  0.004  ppm,  which  corresponds  to  a  dissolution 
rata  of  6.6  tig  hr-1.  Subsequent  runs  were  performed  at  lower 
volumetric  flow  rates  using  both  the  Van  Roy  pump  and  a  Gilson  Mini  2 
pump  for  the  lower  flow  rates.  All  of  these  experiments  were  performed 
using  oxygen-saturated  water.  The  results,  (Table  7)  indicate  that,  In 
general,  the  average  dissolution  rate  decreases  with  decreasing  flow 
rata.  This  Is  expected  because  at  lower  rates  the  hydrodynamic  boundary 
layer  next  to  the  phosphorus  becomes  thicker,  which  retards  the 
diffusion  of  oxygen  and  oxidised  species  to  and  from  the  phosphorus 
surface,  respectively. 


Table  7 

SUMMARY  OF  DISSOLUTION  EXPERIMENTS  IN 
02-SATU RATED  WATER 


Flow  Rata 
(ml/min) 

Avg 

Aqueous 

Cone  P 
(ppm) 

Average 

Dissolution  Rate 
(pg/hr) 

3.00 

0.03 

6.7 

0.54 

0.05 

2.0 

0.42 

0.15 

3.5 

0.34 

0.13 

2.6 

0.14 

0.32 

1.5 

^-purged  water 

0.33 

0.12 

2.4 

0.29 

0.13 

2.5 
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Rate  studies  done  with  nitrogen-purged  water,  fresh  P^  water,  and 
aged  P4  water  gave  the  sane  dissolution  rates  at  steady  state  as  those 
with  oxygen-saturated  water  (compare  the  dissolution  rates  for  the  three 
experiments  with  a  flow  rate  of  approximately  0.3  ml  mln*^).  As  shown 
in  Figure  16,  the  rate  of  orthophosphate  (and  other  water-soluble 
oxidation  products)  dissolution  decreased  rapidly  as  the  oxidized 
species  formed  during  the  loading  of  the  cell  were  dissolved  and  swept 
away  in  the  fresh  P4  water.  The  aged  P^  water,  which  had  been 
conditioned  for  10  days  with  oxygen-saturated  water,  initially  eluted 
more  oxidized  species,  but  the  rate  of  dissolution  decreased  rapidly. 

The  total  amount  of  phosphorus  eluted  with  the  oxide  species  was  about 
25  |ig  for  the  aged  water.  If  the  oxidized  species  were  orthophosphate, 
the  thickness  of  this  oxide  layer  for  the  exposed  area  (1.25  cm2)  can  be 
calculated  to  be: 


4 

V  1 .23  car 


^97  g  H  P0 


or  about  3000  A. 

Thus,  at  long  times,  when  the  oxide  species  are  dissolved,  the 
steady-state  dissolution  rates  are  dominated  by  the  dissolution  of  P^, 
and  therefore  the  fresh  and  aged  waters  gave  the  same  dissolution  rate 
(Figure  17).  Moreover,  at  the  relatively  high  flow  rates  used  in  these 
experiments,  the  dissolved  P^  did  not  oxidize  appreciably  before  being 
collected  in  the  hexane  and  therefore  the  concentration  of  oxide  species 
In  the  aqueous  eluent  became  negligible  at  long  times. 

At  lower  flow  rates,  the  rate  of  dissolution  of  P4  is  expected  to 
diminish  further.  Unfortunately,  0.14  ml  min”*  was  the  lower  limit  of 
volumetric  flow  that  we  could  obtain  with  our  flow  pumps.  However,  the 
dissolution  rate  of  1.5  pg/hr  per  1.25-cra2  exposed  surface  that  we 


51 


FIGURE  17  CUMULATIVE  MASS  OF  PHOSPHORUS  ELUTED  AS  f>4 
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measured  at  this  flow  rata  can  be  used  to  eetimate  an  upper  limit  to  be 
expected  for  dissolution  of  P^« 

For  example,  a  1  co^  spherical  particle  of  P^  In  soil  would  have  an 
exposed  surface  area  of  4.8  cm^.  Thus,  the  amount  of  dissolved  In 
-  one  hour  would  be 

2 

<1.5  pg/hr)  .5.8  pg/hr. 

1.25  am* 


The  amount  dissolved  in  one  year  would  be  51  mg  or  about  2.8%  of  the 
total  mass  of  the  1  cm^  particle.  Thus,  ws  would  expect  that  this  P4 
particle  would  persist  in  the  environment  for  many  years.  The 
persistence  of  this  particle  would,  of  course,  be  orders  of  magnitude 
greater  than  the  persistence  of  dissolved  in  soil  water. 

6.  PH3  Properties-- Solubility  in  Water  and  Henry's  Constants 

PH3  (1010  ppm  in  He)  was  bubbled  through  300  ml  of  Milli-Q  water 
contained  in  the  Henry's  apparatus  (Figure  14)  submerged  in  a  constant-' 
temperature  bath.  Water  samples  were  withdrawn  for  analysis  of  PH3 
content  from  a  sample  port  ~  8"  from  the  tube  bottom  and  ~  4"  beneath 
the  Mllll-Q  water  surface.  Solubility  measurement  were  conducted  at 
three  temperatures  (9.6,  29.7,  and  49.7*C).  Samples  were  taken  until 
the  PH3  analyses  gave  consistent  results,  usually  in  about  60  min. 
Bubbling  was  stopped  just  before  sampling  to  prevent  the  presence  of 
bubbles  in  the  sample.  The  results,  shown  in  Figure  18,  fit  a  good 
straight  line.  From  Equation  2,  the  slope  of  the  line  timns  2.3  R 
equals  the  heat  of  solution  of  PH3  in  water-  The  result,  -5.65  kcal 
mole"*,  suggests  that  PH3  interacts  with  water,  but,  to  a  lesser  degree 
than  does  ammonia  (A^  «  -8.28  for  NHy. 

The  measurement  also  permits  calculation  of  Henry's  constant  (Hc) 
because  at  saturation,  the  partial  pressure  of  PH3  above  the  solution 
equals  the  concentration  of  PH3  in  the  flow  stream.  The  result  at  25°C 
is  a  value  of  6.8  x  10^  L-torr  mole-*-. 
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The  solubility  of  PHj  in  water  has  been  previously  measured  by  A. 
Stock  et  al.  (190$)  and  R.  E.  Weston  (1954) .  Weston  reports  solubility 
data  in  terms  of  the  Ostwald  coefficient  0  (cone •  of  PHj  is 
solutlon/conp.  PH3  is  the  gas  phase).  Values  for  0  are  unitless  because 
the  same  concentration  units  are  used.  Weston's  data  summary  is  poor 
because  he  reports  only  PH3  pressures  and  not  concentration.  However, 
the  concentration  of  dissolved  PH3  (moles/L)  can  be  calculated  by: 

tP«3]g 
tpH3W  "  "IT1 

where  [PH3]g  -  P/RT. 

Henry'a  constant  at  25*0  from  Weston's  results  Is  3419  1-torr  mole™*1,  a 
value  significantly  lower  than  ours.  Furthermore,  from  the  effect  of 
temperature  on  solubility,  Weston  calculates  AHbq1  ■  >2.95  kcal  mola~^, 
which  is  lower  than  our  value  by  a  factor  of  about  two.  The  major 
problem  with  Weston's  results  is  the  uncertainty  of  the  indirect  method 
that  ha  used  to  measure  dissolved  PH3»  Phosphine  was  extracted  from  the 
solution  by  repeated  cycles  of  Pressing,  thawing,  and  stirring,  and  the 
amount  of  extracted  gas  was  measured  in  a  gas  burette. 

Weston  concludes  that  Henry's  law  for  PHj  is  obeyed  from  100  to  700 
mm  and  that  PH3  does  not  show  acidic  or  basic  character  when  dissolved 
in  water.  Van  Waser  1958  reports  that  PH3  in  water  undergoes  gradual 
decomposition  Into  P4,  H2,  and  a  hydride  of  the  approximate  formula 
P2H.  A  solid  hydrate  of  PKj  has  been  reported  (Skinner  1887)  with  the 
probable  formula  PH^»5.9  HgO  (Claussen,  1951). 

The  results,  summarised  in  Tables  8  and  9,  suggest  that  PH3  is  more 
soluble  in  water  and  has  a  higher  heat  of  solution  than  02. 
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Table  8 


HENRY'S  CONSTANT  AND  SOLUBILITY 
FOR  PHOSPHINE.  OXYGEN.  AND  NITROGEN 


Solubility  In  H2O 


Gaa 

*h 

SRI 

(torr-L  rnole”^) 
Weaton* 

Bettinob 

at  1  atm 

(molea/L) 

ph3 

6.80  x  104 

3.62  x  103 

- 

1.12  x  10“2  (SRI) 

®2 

- 

- 

5.99  x  10s 

1.27  x  10"3 

n2 

- 

- 

1.17  x  106 

6.47  x  10‘4 

*Weaton  (1954). 
"Battlno  (1966). 


Gaa 

HC1 

nh3 

°2 

ph3 


Table  9 

HEATS  OF  SOLUTION  OF  GASES 
IN  WATER 

AHioi  (Real  mole  1) 

NBS*  Battinob  SRI 

-17.9 
-8.28 

-3.44 

-5.65 


*NBS  Circular  500. 

*R.  Battlno  (1966). 

H.  Phyaical  Propertlea  and  Tranaport  Studies  Summary 

The  aqueoua  aolubilitiea  for  P^  were  reaaonably  conaiatent  with 
reported  valuea.  The  vapor  preaaure  waa  about  half  that  reported  in  the 
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literature.  The  Henry's  constant  that  we  measured  for  P^  using  a 
gas-purge  method  was  about  twice  what  wa  calculated  from  the  ratio  of 
.  measured  pressure  and  solubility.  This  discrepancy  la  probably  due  to 

an  error  in  our  measured  value  of  the  vapor  pressure.  The  Henry's 
v  constant  value  of  1600  torr  M"1  indicates  that  P^  is  fairly  volatile  in 

water.  The  actual  volatilisation  rate  of  P^  from  water  depends  on  the 
turbulence  of  the  water  body  end  the  air.  For  P^  dissolved  in  a  shallow 
pond,  we  estimate  that  the  volatilization  half-life  is  on  the  order  of 
several  hours  to  several  days. 

The  sorption  of  P^  to  soil  is  complicated  by  solution-phase 
chemistry.  P4  oxidizes  rapidly  in  water,  but  the  process  appears  to  be 
retarded  by  the  presence  of  water-soluble  soil  organic  matter,  p^  is 
not  strongly  sorbed  to  soil  and  is  moderately  voletile.  Therefore,  in 
the  absence  of  chemical  transformation  we  would  expect  P^  to  be 
reasonably  mobile  (leachablllty  and  volatility).  However,  because  the 
oxidation  of  P^  is  very  rapid  relative  to  the  estimated  time  scale  for 
leaching  or  volatilization,  detailed  estimation  of  P4  mobility  is  not 
warranted. 

The  persistence  of  particulate  P^  was  also  estimated  by  measuring 
the  dissolution  rate  of  P4  wafers  in  water  at  25*C.  The  rate  was  found 
to  be  very  low,  therefore,  we  estimate  that  particulate  P4  will  persist 
in  e  terrestrial  environment  for  many  years. 

The  solubility  of  red  phosphorus  in  water  was  too  low  for  us  to 
measure.  Because  of  this  very  low  solubility,  we  expect  that  the  rate 
of  dissolution  of  red  phosphorus  in  water  is  considerably  lower  than 
that  of  white  phosphorus.  The  rate-determining  step  for  dissolution  of 
red  phosphorus  under  aerobic  conditions  is  very  likely  to  be  oxidation 
of  the  outer  layer  of  red  phosphorus  particles. 
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V  CHEMICAL  TRANSFORMATIONS 


The  chemical  transformation  studies  centered  primarily  on  the 
oxidative  transformat lone  of  the  parent  compounds,  RP/BR  and  WP/F,  the 
hydrolysis  of  polyphosphates,  and  the  reactions  of  phosphine  with 
hydroxyl  radical  and  osone.  These  processes  in  soil,  water,  or  air 
environments  are  discussed  below. 

A.  Oxidation  of  Solid  RP/BR  in  Water 

The  potential  exists  for  RP/BR  to  be  deposited  in  aquatic 
environments  from  its  accidental  or  intentional  deployment  into  water 
bodies.  Because  RP/BR  has  essentially  aero  solubility  in  water,  any 
oxidative  tranaformatlons  would  occur  at  the  solid-water  interface, 
eventually  producing  orthophoephorlc  acid  by  hydrolysis  of  phosphorus 
oxides. 

To  study  this  blphaslc  reaction  as  a  function  of  particle  else, 
2.0-g  eamples  of  RP/BF.  pellets  and  powder  were  weighed  into  100-cc 
volumetric  flaaks  and  covered  with  deionised  water.  Synthetic  air  was 
bubbled  into  the  flasks  to  provide  a  continuous  oxygen-saturated 
suspension  at  8  ppm.  Aqueous  safeples  were  withdrawn  periodically  and 
analysed  by  tltrlmetry  for  total  phosphorus  content  (i.e.,  RP  was 
calculated  from  H^PO^).  The  percentages  of  RP/BR  oxidised  as  a  function 
of  time  ere  shown  in  Figure  19.  Orthophosphorlc  acid  was  produced  in  a 
linear  fashion  corresponding  to  sero-order  kinetics  over  the  first 
600  hr  at  1  or  8Z  conversions.  The  projected  lifetime  for  RP/BR  powder 
(<0.01  cm^)  and  pellet  (0.4  cm^)  were  calculated  to  be  1.0  and  8.0 
years,  respectively.  Consequently,  we  project  that  smoke  devices 
deployed  into  aquatic  environments  will  persist,  with  half-lives  greater 
than  4  years. 
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FIGURE  19  OXIDATION  OF  RP/BR  PELLETS  AND  POWER  IN  AIR- 
SATURATED  WATER 
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B.  Oxidation  of  Solid  RP/BR  in  Air  as  a  Function  of  Relative  Humldit; 
and  Oxygen  Partial  Pressure 

RP/BR  can  alao  be  dispersed  on  the  ground  due  to  spillage,  ignition 
aalfunction,  or  incomplete  combustion.  In  these  cases,  air  oxidation 
will  be  a  primary  process  for  transformation. 

To  study  this  process,  2.0-g  samples  of  RP/BR  ware  held  in  a 
desiccator  at  25*C  at  different  relative  humidities  (87Z  and  31%). 
Periodically,  samples  were  removed,  weighed,  and  heated  in  boiling 
water,  and  the  aqueous  phase  was  analysed  by  titrimatry. 

The  pellets  were  found  to  be  very  hydroscopic,  with  weight  gain 
rates  of  0.23  mg  g~*  hr"*  at  87%  relative  humidity  and  0.20  nig  g-1 
hr”*  at  37%  relative  humidity.  This  is  not  unusual  due  to  the  known 
desiccating  properties  of  P2O5  ^A®10^*  *n  expected  oxidative 
tranformatlon  product  (hydrolysis  yields  orthophosphoric  acid). 

Oxidation  of  RP/BR  as  a  function  of  time  appears  to  be  independent 
of  relative  humidity  (Figure  20)  and  follows  zero-order  kinetics  over 
the  first  600  hr.  From  these  data,  a  lifetime  of  3.5  years  waa 
calculated.  At  some  point  we  might  expect  the  oxidation  to  deviate  from 
linearity  and  become  limited  by  oxygen  diffusion  through  the  surface 
oxide  layers.  However,  from  the  length  of  the  projected  lifetime,  we 
would  expect  that  RP/BR  surfaces  would  be  renewed  by  natural 
precipitation. 

The  effect  of  oxygen  pressure  was  examined  to  evaluate  the 
persistence  of  RP/BR  in  environments  where  oxygen  waa  limiting,  such  as 
in  soils  of  various  depths  where  reaction  is  limited  by  diffusion  of 
oxygen  to  the  RP/BR  surface.  For  these  experiments,  2.0-  g  samples  of 
RP/BR  pellets  were  exposed  to  1  atmosphere  of  N2/O2  mixtures  containing 
20,  10,  5,  and  0%  oxygen.  Samples  were  withdrawn  after  300  to  500  hr  of 
exposure  and  analyzed  for  oxidation  products  by  the  standard  titration 
technique.  The  results  are  plotted  in  Figure  21  in  terms  of  the  rate  of 
RP  oxidation  versus  oxygen  partial  pressure.  If  we  assume  that  the 
oxidation  Is  linear  with  time  (at  least  up  to  the  point  where  oxygen 
diffusion  becomes  rate-limiting),  an  initial  half-life  for  each  data 


0  200  400  600 

TIME  (hours) 

FIGURE  20  EFFECT  OF  RELATIVE  HUMIDITY  ON  THE  OXIDATION 
OF  RP/BR  PELLETS 
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point  can  be  calculated  a*  shown  in  Figure  21 .  The  RF/BR  half-life 
versus  oxygen  partial  pressure  is  plotted  in  Figure  22.  This  curve 
allows  us  to  predict  RP/BR  stability  under  partially  anaerobic 
conditions  such  as  in  soils  and  sediments.  It  is  apparent  from  the 
curve  that  RP/BR  will  be  persistent  in  low-oxygen  environments . 

C.  Oxidation  of  Solid  WP/F  in  Water 

The  oxidation  of  solid  WP/F  in  water  was  performed  in  a  similar 
manner  to  that  described  for  RP/BR  in  water  and  was  monitored  by  the 
release  of  orthophosphorlc  acid  to  the  air-saturated  water.  A  plot  of 
the  percent  P4  oxidised  versus  tims  is  shown  in  Figure  23.  By  fitting 
the  data  points  to  a  straight  line  and  assuming  zero-order  kinetics,  a 
half-life  of  2.43  years  (2.35  x  10^X  hr-1)  was  calculated  for  a  1  cm3 
chunk  of  white  phosphorus.  During  the  course  of  this  and  subsequent 
experiments,  the  surface  whitened  with  time  due  to  oxidative  and/or 
hydrolytic  changes.  This  may  limit  oxidation  by  establishing  a  barrier 
at  which  the  tranformation  becomes  limited  by  oxygen  diffusion  through 
the  barrier.  Thla  surface  skin  was  difficult  to  characterize  due  to  its 
poor  solubility  in  water  and  organic  solvents.  These  findings  suggest 
that  the  barrier  is  polymeric  in  nature.  Thus,  we  project  that  the 
half-life  of  2.43  years  is  minimum  for  a  1  cm^  piece  of  white  phosphorus 
and  that  longer  half-lives  could  result  when  diffusion  of  oxygen  becomes 
the  rate-limiting  process. 

D.  Oxidation  of  Solid  WP/F  in  Air  and  Soil 

WP/F  reacts  rapidly  in  air,  beginning  with  wisps  of  white  smoke 
followed  by  inflammation  in  less  than  5  min.  Dainton  and  Bevlngton 
(1945)  have  studied  this  reaction  as  a  function  of  temperature  and  air 
partial  pressure.  Their  data  suggest  that  phosphorus  will  inflame  under 
most  environmental  temperatures  (>5° C)  and  pressures.  Therefore,  no 
further  studies  were  conducted  to  evaluate  this  process  because 
environmental  persistence  is  expected  to  be  extremely  short  for  WP/F  on 
terrestrial  surfaces. 


63 


0 


50  100 

OXYGEN  PARTIAL  PRESSURE 


150 


FIGURE  22  OXIDATION  HALF-LIFE  OF  RP/BR  PELLETS  AT  1  ATMOSPHERE 
AND  100%  RELATIVE  HUMIDITY  AS  A  FUNCTION  OF  OXYGEN 
PARTIAL  PRESSURE 
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From  che  two  examples  of  WP/F  oxidation  in  air  ( 20Z  O2, 
spontaneous)  and  water  (8  x  10"^Z  02,  tj^2  “  2*46  years),  the  dependence 
on  oxygen  concentration  Is  obvious.  Thus,  for  WP/F  burled  In  soil, 
oxygen-diffusion  rates  through  the  soil  column  and  depth  will  be 
critical  factors  in  the  estimation  of  persistence.  Because  most  soils 
become  anaerobic  short  distances  below  the  soil  surface,  we  project  that 
the  half-life  of  WP/F  burled  in  soil  will  be  extremely  long. 

E.  Hydrolysis  of  WP/F  and  RP/BR  Bum  Residues 

The  residues  of  WP/F  and  RP/BR  burns  contain  the  same  profile 
of  linear  condensed  polyphosphates  as  observed  in  the  aerosol  burns 
with  the  addition  of  small  amounts  of  higher  homologs  up  to  P22- 
polyphosphates.  The  pH  of  the  residue  is  low  (2.1)  and  results  in  a 
self-catalyzed  acid  hydrolysis  of  the  polyphosphates.  This  is  readily 
observed  in  the  chromatographic  profiles  of  UP  burn  residues  at  0,  4, 

24,  and  96  hr  (Figure  24).  If  the  rate  of  formation  of  orthophosphorlc 
acid  is  used  as  a  measure  of  the  rate  of  hydrolysis  of  the  mixed  linear 
polyphosphates,  then  a  first-order  rate  constant  of  0.014  hr”1  can  be 
calculated  from  the  0-  and  96-hr  data.  This  corresponds  to  a  residue 
half-life  of  49  hr  for  conversion  to  orthophosphorlc  acid.  In  this  time 
frame,  the  concentrations  of  ?2  and  P3  homologs  increase  due  to  the 
hydrolysis  of  higher  homologs,  but  their  hydrolysis  will  essentially  be 
complete  in  96  hr.  Thus,  the  persistence  of  burn  residues  is  expected 
to  be  short  as  long  as  low  pH  conditions  are  maintained. 


F.  Generation  of  Phosphine  from  WP/F  and  RP/BR 

The  initial  analysis  of  the  vapor  phase  from  WP/F  and  RP/BR  burns 
showed  the  presence  of  small  amounts  of  phosphine  (PHj).  The  finding  of 
a  reduced  product  (oxidation  state  -3)  in  an  oxidizing  environment 
probably  results  from  the  disproportionation  of  lower  phosphorus  oxides 
in  the  presence  of  water,  as  shown  below. 
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FIGURE  24  CHROMATOGRAPHIC  PROFILES  OF  WP  BURN  RESIDUES  AS  A  FUNCTION  OF  TIME 
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To  study  this  effect.,  s  1.67-g  sample  of  WP/F  was  allowed  to  react 
with  1  atm  of  air  at  100%  relative  humidity  in  a  492-cc  Monel  reactor. 
After  the  initial  reaction,  the  pressure  in  the  reactor  was  readjusted 
to  l  atm  with  air.  Vapor  samples  were  removed  periodically  and  analyzed 
for  phosphine  by  GC  the  following  conditions: 

Instrument:  Hewlett-Packard  5730A  Gas  Chromatograph 

Column:  4  mm  x  2.5  m  glass  column  packed  with 
80/100  mesh  Porapak  P 

Column  Temp:  75*C  isothermal 

Flow  rate:  35  cc  min"1  N2 

Detector:  Nitrogen-phosphorus  (NP) 

Instrument  calibration  was  psrformed  using  10-  and  1000-ppm  standards  in 
helium  obtained  from  Matheson  Gas  Products. 

The  formation  of  phosphine  as  a  function  of  time  is  shown  in 
Figure  25.  Nearly  500  ppm  of  phosphine  (0.1%  of  Initial  P^)  was 
generated  after  300  hr.  Although  the  rate  of  oxidation  was  not 
determined,  the  oxidation  of  WP/F  was  91%  complete,  as  measured  by 
titrlmetry,  at  the  end  of  the  experiment. 

The  effect  of  oxygen  on  the  generation  of  phosphine  was  inves¬ 
tigated  by  weighing  1-g  samples  of  WP/F  into  the  Monel  reactor  and 
condensing  various  oxygen  concentrations  into  the  reactor  at  -196aC. 

The  reactor  was  allowed  to  warm  to  room  temperature,  at  which  a  vigorous 
reaction  occurred.  Vapor  samples  were  removed  approximately  l  hour 
after  combustion  and  analyzed  for  phosphine.  The  concentration  of 
phosphine  versus  oxygen  concentration  is  shown  in  Figure  26.  The 
results  demonstrate  that  in  the  presence  of  water,  phosphine  is  formed 
in  a  logarithmic  fashion  that  is  inversely  related  to  the  oxygen 
concentration.  This  is  a  reasonable  expectation  for  the  decomposition 
of  reactive  intermediate  (X)  by  competitive  pathways  (Equation  8). 
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FIGURE  26  GENERATION  OF  PHOSPHINE  FROM  1 ,0-fl  OF  WP/F  AS  A  FUNCTION 
OF  OXYGEN  CONCENTRATION 
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Thus,  phosphine  formation  is  favored  at  low  oxygen  pressures  such  as 
night  be  found  in  soils  or  storage  tanks  when  moisture  is  presents 

The  generation  of  phosphine  from  RP/BR  burns  was  also 
Investigated.  Analogous  to  the  WP/F  study,  2-g  samples  of  RP/BR  pellets 
were  ignited  with  a  Nlchrom  fuse  wire  in  3  atm  of  oxygen  at  100% 
relative  humidity  in  the  492-cc  Monel  reactor.  Following  combustion, 
the  reactor  was  adjusted  to  1  atm  with  air.  Vapor  aliquots  were 
periodically  removed  and  analysed  for  phosphine  by  GC.  After  time 
periods  up  to  200  hrs,  no  phosphine  was  detected  at  30,  50,  or  60% 
conversions  of  RP/BR.  Although  we  had  difficulty  in  maintaining  the 
burn  of  RP/BR  in  the  reactor,  we  conclude  that  the  partial  oxidation  of 
RP/BR  yields  a  lower  fraction  of  lower  phosphorus  oxides  than  was 
observed  with  WP/F.  This  result,  however,  is  inconsistent  with  the 
similarity  in  aerosol  products  between  the  two  smoke  agents. 

In  studies  designed  to  determine  products  arising  from  uncombusted 
RP/BR  (euch  as  materials  that  might  be  discarded  or  partially  combusted) 
we  were  able  to  detect  the  generation  of  phosphine.  In  this  study, 

RP/BR  pellets  were  maintained  in  the  reactor  in  atmospheres  containing 
5,  10,  and  100%  oxygen  at  1  atm  and  100%  relative  humidity.  The 
formation  of  phosphine  at  each  oxygen  concentration  versus  time  is  shown 
in  Figure  27.  In  thia  case,  phosphine  concentration  increases  with  time 
and  oxygen  concentration  (note  Y-axis).  This  is  opposite  to  what  was 
observed  for  WP/F  and  is  a  result  of  the  nonspontaneity  of  RP/BR 
reactions  in  air  (i.a.,  higher  energy  of  activation  for  polymeric  P^  vs 
monomeric  P^).  The  expected  initial  rate  of  phosphine  formation  (up  to 
200  hr)  vs  oxygen  concentration  is  shown  in  Figure  28.  For  normal 
environmental  conditions  under  high  humidity,  we  project  a  rate  of 
approximately  0.3  ppm  hr-*1. 

If  the  generation  of  phosphine  is  dependent  on  water  as  a  reactant, 
its  rate  of  formation  should  be  a  function  of  relative  humidity.  To 
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FIGURE  27  GENERATION  OF  PHOSPHINE  FROM  RP/BR  PELLETS  AT  100%  RELATIVE 
HUMIDITY  AS  A  FUNCTION  OF  TIME  AND  PERCENT  ATMOSPHERIC 
OXYGEN 
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study  this  effect,  2-g  samples  of  RP/BR  were  placed  In  a  glass  beaker  in 
the  Monel  reactor  at  1  atm  air  above  constant  humidity  solutions 
(saturated  CaCl2*6H20  -  32%  RH;  saturated  Ca(N0j)2»4H20  ■  52%  RH; 

H2O  ■  100%  RH) .  Vapor  aliquots  were  removed  periodically  and  analysed 
for  phosphine  by  GO. 

Phosphine  concentration  vs  time  for  each  relative  humidity  are 
plotted  in  Figure  29.  These  data  Indeed  show  a  phosphine  dependence  on 
relative  humidity  for  the  initial  formation  rates.  The  leveling-off 
effect  is  probably  due  to  surface  coatings  that  limit  the  diffusion  of 
oxygen  to  the  RP/BR  surface.  A  semilog  plot  of  initial  rate  vs  relative 
humidity  la  shown  in  Figure  30  and,  except  for  one  data  point,  gives  a 
nearly  linear  fit  over  humidities  ranging  from  32  to  100%. 

The  fraction  of  mass  of  RP/BR  converted  to  phosphine  is  small. 
However,  because  of  the  toxic  nature  of  the  latter  component,  the 
conditions  that  favor  its  formation  should  be  noted,  especially  in 
environments  where  RP/BR  may  be  concentrated. 

G.  Tranformations  of  WP  in  Water— ‘Oxidation  and  Hydrolysis 

The  stability  of  WP  in  water  has  been  reviewed  by  Sullivan  (1979); 
variable  half-lives  ranging  from  0.85  to  700  hr  were  reported.  Our 
review  of  these  studies  revealed  that  the  experimental  design  was 
inadequate  in  many  cases  because  losses  due  to  volatilisation  were  not 
considered.  We  therefore  performed  experiments  to  identify  the 
transformation  processes  responsible  for  the  loss  of  WP  in  water  and  to 
estimate  the  rate  of  each  loss  process.  Losses  due  to  volatilisation 
were  discussed  in  Section  III.E. 

The  loss  of  P^  in  water  was  determined  in  sealed  glass  reaction 
flasks  with  high-purity,  metal-free,  deoxygenated  water  prepared  using 
vacuum-line  techniques.  This  precaution  was  taken  because  of  the 
reported  catalytic  effects  of  motal  ions  on  WP  transformations 
(Sullivan,  1979).  The  prepared  water  was  saturated  with  either  air  or 
argon  for  2  hr  and  filled  the  flask  to  90%  of  its  volume.  WP  in  ethanol 
was  then  Introduced  into  the  aerated  or  argonated  reaction  flasks  by 
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FIGURE  29  EFFECT  OF  RELATIVE  HUMIDITY  ON  THE  GENERATION 
OF  PHOSPHINE  FROM  RP/BR  PELLETS  AT  1  ATMOSPHERE 
AIR 
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FIGURE  30  GENERATION  OF  PHOSPHINE  FROM  RP/BR  PELLETS  IN  AIR 
AS  A  FUNCTION  OF  RELATIVE  HUMIDITY 


77 


syringe  through  a  Teflon-rubber  septum  to  yield  final  WP  concentrations 
of  0.5  to  1.0  ppm.  Upon  the  introduction  of  WP/ethanol,  white  vapor  was 
observed  above  the  solution.  The  vapor  dissipated  slowly— more  rapidly 
with  air  than  with  argon.  An  analysis  of  the  vapor  by  GC  showed  P4 
present  in  the  headspace. 

Aliquots  of  the  aqueous  solutions  were  removed  at  selected  time 
intervals,  extracted  with  hexane,  and  analysed  for  WP  by  GC.  A  first- 
order  plot  of  the  loss  of  P^  in  argon-  and  air-saturated  water  is  shown 
in  Figure  31.  The  data  show  good  correlation  for  a  first-order  kinetic 
plot  and  suggest  that  both  hydrolysis  and  oxidation  occur  in  water.  The 
apparent  first-order  rate  constant  for  hydrolysis  is  Vo"  8’27 
(±  1.21)  x  10"3  hr”1  (tj_/2  ■  84  hr)  and  for  oxidation  p£ua  hydrolysis 
is  k^  Q  +  0  -  1.66  <±  0.47)  x  10“2  hr”1  (tj^2  “  42  hr).  !Jy 
difference, 2 the  oxidation  rate  constant  is  k^g  -  8.33  x  10“3  hr”1, 
nearly  identical  to  the  hydrolysis  rate  constant. 

If  hydrolysis  occurs,  reduced  P4  species  must  be  produced  along 
with  phosphorus  oxides  (H^POj,  HjPOj,  or  H3PO4).  The  only  reduced 
species  that  result  from  phosphorus  hydrolysis  are  hydrogen  and 
phosphine.  To  test  for  the  production  of  hydrogen  gas,  P4  was  placed  in 
water  in  an  N2  atmosphere,  and  headspace  gases  were  analyzed  for 
hydrogen  by  GC  and  a  thermal  conductivity  detector.  After  94  hr,  no 
hydrogen  gas  could  be  detected.  A  similar  experiment  was  performed  to 
detect  phosphine.  The  results  of  these  experiments  are  detailed  in 
Tables  10  (absence  of  oxygen)  and  11  (presence  of  oxygen).  The  rate 
constants  for  loss  of  P4  calculated  from  these  data  are  8.27  (±  0.15)  x 
10“3  hr"1  for  hydrolysis  (t1y2  ■  84  hr)  and  1.23  (±  0.05)  x  10"2  hr”1 
for  oxidation  plus  hydrolysis  (t^2  “  36  hr).  The  data  show  that  a 
small  percentage  of  P4  lost  is  converted  to  phosphine,  which  appears  in 
the  headspace.  Analyses  of  phosphate  (PO4"3),  phosphite  (PO3"2), 
and  hypophosphite  by  indirect  photometric  chromatography  and 
spectrophotometry  (molybdate  complex)  show  erratic  results  during  the  P4 
transformation,  as  shown  in  Table  12.  This  may  be  a  result  of  unknown 
Interfering  ions  in  the  profile  or  poor  peak  integration  of  small  peaks 
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•luting  close  to  ««ch  oth«r.  However,  the  production  of  phosphine  does 
suggest  that  internal  oxidation-reduction  reactions  are  occurring  and 
that  they  nay  be  hydrolytic  in  nature  or  may  involve  oxygen  after  the 
generation  of  a  hydrolytic  intermediate. 

A  possible  hydrolytic  mechanism  that  could  generate  phosphine  is 
shown  in  Scheme  1. 


SCHEME  1 


The  overall  raactlon  la  shown  In  Equation  9. 


P4  +  5H20  - ►P^  +  2H2P0H  +  H3P03  (9) 


Tabic  12 

THE  FORMATION  OF  PHOSPHATE,  PHOSPHITE,  AND 
HYPOPHOSFHITE  IN  P4  TRANSFORMATION  STUDIES 
(PP®) 


Time 

P4  x  iO2 
mmoles 

H3PO2  X  102 
mmoles 

H3PO3  x  102 
mao las 

H3PO4  x  102 
mmoles 

(hr) 

gar 

m. 

Oxid 

Hyd 

Ox  Id 

521 

Oxld 

hj!_ 

0 

2.2 

2.0 

0.9 

0.2 

0.4 

1.0 

1.1 

1.4 

17 

2.1 

1.6 

2.6 

3.1 

0.8 

1.8 

0.3 

3.0 

40 

1.4 

1.3 

1.0 

0.4 

1.2 

1.3 

3.6 

1.3 

161 

0.3 

0.1 

1.8 

6.0 

1.2 

5.6 

2.4 

14.4 

Although  thcac  data  arc  Inconclusive  as  to  the  exact  mechanism  of 
loss  in  water,  they  do  indicate  that  the  loss  does  not  depend,  to  any 
large  degree,  on  oxygen  and  that  half-lives  are  generally  less  than 
100  hr. 

Contrary  to  the  above  transformation  results,  we  found  that  P4  was 
stable  in  water  for  up  to  21  days  in  the  presence  of  cysteine  and  sodium 
thioglycollate,  which  were  used  as  oxygen  scavengers  (see 
Biotransformation,  Section  VI. A).  We  reported  that  P4  solutions  were 
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also  stabilized  in  natural  watara  in  the  absence  of  oxygen  (see  Physical 
Transport,  Section  IV *B).  In  an  attempt  to  understand  the  stabilization 
afforded  by  these  solutions,  we  repeated  the  cysteine/ sodium  thio- 
glycollata  experiment  under  conditions  Identical  to  those  in  the 
hydrolysis  studies,  using  0.25  g  each  of  cysteine  and  sodium 
thloglycollate  dissolved  in  500  ml  of  argonated  Hllllpore  water.  The 
results  of  three  trials  are  shown  in  Table  13. 


Table  13 

P4  HYDROLYSIS  STUDIES  IN  THE  PRESENCE  OF 
CYSTEINE  AND  SODIUM  THIOGLYCOLLATE 


Trial  1 

Trial  2 

Trial 

3 

Time  (hr) 

*4  (PPo) 

Time  (hr) 

P4  (ppm) 

Time  (hr) 

P4  (Ppn) 

0 

0.42 

0 

1.00 

0 

0.29 

16.5 

0.39 

16.5 

0.89 

16.5 

0.30 

89 

0.23 

89 

0.57 

89 

0.14 

The  rate  constants  for  the  loss  of  P4  are  Trial  1  -  5.69  (±  0.1)  x 


10’3  hr’1  (t1/2  "  121  hr);  Trial  2  -  1.12  (±  0.17)  x  10’J  hr~A  (t1/2 


100  hr),  and  Trial  3  -  4.99  ±  (1.3)  x  10’3  hr"1  (t1/2  -  139  hr) 
average  for  these  trials  is  3.93  (±  0.52)  x  10"3  hr’l  (t^/2  ■  120  hr). 
These  results  Indicate  some  degree  of  stability  relative  to  the  absence 
of  cysteine  and  sodium  thloglycollate;  however,  the  stability  is  not  as 
great  as  that  observed  in  the  blotransformatlon  studies. 

The  interaction  of  WP  with  other  solutes  wes  further  investigated 

by  UV  absorption  spectroscopy  and  by  nuclear  magnetic  resonance  (NMR- 

31P)  spectroscopy.  We  measured  the  UV-VIS  spectra  in  10-cm  cells  using 

3  x  10’5  M  phosphorus  solutions  in  water  and  in  benzene.  The  \  and 

max 

molar  absorptlvities  appear  in  Table  14. 


The 
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Table  14 


MOLAR  ABSORPTIVIT1ES  AND  ABSORPTION  MAXIMA 
FOR  WHITE  PHOSPHORUS  IN  WATER  AND  BENZENE 


P4  In  water 


P4  In  benzene 


e 

(mol'*1  ca~^) 


320  480 

335  455 

325  4.88  x  103 

340  3.89  x  103 


The  shift  In  X.  _  is  not  unusual  for  different  solvents,  but  the 
increase  in  the  molar  absorptivitlee  by  a  factor  of  10  suggests  a 
solvent-solute  interaction.  This  type  of  association  nay  ninlc  the 
charge-transfer  complexatlon  observed  with  iodine  and  benzene. 

The  Interaction  of  WP  with  benzene  was  further  studied  by 
31P-NMF.  WP  was  dissolved  in  deuterated  cyclohexane  containing  100  ppm 
of  trlbutylphosphata  (as  an  Internal  standard).  The  WP  resonance  was 
observed  at  -526  ppm.  The  chemical  shift  reported  In  the  literature  for 
WP  is  -460  ppm  in  carbon  disulfide  and  is  stated  to  be  concentration- 
dependent  (Crutchfield  et  al. ,  1967).  Benzene  was  then  added  to  the  WP 
solution  and  the  spectrum  was  rerun.  Chemical  shifts  of  up  to  l  ppm 
were  observed.  From  a  plot  of  chemical  shift  differences  vs  added 
benzene,  the  energy  of  interaction  was  calculated  by  the  method  of 
Lynden-Bell  (1969)  to  be  1  Kcal/mole  for  the  WP-benzene  Interaction  in 
cyclohexane.  Although  this  observed  Interaction  appears  to  be  very 
weak,  It  doett  indicate  that  such  interactions  exist  and  could  possibly 
offer  Nome  degree  of  stabilization  of  WP  in  natural  water. 

Previous  to  these  studies,  the  kinetics  and  mechanism  of  oxidation 
of  WP  in  water  were  believed  to  parallel  vapor-phase  reaction  pathways 
where  the  oxidation  proceeds  via  a  radical  chain  mechanism.  The  results 
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of  this  study  demonstrate  that  the  chemistry  is  not  straightforward 
because  the  WP  loss  in  water  was  approximately  the  same  with  or  without 
oxygen.  For  environmental-fate  analysis,  a  rate  constant  of  1  x  10~2 
hr"^  will  adequately  describe  the  chemical  transformation  processes 
(oxidation  plus  hydrolysis);  however,  it  should  be  remembered  that 
stabilisation  of  P4  under  selected  conditions  was  observed.  Although 
the  loss  of  HP  followed  good  pseudo-first-order  behavior,  the  lack  of  a 
complete  material  belance  between  WP  and  its  products  Indicates  that  the 
loss  processes  are  not  completely  understood. 


H.  Hydrolysis  and  Photolysis  of  Selected  Polyphosphates 

Osterheld  (1972)  reviewed  the  literature  on  hydrolysis  of  poly¬ 
phosphates  (discussed  in  the  Phase  I  report),  and  extensive  literature 
exists  lndlcstlng  that  the  hydrolysis  is  highly  pH-dependent ,  showing 
both  acid  and  base  catalysis.  At  neutral  pH,  the  hydrolysis  of  any 
linear  long-chain  polyphosphate  to  shorter  chains  has  a  half-life  around 
20  days.  The  cyclic  polyphosphates  decompose  at  a  much  slower  rate;  at 
neutral  pH,  their  half-life  is  around  3  years. 

Because  these  hydrolysis  reactions  have  been  so  well  studied,  we 
performed  studies  to  confirm  these  numbers  and  to  determine  whether 
photolysis  could  enhance  these  hydrolysis  rate  constants. 

In  the  Initial  studies,  we  measured  the  UV-VIS  spectra  and  NMR 
spectra  for  orthophosphate,  phosphorus  acid,  hypophosphorus  acid, 
pyrophosphate,  trlpolyphoaphate,  hexametaphosphate,  and  the  hydrolytic 
products  of  phosphorus  pentoxide  in  water. 

The  UV-VIS  spectra  were  measured  over  a  concentration  range  of  10“^ 
to  10“l  M  in  water  in  1-cm  cells.  Over  the  concentration  ranges 
studied,  no  systematic  deviations  from  Beer's  Law  were  observed.  The  pH 
of  each  solution  and  ^P-NMR  of  selected  samples  were  measured  to 
confirm  the  species  present.  For  (NaPOj)g,  the  pH  decreased  from  8.1  to 
7.1  with  increasing  concentration,  Indicating  the  development  of  an 
acidic  species-  Only  one  species  was  detected  for  (NaPO^g  in  water 
by  ^P-NMR.  Phosphorus  (V)  oxide  in  water  has  three  species  detectable 
by  ^1P-NMR;  the  mixture  contains  10X  phosphate,  50%  and  40* 

P^0i3®“.  All  other  compounds  were  as  expected. 
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All  the  spectra  contain  a  strong  UV  absorption  with  a  tail 
extending  toward  the  visible  and  in  most  cases,  a  shoulder  superimposed 
on  the  tail.  The  results  are  tabulated  for  all  compounds  in  Table  15. 
Halmann  and  Platzer  (1965)  previously  studied  the  UV-VIS  spectra  of 
Ma2HP0^,  and  our  results  are  in  good  agreement  with  theirs. 

We  concluded  that  photolysis  screening  studies  for  P^Ojq,  (NaP03)g, 
and  Wa^P^O^Q  vers  necessary.  The  extinction  coefficient  values  for 
Na2HPO^  and  Na^PjOy  are  smaller,  but  photolysis  may  also  be  important 
for  these  compounds  if  the  quantum  yield  is  very  high.  Therefore, 
photolysis  screening  studies  were  conducted  for  all  five  compounds. 

We  performed  both  the  photochemical  screening  study  and  the 
preliminary  hydrolysis  screening  study  by  31P-NMR.  The  chemical  shift 
of  selected  compounds  as  measured  on  our  Jeol  7X90Q  or  reported  in 
Crutchfield  et  al.  (1967)  are  given  in  Table  16.  The  hydrolysis  of 
these  chemicals  was  followed  by  comparing  the  Integrated  intensity  due 
to  phosphoric  acid  to  all  other  resonances. 

Each  compound  was  dissolved  in  a  1:1  mixture  of  H20  and  D20  to 
yield  5  x  10"3  M  solution.  The  pH  of  the  solution  was  adjusted  to  7 
with  phosphoric  acid.  Each  solution  was  transferred  to  two  NMR  tubes, 
sealed,  and  analysed  by  3^P-NMR.  One  tube  of  each  compound  was  placed 
In  a  rack  and  exposed  to  sunlight;  the  second  tube,  to  be  used  as  n 
control,  was  covered  with  aluminum  foil  and  also  placed  in  the  rack. 

The  tubes  were  removed  from  the  rack  and  analysed  by  31p-NMR  every  7 
days  up  to  a  total  of  28  days.  No  differences  between  the  blanks  and 
photolyzed  solutions  were  observed,  suggesting  that  photolysis  will  not 
be  an  important  fate  process  for  these  compounds.  No  change  in  the  NMR 
spectra  of  orthophosphate  or  hexametaphosphate  was  observed.  Small 
changes  were  observed  for  pyrophosphate  and  tripolyphosphate,  but  the 
reaction  had  not  reached  the  first  half-life  after  28  days. 

To  further  qualify  the  hydrolysis  rate  of  pyrophosphate  and 
tripolyphosphate,  we  ran  the  reaction  at  pH  3,  7,  and  12  and  analyzed 
the  reaction  by  HPLC.  Solutions  of  1  x  10"^  M  polyphosphate  were  made 
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MOLAR  ABSORPTIVITIES  e  (L  aol  c*)  OF  PHOSPHORUS 
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Table  16 


CHEMICAL  SHIFTS  OF  SOME  PHOSPHORUS  CONTAINING  COMPOUNDS* 


Compound 

JBES. 

b3po4 

0 

Ve 

-112 

*4* 

-460 

H3PO3 

4  <  5 1  Jj]_p  ■  700  Hs 

ph3 

238,  Jy_p  -  180  Ha 

H3PO2 

12,  J^p  -  535  He 

(NaP03)6c 

-22.5 

(NaP03)4c 

-23 

(NaP03)3c 

-21 

VT* 

-11 

p3°lo"3 

.--11.5,  P--S3.9.  J„p 

■  16.7  He 

p4o13-‘ 

a  -  -11.5,  0  ■  -23.9,  J0p 

-  16.7  He 

J8P 

■  15.9  He 

^Reported  In  Crutchfield  at.  *1.  (1967). 

"Measured  In  CS2'  In  hsxane  we  measured  the  chemical  shift  to  be  -526 
opm. 

cThess  resonances  are  broad  and  virtually  indistinguishable  on  our  NMR. 
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In  pH  3,  7,  and  12  buffers*  These  solutions  were  loaded  into  sealed 
glass  ampoules  such  that  there  were  ~  10  ampoules  for  each  chemical  at 
each  pH.  The  ampoules  were  then  submerged  in  an  oil  bath  regulated  at 
25*C;  at  regular  intervals  they  were  removed  and  analyzed  by  HPLC.  At 
pH  3,  a  pseudo-first-order  rate  constant  of  1.39  x  10"  3  hr-1  was 
measured,  yielding  a  calculated  half-life  of  20  days  for  either 
phosphate. 


We  were  unable  to  obtain  rate  data  for  the  pH  7  and  12  solutions 
because  of  very  slow  transformation  rates.  We  suspect  that  at  the 
higher  pHs,  the  high  negative  charge  associated  with  the  ionized 
polyphosphate  limits  a  base  hydrolysis. 


Although  these  findings  may  be  true  for  pyro-  and  tripoly- 
phosphatns,  Griffith  and  Buxton  (1967)  showed  hydrolysis  rates  to 
approach  a  maximum  at  pH  7,  with  increasing  chain  length.  The  first- 
order  rate  constants  ranged  from  4.6  x  10"®  min"*  (tj.^  "  105  days)  for 
tetrapolyphosphate  to  7.5  x  10"3  min"*  (tj^  ■  $.4  days)  for 
heptapolyphosphate  at  pH  7  and  30°C.  These  investigations  suggested  a 
trend  toward  a  maximum  hydrolysis  rate  constant  of  8.0  x  10”3  nin"*  with 
increasing  chain  length.  However,  Gill  and  Rlaz  (1969)  did  not  find 
this  to  be  true;  they  measured  a  first-order  rate  constant  for  a  54 
phosphorous  atom  polyphosphate  of  3.0  x  10"^  rain"*  (ti/2  m  *6°  d*y«)  at 
pH  7  and  40°C. 


The  above  rate  constants  are  the  result  of  two  competitive 
processes,  one  being  an  end-group  clipping  of  orthophosphorlc  acid  and 
the  other  being  the  splitting  of  a  trimetaphosphate  (cyclic 
tripolyphosphate)  from  the  chain.  The  latter  process  is  probably 
influenced  by  a  coiled  structure  of  the  polyphosphate  in  solution.  The 
data  of  Gill  end  Rlaz  show  that  end-group  clipping  is  three  times  more 
rapid  than  trimetaphosphate  production  at  pH  7,  whereas  the  data  of 
Griffith  and  Buxton  show  the  same  ratio  of  rate  constant  to  vary  from  3 
to  5  to  7  for  penta-,  hexa-,  and  heptapolyphosphates .  Although  there  is 
no  obvious  rationale  for  the  discrepancy  in  these  data,  it  should  be 
pointed  out  that  both  processes  do  occur  and  are  environmentally 


relevant. 
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Associated  with  the  biotransformation  studies  of  polyphosphate 
smoke  aerosols  (P^-P^g),  we  observed  complete  stability  of  all 
components  in  pure  water  at  pH  7  over  a  16-day  period*  However,  in 
sterilized  natural  water,  substantial  losses  of  the  P^-P^g 
polyphosphates  were  observed  within  six  days,  with  concomitant  increases 
in  tripolyphosphate  and  trimetaphosphate,  the  cyclic  analog  of 
tripolyphosphate  (Figure  32).  Thus,  catalysis  by  natural  water 
components  (probably  metal  ions,  i.e.,  Mg+2)  appears  to  promote  a 
tripolyphosphate  chain  cleavage  that  assists  the  hydrolysis  of  long- 
chained  polyphosphates  in  the  environment  at  neutral  pHs. 


1.  Transformations  of  Phosphine  (PH^) 


Because  of  the  toxic  nature  of  PH^,  studies  were  performed  to 
evaluate  its  persistence  in  the  environment.  PH^  i®  a  minor  component 
of  the  reactions  of  in  water  and  soil.  Owing  to  its  higher  vapor 
pressure,  (1  x  10^  torr),  almost  all  PHj  formed  will  volatize  to  the 
atmosphere,  where  oxidation  and  photolysis  are  the  important  loss 


processes. 


The  UV  spectrum  of  PH^  shows  no  absorption  beyond  240  nm  (Calvert 
and  Pitts,  1966),  indicating  that  no  photolysis  will  occur  in 
sunlight.  The  other  reactions  of  possible  Importance  are  direct 
reaction  with  0£  and  oxidation  by  ozone  and  OH  radical. 

We  examined  the  stability  of  PH3  in  the  dark  in  the  presence  of  air 
and  found  no  significant  loss  over  a  period  of  several  days.  The 
stability  of  mixtures  of  PH3  and  NO  in  air  also  was  studied;  at  40  ppm 
of  NO  in  air,  we  found  that  the  rate  constant  for  loss  of  PHg  was  4.44  x 
10"21  cc/molec-s,  which  corresponds  to  half-life  in  the  troposphere  of 
about  3000  yr.  Therefore,  we  Initiated  both  literature  and  laboratory 
studies  to  develop  reliable  values  for  rate  constants  for  oxidation  by 


ozone  and  OH  radical. 


We  found  no  literature  values  for  oxidation  of  PH3  by  ozoue, 
but  did  find  one  report,  by  Lorenz  and  Zellner  (1982),  of  the  flash 
photolysis  oxidation  of  PH3  with  OH.  They  reported  a  value  of 
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FIGURE  32  HPLC  PROFILES  OF  CONTROL  WATER  FOR  THE  BIOTRANSFORMATION 


OF  POLYPHOSPHATES 


k(OH)  *  1.6  x  10""11  cc/molec-s,  which  corresponds  to  a  half-life  in  the 
troposphere  of  5  hr,  based  on  an  averaged  OH  concentration  of  2  x  10® 
molec/cc.  For  comparison,  toluene  has  a  value  of  k(0H)  of  5.1  x  10“^ 
cc/molec-s  (Atkinson,  1975). 

Although  k(0H)  values  for  organics  often  are  correlated  with  bond 
strengths  within  families  of  compounds  with  reported  bond  strengths,  no 
such  data  are  available  for  the  phosphine  family  and  we  found  only  one 
reference  to  the  P-H  bond  strength  in  PH3,  reported  by  Sanderson  (1971) 
as  between  70  and  80  kcal/mol.  This  low  value  implies  that  PH3  should 
react  readily  with  a  variety  of  oxidants  in  the  atmosphere,  including 
OH,  NO3,  and  O3. 

1.  Oxidation  with  Ozone 

Several  experiments  were  performed  with  mixtures  of  ozone  and  PH3 
in  air  in  the  dark  at  25*C,  using  1  to  2  ppm  PH3  and  11  to  4000  ppm  O3 
to  measure  the  rate  constant.  These  experiments,  summarized  in  Table 
17,  show  that  the  rate  constant  is  about  0. 5-1.0  x  10”^  cc/molec-s. 
Experiment  4  showed  good  linearity  over  three  orders  of  magnitude  in 
concentration,  an  usually  wide  range  for  a  kinetic  experiment.  Ozone 
typically  is  about  1  x  10®  times  more  abundant  than  OH  radical  in  the 
troposphere  and,  when  concentration  differences  are  taken  into  account, 
lc( O3 )  would  be  comparable  to  a  k(0H)  value  of  1  x  10“^  cc/molec-s,  very 
close  to  the  reported  value. 


Table  17 


OXIDATION  OF  PH3  BY  OZONE 


Run  # 

Phosphine 

<PPa? 

°3 

i£E2l 

1 

2.0 

0 

2 

1.92 

4120 

3 

1.68 

368 

4 

1.32 

21.6 

5 

1.34 

11.5 

_ Results _ 

No  measurable  reaction  after 
86  min 

1002  PHo  reacted  in  less  than 
2  min 

1002  PHo  reacted  In  less  than 
2  min 

-17  3  -1 

k_  ■  1.0  x  10  cm  molec  tec 
°3 

J  -17  3  -' 

k  *  0.36  x  10  cm  molec  *  sec 


2 1  Oxidation  of  PH3  with  OH 

We  initiated  some  experiments  to  measure  k(0H)  for  PH3  to  compare 
with  the  one  literature  value,  which  was  measured  at  lower  pressure. 
Before  doing  any  experiments  with  OH-  generating  system  (Atkinson  et 
al.,  1981),  we  performed  a  few  control  experiments  with  mixtures  of  PH3 
in  air  with  and  without  NO,  MeONO,  and  CCl^— all  constituents  of  the  0H- 
generating  system.  Dark  controls  showed  no  loss  of  FH3  in  any  of  the 
mixtures.  However,  we  found  rapid  loss  of  PH3  (1  to  102/min)  In  several 
mixtures  when  they  were  exposed  to  afternoon  sunlight.  Results  vers 
quite  reproducible  from  one  day  to  the  next  in  a  given  mixture,  but 
showed  wide  variation  (tenfold)  from  one  mixture  to  another.  Thus, 
removal  of  NO  from  mixtures  in  air  markedly  slowed  the  loss  of  PH3,  but 
photolysis  of  PH<j  in  ultrapure  N3  gave  a  rate  ten  times  larger  than  that 
in  air.  In  helium  the  results  are  opposite  to  those  in  nitrogen.  These 
results  are  summarized  in  Tables  18  and  19. 


Because  none  of  the  components  of  these  PHj  mixtures  should  absorb 
any  significant  solar  photons,  we  have  to  seek  an  explanation  for  these 
rapid  reactions  in  some  other  directions  or  at  least  exclude  other 
reasonable  possibilities. 


Table  18 

PHOTOLYSIS  OF  PH3  IN  AIR  OR  NITROGENa,b 


PH3 

(ppm) 

n2 

(ppm) 

°2 

(ppm) 

%  PH3  Reacted 

%  tfH3/min* 

1.45 

1  x  106 

— 

83c»d 

11 

1.45 

9.7  x  105 

2.2  x  105 

uc,d 

1.4 

*Photolyzed  from  3:09  to  3:19  pm  in  cloudy  weather. 
bPH3  derived  from  1100  ppm  PH3  in  N2*  all  concentrations  in  ppm. 
cData  corrected  to  clear-sky  conditions  using  a  power  mater. 
dDark  controls  gave  nc  change  in  PH3 • 


Table  19 

PHOTOLYSIS  OF  PHj* 

PH3  He  O2  N2 

(ppu)  (ppm)  (ppm)  (ppm) 

11  10%  —  — 

1.08  9.8  x  104  7.0  x  105  1.98  x  10* 


aPH3  derived  from  11  ppm  PH3  in  He. 

^hotolyzed  10  min  in  sunlight  from  11:01  to  11:11  am;  dark  control 
stable. 

cPhotolyzed  13  min  in  sunlight  from  5:03  to  5:17  pm;  dark  control 
stable. 


X 

%  /5PH3/min 

3.4^ 

0.34 

19 .8C 

1.46 

We  have  no  clear  explanation  for  the  photoreactivity  of  these  mixtures; 
however,  some  possibilities  are  discussed  below* 

Surface-promoted  photolysis  might  occur  if  PHj  or  some  other 
species  were  sorbed  to  the  glass,  thus  creating  a  surface-sorbed  species 
that  absorbs  light  of  >  300  nm.  To  test  for  this  possibility,  and  using 
carefully  clean  bulbs  (5%  HF,  distilled  water,  baking  at  550°C),  we 
measured  rates  of  loss  of  PH3  in  air  where  the  light  source,  a  filtered 
Hg  light,  either  was  collimated  to  a  narrow  spot  on  the  bulb  or  remained 
uncollimated  to  expose  the  entire  bulb  surface.  No  significant 
differences  in  rates  were  noted,  indicating  that  photolysis  is  a 
homogeneous  process. 

Trace  photocatalysts  in  the  PH3  are  likely  sources  of  the 
phenomenon.  Pluck  (1973)  describes  diphosphine  (P2H^)  as  a  common 
contaminant  in  PH3  that  causes  the  PH^-air  mixture  to  become 
spontaneously  flammable.  Possibly  P2H^  also  is  a  photocatalyst  to 
initiate  chain  oxidation  of  PH3.  However,  P2H^  cannot  explain  the  very 
fast  loss  of  PH3  in  pure  N2«  Only  the  presence  of  another  trace 
oxidant— such  as  N02  or  O3  can— lead  to  loss  of  PH3.  However,  ozone 
will  lead  to  rapid  dark  oxidation  and  N02  needs  to  be  present  at 
stoichiometric  levels  to  effect  significant  conversions  of  PH3.  We 
cannot  rule  out  residual  light  absorption  by  PH3  itself  as  the  cause  of 
the  slower  rates  in  helium  and  air. 

A  practical  approach  to  estimating  the  lifetime  of  PH3  in  the 
troposphere  is  to  expose  PH3  to  sunlight  in  outdoor  air  and  use  the 
measured  rate  constant  to  estimate  the  half-life  or  lifetime.  One 
experiment  with  outdoor  air  with  1.8  ppm  PH3  and  a  15-min  exposure  to 
sunlight  gave  a  rate  of  loss  of  1,6 X  min-1  at  23%  conversion, 

essentially  the  same  as  reported  in  Tables  18  and  19.  If  we  take  0.5% 

min“*  as  the  diurnally  average  rate  of  loss  of  PH3,  this  is  equivalent 
to  a  half-life  in  the  troposphere  of  2.3  hr.  We  can  compare  this  value 

of  t^/2»  with  the  value  of  t^/2  calculated  at  5  hr  for  OH  radical  and  at 

8  hr  for  ozone.  Thus,  ozone  and  OH  radical  have  nearly  the  same 


reactivity  and  both  are  somewhat  slower  than  photolysis.  We  conclude, 
therefore,  that  the  persistence  of  PH3  in  the  environment  will  be  short 
even  Chough  the  roles  of  sunlight  and  other  tropospheric  agent 3  in  PH3 
transformation  are  not  completely  understood. 

J.  Chemical  Transformation  Summary 

The  chemical  fate  of  RP/BR  released  into  the  environment  will  be 
controlled  by  oxidative  transformation.  The  expected  half-life  of  RP/BR 
suspended  in  aerated  aquatic  or  soil  environment  is  expected  to  range 
from  1.7  to  3.0  years.  In  environments  where  transformation  is  limited 
by  oxygen  diffusion  (such  as  in  soil)  or  the  buildup  of  surface  oxide 
layers,  the  half-life  is  expected  to  be  significantly  longer. 

The  chemical  fate  of  WP/F  is  strongly  dependent  upon  the 
environment  in  which  it  is  released.  WP/F  will  oxidize  rapidly  in  air 
at  temperatures  above  5*C.  However,  the  rate  of  oxidation  is  dependent 
on  oxygen  concentrations.  Under  anaerobic  conditions,  such  as  buried  In 
soil,  the  half-life  is  projected  to  be  greater  than  2  years. 

There  is  evidence  to  suggest  that  WP/F  is  lost  by  other  mechanisms 
other  than  oxidation.  Hydrolytic  processes  appear  to  occur  to  soluble 
WP  in  water  at  equivalent  rates  to  that  of  oxidation.  The  half-life  of 
soluble  WP  in  water  is  estimated  at  A2  hours. 

The  products  of  the  chemical  transformation  of  RP/BR  or  WP/F  are 
orthophosphoric  acid,  polyphosphoric  acids,  meta-phoaphoric  acids,  and 
phosphine.  The  linear  polyphosphates  are  rapidly  converted  to 
orthophosphoric  acid  at  pH  3  and  below.  At  environmental  pHs,  their 
stability  is  enhanced  and  half  lives  of  greater  than  100  days  can  be 
expected. 

Phosphine  is  expected  to  partition  to  the  atmosphere  rapidly  where 
reactions  with  ozone  and  hydroxyl  radical  will  limit  it's  half-life  to 
several  hours.  Laboratory  measurements  of  the  rate  constants  for 
oxidation  of  PH3  with  ozone  and  HO  radical  use  much  higher 
concentrations  of  reactants  than  are  found  in  the  atmosphere.  However 


the  calculated  half-lives  for  PHj  in  the  troposphere  are  based  on  space 
and  time-averaged  concentrations  of  ozone  and  HO  in  the  lower 
troposphere  and  are  independent  of  the  method  used  for  measuring  the 
rate  constant. 
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VI  BIOTRANS FORMATION 


The  chemical  transformation  studies  indicated  that  phosphorus 
compounds  would  be  stable  in  environments  low  in  oxygen  and  at  neutral 
pH.  Under  these  conditions,  microbial  populations  offer  an  alternative 
pathway  to  chemical  transformation.  Therefore,  the  biotransformation  of 
several  phosphorous  species  was  investigated. 

A.  WP/F  -  Water 

Because  of  the  rapid  volatility  of  VP  and  its  reactions  with  oxygen 
in  water  solutions,  the  biotransformation  studies  were  conducted  in 
sealed  containers  under  anaerobic  conditions  to  enhance  stability. 

In  the  initial  studies  we  evaluated  the  toxic  effect  of  V?  on 
anaerobic  microorganisms. 

We  collected  soil  at  SRI  and  Inoculated  the  supernatant  from  a 
soil-water  mixture  into  Dlfco  f luid-thioglycolate  medium  or  into  a 
medium  containing  glucose  (500  ppm),  yeast  extract  (100  ppm),  and  basal 
salts.*  The  flasks  were  flushed  with  N2  gas,  stoppered  with  a  Bunsen 
stopper,  and  incubated  at  25aC. 

The  grown  anaerobic  microorganisms  were  inoculated  into  flasks 
containing  medium  and  white  phosphorus  that  had  been  added  to  benzene  to 
give  a  final  concentration  of  1-3  ppm  (near  saturation).  A  control 
flask  was  inoculated  with  solvent  only.  The  turbidity  of  the  madia  was 
measured  after  24  hr;  no  growth  inhibition  was  noted  in  the  white 


ft 

The  phosphate-buffer  basal-salts  medium  contained,  per  liter:  1.8  g  of 
K2HP04,  0.2  g  of  KH2P04,  0.5  g  of  (NH/)2S0a,  0.05  g  of  NaCl,  0.05  of 
MgSO^.^I^O,  0.1  g  of  CaCl2.2H20,  0.003  g  of  FeS0^.7H20,  and  1  ml  of 
trace  elements  solution.  The  trace  elements  solution  contained,  per 
liter:  0.1  g  of  H3BO3  and  0.05  g  each  of  CuS0^.5H20,  MnS0^.H20, 

ZnSO^ .7H20,  Na2Mo0^ ,  and  CoCl2»6H20. 
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phosphorus-containing  solutions  compared  with  the  controls,  indicating 
that  white  phosphorus  will  not  be  tcxic  to  anaerobic  microorganism 
populations  at  aqueous  saturation  levels. 

The  biodegradability  of  WP  was  screened  in  soils  and  sediments 
collected  from  PBA.  For  the  anaerobic  hiotrr ic*rrmation  test,  500-ml 
screw-capped  erlenmeyer  flasks  were  filled  1th  PBA  creek  water  (WP 
Creek  leading  to  Yellow  Lake)  and  0.05  M  tris  buffer  was  added.  The 
flasks  were  placed  in  an  anaerobic  chamber,  purged  with  N2  gas  for  2  hr, 
and  preincubated  for  4  days.  To  one  flask  was  added  50  ppm  of  yeast 
extract  and  enough  WP  from  a  600-ppm  benzene-ethanol  stock  solution  to 
yield  a  final  WP  concentration  of  3  ppm.  Autoclaved  creek  creek 
containing  solvent  only  was  used  as  a  sterile  control.  Our  original 
attempts  to  maintain  sterility  of  autoclaved  cfeek  water  with  mercuric 
chloride  solutions  failed  when  WP  was  observed  to  transform  in  less  than 
one  day,  probably  due  to  a  catalytic  effect  of  mercuric  ion.  Also, 
significant  losses  of  WP  occurred  if  rubber-stoppered  flasks  were  used, 
probably  due  to  sorption  from  the  headspace. 

All  flasks  were  capped  with  Teflon  liners  and  Incubated  in  the 
anaerobic  chamber  in  the  dark.  Aliquots  were  removed  periodically, 
extracted  with  hexane  in  the  anaerobic  chamber,  and  analyzed  for  WP  by 
OC.  Under  the  above  conditions,  WP  was  not  biologically  or  chemically 
tranformed  after  63  days  of  incubation  with  or  witnout  50  ppm  of  added 
yeast  extract.  WP  added  to  a  local  pond  water  was  found  to  be  stable 
after  14  days  of  incubation  under  anaerobic  conditions. 

That  WP  was  stable  in  water  for  33  days  and  in  pond  water  for  14 
days  in  the  anaerobic  chamber  was  at  variance  to  the  chemical 
transformation  studies  on  VP,  which  showed  both  hydrolytic  and  oxidative 
transformatlono  to  occur,  with  half-lives  of  80  and  42  hr 
respectively.  When  the  study  was  conducted  in  PBA  well  water,  WP 
decreased  from  3.9  to  0.6  ppm  in  53  days.  A  sterile  control  decreased 
from  3.6  to  2.6  ppm  during  the  same  period.  However,  when  the  PBA  well 
water  was  inoculated  into  media  containing  TBSM  (0.05  M  tris  buffer 
based  salts  medium  lacking  phosphate),  thloglycolate,  cysteine,  glucose, 


I 
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and  casamino  acids,  no  transformation  occurred  after  32  days.  This 
difference  in  stability  is  puzzling  but  may  be  related  to  WP 
interactions  with  media  components,  which  are  not  well  understood.  We 
also  noted  enhanced  WP-soll  extract  stability  over  spiked  blank 
controls.  Regardless  of  what  these  interactions  are,  we  conclude  that 
biotransformation  of  WP  will  not  occur  rapidly  in  the  environment. 

B.  WP/F-Soll 

A  series  set  of  experiments  similar  to  those  described  above  was 
performed  using  soil  and  creek-bottom  sediments  collected  from  PBA. 

The  creek-bottom  sediment  was  sieved  through  a  2-mm  screen  and  26.7 
g  of  the  wet  sediment  (20  g  dry  weight)  was  put  into  each  of  several 
screw-capped  test  tubes  (25  x  200  mm).  Nitrogen-purged  water  was  added 
to  give  a  total  water  content  of  10  ml  in  each  tube.  The  tubes  were 
placed  in  a  Forma  Scientific  Anaerobic  System  Chamber.  After  5  days  of 
preincubation,  white  phosphorus  in  benzene/ethanol  (1/4  600  ppm  stock 
solution)  was  added  to  yield  6  ppm  in  dry  soil,  and  the  samples  were 
capped  and  incubated  in  the  anaerobic  chamber. 

Sterile  controls  were  prepared  by  autoclaving  the  sediment  twice  in 
successive  days  followed  by  the  addition  of  white  phosphorus.  Previous 
attempts  to  maintain  the  sterility  of  soil  by  autoclaving  followed  by 
mercuric  chloride  addition  resulted  in  a  90X  loss  of  white  phosphorus 
after  one  day  of  incubation.  Therefore,  mercuric  chloride  was  not  added 
in  the  control. 

Periodically,  duplicate  tubes  were  removed  from  the  anaerobic 
incubator,  extracted  with  hexane  anaerobically,  and  analyzed  the  WP. 

The  concentration  of  white  phosphorus  in  PBA  creek  sediment  gradually 
decreased  from  6.1  to  3.1  ppm  during  44  days  of  incubation.  However, 
the  concentration  in  the  sterile  control  also  decreased — from  5.4  to 
3.7  ppm  during  39  days.  These  data  indicate  that  there  is  no 
significant  evidence  for  biotransformation  in  this  sediment. 
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We  had  also  Initiated  a  biotranaformation  screening  test  in  PBA 
soil  obtained  from  a  white  phosphorus  dump  site.  During  35  days  of 
incubation,  the  white  phosphorus  concentration  decreased  from  5.2  to 
0.2  ppm.  The  concentration  in  the  sterile  control  decreased  at  a  slower 
rate-- from  5.2  to  1.7  ppm  over  the  same  period. 

We  observed  the  same  phenomenon  in  SRI  soil;  the  white  phosphorus 
concentration  decreased  from  5.0  to  0.02  ppm  in  4  days  but  the 
concentration  in  the  sterile  control  decreased  from  5.0  to  3.2  ppm. 

These  results  show  that  (1)  white  phosphorus  transformed  in 
flooded,  autoclaved  soila  under  anaerobic  conditions  and  that  the  rate 
varied,  depending  on  the  soil,  and  (2)  white  phosphorus  transformed 
faster  in  some  nonsterile  soils  than  in  autoclaved,  sterile  controls. 
Although  this  transformation  could  be  biological,  some  change  in  the 
chemical  properties  of  the  soil  due  to  autoclaving  could  be  responsible 
for  the  differences  in  transformation  rates. 

To  determine  whether  these  WP  transformations  in  nonsterile  soil 
are  biological,  small  amounts  of  these  WP-transformed  soil-water 
mixtures  were  Inoculated  into  TBSM  with  thloglycolate  and  cysteine, 

100  ppm  glucose,  100  ppm  casamino  acids,  and  3  ppm  WP,  contained  in 
Teflon-lined,  screw-capped  flasks.  After  32  days  of  incubation  no  WP 
biotranaformation  was  observed. 

Because  the  pH  of  the  dump-site  soil  was  3.2,  we  also  prepared 
citric  acid  buffer  (pH  3.6)  to  replace  the  above  trie-buffered  medium, 
and  the  microorganism  from  the  test  tube  containing  WF-tranformed  dump- 
site  soil  were  inoculated.  No  biotranaformation  was  observed  during  12 
days  of  iff  nation. 

From  this  result  it  appears  that  the  difference  in  WP-transf ormation 
rate  thet  we  observed  in  PBA  and  SRI  soils  bstween  autoclaved  and 
nonautoclaved  soil  ia  due  to  chemical  changes  in  soil,  and  hnnce  the 
transformation  was  nonbiologictl,  or  that  biotranaformation  microorganisms 
ars  difficult  to  grow  under  our  test  conditions. 


C.  RF/BR 


Because  of  the  extremely  low  solubility  of  RP  in  water  and  the  lack 
of  evidence  for  the  biotransformation  of  WP  In  water  and  soli,  no 
screening  studies  for  the  biotransformation  of  RP/BR  were  conducted. 

D.  Condensed  Linear  Polyphosphates 

Because  linear  condensed  phosphates  are  present  in  the  smoke 
product,  screening  tests  for  biotransformation  were  conducted  with  local 
water  and  soil.  Tripolyphosphate  (TPP)  was  used  as  a  representative  of 
this  chemical  group. 

For  aqueous  aerobic  biotransformation  screening,  test  water  from  a 
pond  near  Searsville  Lake  was  collected  and  4  or  20  ppm  of  TPP  (1  and  5 
ppm-P,  respectively)  was  added  to  2  liters  of  pond  water  with  and 
without  100  ppm  glucose  and  20  ppm  yeast  extract  as  supplemental  organic 
nutrients.  The  amount  of  condensed  phosphate  was  analyzed  by  the 
difference  in  orthophosphosphate  content  before  and  after  heat-acid 
hydrolysis  of  the  water  sample. 

In  the  presence  of  organic  nutrients,  both  the  4-  and  20-ppm  TPP 
solutions  ware  completely  transformed  (>95%)  to  phosphate  in  5  days.  In 
the  absence  of  organic  nutrients,  30%  of  the  TPP  was  transformed  (0.3 
ppm-P  consumed)  at  4  ppm  TTP  and  only  10%  (0.5  ppm-P  consumed)  at  20  ppm 
TPP.  No  transformation  of  TPP  was  observed  in  autoclaved  lake  water  at 
4  or  20  ppm.  These  results  indicate  that  TPP  can  be  transformed  rapidly 
(compared  with  chemical  hydrolysis)  in  nonsterile  waters  under  aerobic 
conditions.  It  is  apparent  that  organic  nutrients  accelerate  the 
transformation,  and,  they  may  be  the  limiting  factor  in  the  above  pure 
natural  water  experiment. 

To  confirm  that  the  transformation  is  biological,  the  pond 
microorganisms  were  inoculated  into  trls  buffer-basal  salt  medium 
containing  various  amounts  of  glucose  and  yeast  extract  and  20  ppm  of 
TPP.  The  percent  transformation  of  TPP  after  5  and  7  days  of  incubation 
is  shown  in  Table  20. 


Table  20 


AEROBIC  CONVERSION  OF  TPP  TO  P04-3  AS  A 
FUNCTION  OF  EXTRA  ORGANIC  NUTRIENTS 


Percent 


Nutrients  Added  (ppm) 

Conversion 

to  P04" 

Glucose 

Yeast  Extract 

5  days 

7  days 

100 

20 

95 

95 

50 

10 

85 

100 

25 

5 

60 

80 

50 

10  (noninoculated 
medium 

0 

0 

The  blotrene format ion  of  TPP  was  also  studied  in  anaerobic  water. 
Sear svi lie  Pond  water  was  placed  in  the  anaerobic  chamber,  purged  with 
nitrogen  gas,  and  incubated  overnight.  Tria  buffer  was  added  to 
maintain  the  pH  at  7.  TPP  (20  ppm)  was  added  to  the  water  with  an 
without  the  addition  of  glucose  and  yeast  extract.  After  6  and  10  days 
of  incubation,  TPP  was  transformed  in  all  cases,  as  shown  in  Table  21. 


Table  21 

ANAEROBIC  TRANSFORMATION  OF  TPP  AS  A  FUNCTION 
OF  ORGANIC  NUTRIENT  CONCENTRATION 


Percent  TPP 


Nutrients  Added  (ppm) _  Transformed 


Glucose 

Yeast  Extract 

6  days 

10  days 

100 

20 

72 

81 

50 

10 

48 

62 

0 

0 

20 

38 

Sterile  Control 

— 

16 

16 

These  results  show  that  TPP  can  be  biotransformed  anaerobically. 
This  transformation  Is  slightly  slower  than  that  observed  under  aerobic 
conditions,  and  the  transformation  rate  Is  dependent  on  the  level  of 
organic  nutrients  present.  The  losses  observed  In  the  sterile  control 
are  probably  artificial. 

The  aerobic  biotransformation  of  TFP  also  occurred  in  soil.  TFP 
(500  ppm)  was  added  to  10  g  of  soil  (dry  weight)  in  test  tubes  and 
moistened  with  3  ml  of  water;  500  ppm  of  TPP  (125  ppm-P)  was  used 
because  the  soil  contained  a  large  amount  (~  350  ppm-P)  of  extractable 
phosphate.  Soil  autoclaved  twice  in  2  days  was  used  as  a  sterile 
control.  Periodically,  test  tubes  were  removed  from  the  incubator  and 
the  soil  was  extracted  with  50  ml  of  0.1N  H2SO4  and  filtered.  The 
residue  was  washed  with  150  ml  of  IN  H2SO4.  The  filtrates  were  combined 
and  analyzed  for  orthophosphate  before  and  after  heat-acid  hydrolysis  of 
the  sample. 

TPP  In  nonsterile  soil  was  hydrolyzed  rapidly,  with  a  half-life  of 
2  days.  In  sterile  soil,  the  rate  of  hydrolysis  was  much  slower,  with  a 
projected  half-life  of  more  than  8  days.  Thus,  it  appears  that  the 
biotransformation  of  TPP  will  occur  readily  in  soil  under  aerobic 
conditions. 

The  biotransformation  of  TPP  was  also  studied  in  anaerobic  soil  in 
the  following  manner.  Anaerobic  soil  was  prepared  by  placing  10  g  of 
soil  in  test  tubes,  flooding  the  soil  with  water  in  an  anaerobic 
chamber,  and  incubating  overnight.  Soil  autoclaved  twice  on  successive 
days  and  placed  in  the  anaerobic  chamber  was  used  as  a  eterlla 
control.  To  the  nonsterile  soil  was  added  1000  ppm  of  TPP  (250  ppm-P) 
and  the  soil  was  Incubated.  Test  tubes  were  removed  from  the  chamber 
periodically  and  analyzed  for  increase  in  phosphate  production.  After 
13  days  of  incubation,  90%  of  TPP  was  biotranformed,  whereas  65%  was 
transformed  in  sterile  soil. 

Under  aerobic  and  anaerobic  conditions  in  waters  and  soils,  TPP  was 
evidently  biologically  hydrolyzed.  Linear  polyphosphates  were  found  to 


be  synthesized  and  hydrolyzed  in  microbial  cells.  Polyphosphates  serve 
as  a  phosphate  reserve  and  storage  center  for  phosphorus  within  cells, 
from  which  rapid  biosynthesis  of  nucleic  acid  and  phospholipids  can  take 
place.  A  number  of  enzymes  that  have  been  isolated  from  living  tissues 
will  hydrolyze  pyrophosphate,  tripolyphosphate,  tetrapolyphosphate,  and 
higher  polyphosphates  (Hooper,  1973) 

Sutton  and  Larson  (1964)  and  Sutton  et  L.  (1966)  found  that  half 
of  pyrophosphate  added  to  samples  of  over  200  soils  hydrolyzed  to 
orthophosphate  in  a  few  days  to  a  few  weeks.  They  found  that  the 
general  level  of  biological  activity  in  the  soil  was  the  factor  exerting 
the  moot  Influence  on  the  rate  of  hydrolysis.  Gilliam  and  Sample  (1968) 
found  appreciable  hydrolysis  of  pyrophosphate  in  soils  that  had  been 
steam-sterilized.  Hydrolysis  in  the  same  soils  without  sterilization 
was  severalfold  faster,  presumably  because  of  microbiological 
activity.  Our  findings  for  TPP  parallel  closely  the  reported  fate  of 
pyrophosphate. 

The  smoke  aerosol  and  residue  not  only  contain  TPP  and 
pyrophosphate,  but  also  higher  homologs  up  to  P22  polyphosphates.  The 
biotransformation  of  these  components  was  Investigated  using  a  WP  burn 
residue  mixture  of  polyphosphates.  These  studies  were  complicated  by 
chemical  transformations,  as  described  in  Section  V.H;  however, 
transformation  rates  far  exceeded  those  expected  for  chemical 
hydrolysis. 

When  the  smoke  residue  concentrate  (0.8  mM  as  phosphate)  was  added 
to  pond  water,  the  water  became  highly  turbid,  and  HPLC  analysis  of  the 
water  showed  that  part  of  the  polyphosphates  had  precipitated.  The 
addition  of  1  mM  EDTA  prevented  precipitation,  indicating  that  metals 
were  Involved  in  complexlng  the  polyphosphates.  The  chemical 
transformations  that  occur  in  sterile  water  are  shown  in  Figure  32. 

(P8.  92) 

In  nonsterlle  pond  water  with  added  glucose  (50  ppm)  and  yeast 
extract  (10  ppm),  the  linear  polyphosphates  with  n  >  3  undergo  rapid 
transformation  (t^j  "  2  days),  with  observable  increases  in  TPP, 
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trimetaphosphate,  and  hexametaphosphate  during  the  first  week  of 
incubation.  The  KPLC  profiles  of  the  smoke  concentrate  at  Day  0  and  Day 
6  are  shown  in  Figure  33.  In  the  pond  water  without  added  nutrients, 
the  results  were  similar  to  those  for  the  sterile  control,  showing 
decreases  in  the  concentration  of  higher  polymers  while  the  lower 
polymers  concentrate  after  2  weeks. 

The  above  microorganisms  were  enriched  and  transferred  into  tris- 
buffered  basal-salts  medium  with  smoke  concentrate,  SO  ppm  glucose,  and 
10  ppm  yeast  extract.  In  this  medium,  pyrophosphate  was  readily 
hydrolysed  (t^  “  3  days),  TPP  concentration  increased  in  the  first 
week  then  decreased  slowly,  the  P^-p^g  polymers  decreased  at  variable 
rates  (tjy2  ~  7  days),  and  trlmetaphosphate  did  not  increase  but 
remained  relatively  persistent.  When  the  glucose  and  yeast  extract 
concentrations  were  increased  to  200  ppm  and  40  ppm,  respectively,  the 
hydrolytic  processes  were  greatly  accelerated  and  the  TPP  and 
trlmetaphosphate  concentrations  were  reduced  90%  and  40Z  of  their 
initial  values,  respectively,  In  three  days.  From  these  results,  it 
appears  that  many  processes  are  occurring  in  natural  water  that  affect 
the  persistence  of  polyphosphates: 

(1)  Metals  can  complex  certain  polymers  and  precipitate  them  from 
solution. 

(2)  Nonbiologlcal  hydrolytic  processes  occur,  leading  to  the 
production  of  trlmetaphosphate. 

(3)  The  biotransformation  of  polyphosphates  is  nutrient-dependent, 
showing  Increasing  rates  with  increasing  nutrient 
concentrations. 

(4)  The  biotransformation  of  pyrophosphate  is  relatively  faster 
than  that  of  the  higher  linear  homo logs. 

(5)  Biotransformation  in  the  smoke  condensate  is  slower  for  the 
cyclic  metaphosphates  than  for  the  linear  polyphosphates. 

Polyphosphates  can  be  a  stored  form  of  phosphorus  nutrient  in 
living  cells  and  can  even  be  synthesized  by  microorganisms  (Hooper, 
1983).  Many  microorganisms  therefore  possess  specific  enzymes  to 
synthesize  and  hydrolyze  polyphosphates.  This  specificity  was 
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Hexamata 


Trimeta 


demonstrated  by  Peck  et  al.  (1973),  who  showed  that  bacteria  that 
readily  utilized  TP?  were  poor  utilizers  of  tetrapolyphosphate  and  those 
bacteria  that  readily  utilized  tetrapolyphosphate  were  poor  utilizers  of 
TPP.  Thus,  in  the  smoke  condensate,  It  is  not  surprising  to  observe 
that  the  components  hydrolyze  at  variable  rates  due  to  the  concentration 
and  availability  of  specific  enzymes. 

For  a  natural  water  body,  we  can  assume  under  the  worBt  case 
condition  that  the  transformation  of  polyphosphates  (P  >  4)  is  limited 
by  chemical  hydrolysis  if  no  polyphosphate-specific  bacterial  enzymes 
ate  present.  The  pseudo-first-order  rate  constants  would  range  from  3  x 
10"®  min"1  to  3  x  10"®  min"1  (tjy^  *  160  t0  1®  days)  (sea  Section 
V.G).  Because  most  bacterial  systems  have  tome  ability  to  hydrolyze 
polyphosphetae  (for  example,  TPP  t^j>  ■  13  days,  Table  21),  we  expect 
that  most  polyphosphates  will  have  half-lives  around  13  days  (kb  -  3.7  x 
10"®  min"1),  with  faster  rates  predicted  based  on  the  level  of  organic 
nutrients. 


E.  Hexametaphosphate 

Because  mataphosphatas  resulting  from  the  hydrolysis  of  the  smoke 
condensate  ware  observed  by  NMR  and  HPLC,  we  evaluated  the  ability  of 
microorganisms  to  hydrolyze  this  class  of  cyclic  polyphosphates. 

Searsvllle  Pond  water  wee  collected  and  buffered  with  trls 
buffer.  To  aliquota  of  the  pond  water  we  added  15  ppm  of  sodium 
hexametaphosphate  (HMP)  alone,  or  15  ppm  HMP  (4.5  ppm-P)  plus  50  ppm 
glucose  (Glu)  and  10  ppm  yeast  extract  (YE),  or  15  ppm  HMP  plus  100  ppm 
glucose  and  20  ppm  yeast  extract.  Autoclaved  pond  water  plus  50  ppm 
glucose  and  10  ppm  yeast  extract  was  used  as  a  sterile  control.  The 
waters  were  incubated  at  25*0  in  the  dark.  Periodically,  samples  were 
removed  and  HPM  was  analyzed  by  the  orthophosphate  colorimetric  method 
before  and  after  heat-acid  hydrolysis  of  the  polymer. 

The  amount  of  HMP  transformed  after  7  and  12  days  of  incubation  is 
shown  in  Table  22. 


Table  22 


HMP  HYDROLYZED  AFTER  AEROBIC  INCUBATION 

WITH  SEARSVILLE 

POND  WATER 

HMP  Hydrolyzed  (%) 

Pond  Water  Sample 

12  days 

Water  alone 

24 

34 

Water  +  50  ppm  Glu  and  10  ppm  YE 

38 

52 

Water  +  100  ppm  Glu  and  20  ppm  YE 

65 

79 

Sterile  water  +  SO  ppm  Glu  and  10  ppm  YE  0  0 


The  above  reeulte  Indicate  that  HMP  la  stable  to  chemical 
hydrolysis  at  neutral  pH,  but  that  biological  transformation  can  occur 
[kh  ■  3.6  x  io"7  min”1,  tjy2  "  3*6  Griffith  and  Buxton,  1965). 

The  rate  of  biotransforaatlon  increased  as  a  function  of  the  amount  of 
extra  organic  nutrients  present  In  the  solution.  To  confirm  that 
microorganisms  were  responsible  for  the  transformation,  the  pond  water 
that  showed  HMP  loss  was  inoculated  Into  sterile  media  containing  trie 
buffer  basal  aalts  (TBSM) ,  15  ppm  HMP,  and  50  ppm  glucosa-10  ppm  YE  or 
100  ppm  glucose-20  ppm  YE.  The  amount  of  HMP  hydrolyzed  after  8  and  12 
days  of  incubation  is  shown  in  Table  23. 


Table  23 

HMP  HYDROLYZED  IN  TBSM  AFTER  INOCULATION 
WITH  POND  WATER  ORGANISMS 


Nutrients  Added 


HMP  Hydrolyzed  (Z) 

8  Days  li  Days 


50  ppm  Glu-10  ppm  YE  68  82 

100  ppm  Glu-20  ppm  YE  81  92 


The  results  show  that  the  transferred  microorganisms  grew  in  the 
media  and  hydrolyzed  HMP.  As  we  have  observed  in  the  past,  the 
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transformation  rata  vaa  fastest  in  the  medium  that  was  highest  in 
organic  nutrients. 

An  experiment  similar  to  that  described  above  was  performed  under 
anaerobic  conditions  in  which  the  buffered  pond  waters  were  purged  with 
N2  ga«  and  placed  in  an  anaerobic  chamber  overnight.  Autoclaved  pond 
water  was  used  as  a  sterile  control.  The  amount  of  HMP  hydrolyzed  after 
5  and  9  days  of  incubation  is  shown  in  Table  24. 


Table  24 

HMP  HYDROLYZED  AFTER  ANAEROBIC  INCUBATION  WITH 
SEAR9VILLE  POND  HATER 


Pond  Hater 


HMP  Hydrolyzed  (X) 
5  Days  9  Days 


Pond  water  alone  22  68 
Pond  water  +  50  ppm  Glu  and  10  ppm  YE  79  100 
Pond  water  (sterilized)  5  6 


These  results  show  that  the  HMP  can  be  biotransformed  under 
anaerobic  conditions  and  that  the  transformation  is  accelerated  by  extra 
organic  nutrients. 

Like  condensed  linear  phosphates,  the  metaphosphates  are  chemically 
stable  at  neutral  pH.  However,  the  letter  appear  to  be  biologically 
hydrolyzed  in  eerobic  and  anaerobic  waters.  The  amount  of  organic 
nutrients  affected  the  biological  activity  of  the  water,  which  in  turn 
affected  the  rate  of  hydrolysis. 

F.  Phosphite 

1.  Screening  Studies 

A  biotransformation  screening  test  of  phosphorus  acid  (H3PO3)  was 
conducted  because  this  compound  was  identified  in  the  smoke 
condensate.  Phosphite  at  5  or  15  ppm-P  (13  or  39  ppm  H3PO3)  was  added 
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to  natural  water  collected  from  Searsvllle  Pond  and  Coyote  Creek, 
California  with  and  without  100  ppm  glucose  and  20  ppm  yeast  extract* 

The  waters  were  incubated  at  25*C  in  the  dark.  Samples  were  withdrawn 
periodically  and  analyzed  for  phosphite  color imetrlcally  after  an 
aliquot  of  sample  was  oxidized  by  bromine  to  form  orthophosphate.  The 
orthophosphate  was  analysed  by  the  vanadomolybdophosphoric  acid 
colorimetric  method  before  and  after  oxidation.  Excess  bromine  was 
removed  before  color  reaction  by  heating  the  sample  in  the  autoclave  or 
In  boiling-water  baths. 

Phosphite  remained  stable  for  15  days  in  Searsvllle  Pond  water; 
then  the  concentration  decreased  from  5  to  1  ppm  at  22  days  and  was  non- 
detectable  (<0.1  ppm)  at  29  days.  In  pond  water  with  added  glucose  (50 
or  100  ppm)  and  yeast  extract  (10  or  20  ppm),  phosphite  remained  stable 
for  22  days,  with  the  concentration  decraasing  to  <0.1  ppm  at  36  days. 
These  results  are  shown  in  Figure  34. 

A  similar  study  was  conducted  in  Coyote  Creak  water  (Figure  35) • 
Phosphite  concentration  decreased  from  15  ppm  to  1  ppm  in  35  days  after 
an  acclimation  period  of  28  daya.  In  this  water,  in  the  Searsvllle  Pond 
water,  the  transformation  wee  delayed  as  a  function  of  the  amount  of 
extra  organic  nutrients  added.  In  fact,  when  the  above  transforming 
organisms  ware  inoculated  into  TBSM  containing  phosphite  and  various 
levels  of  glucose  and  yeast  extract,  no  transformation  of  phosphite  was 
observed  after  23  days  with  either  Searsvllle  Pond  or  Coyote  Creek 
organisms.  Similar  results  were  observed  when  neopeptone  or  nutrient 
broth  was  added  to  the  TSBM-glucose-yeast  extract  medium  and  Incubated 
for  30  days. 

The  Coyote  Creek  water  that  showed  biotransformation  of  phosphite 
was  centrifuged  and  the  sediment  was  divided  into  three  portions t  one 
was  resuspended  in  Coyote  Creek  water;  the  second  wan  suspended  in  TBSM; 
and  the  third  was  suspended  in  TBSM  plus  5  ppm  of  yeast  extract.  To  all 
three  suspensions  was  added  10  ppm-P  of  phosphite.  The  organisms  in  the 
resuspended  Coyote  Cruek  water  transformed  phosphite  to  <  1  ppm-P  in  5 
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PHOSPHITE  (ppm-P) 


Phoiphlte: 

fH”  O  1  -  Starlit  Watar  5  ppm-P 

•  2  -  Nonitarlla  Watar  6  ppm-P 

L-  A  3  -  Nonstarila.  Glucoia  50  ppm  Yaait  Extract  10  ppm  6  ppm-P 

|—  A  4  -  Nonstarila,  Glucoia  100  ppm  Yaait  Extract  20  ppm  6  ppm-P 

l_  □  5  -  Nonatarila,  Glucoia  100  ppm  Yaait  Extract  20  ppm  20  ppm-P 


20 


10  h~ 


0.1  L: 
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FIGURE  34  BIOTRANSFORMATION  OF  PHOSPHITE  IN  SEARSVILLE  POND  WATER 
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PHOSPHITE  (ppm-P) 


Photphlte: 


FIGURE  3B  BIOTRANSFORMATION  OF  PHOSPHITE  IN  COYOTE  CREEK  WATER 


114 


*r**ri 


*#  w 


*4>U*d**»±  ^  V-*>  \ 


days.  No  transformation  was  observed  In  the  other  two  media  after  12 
days  of  Incubation. 

When  the  organisms  from  the  first  suspension  were  reinoculated  into 
Coyote  Creek  or  Searsvllle  Pond  water  along  with  10  ppm-P  of  phosphite! 
phosphite  was  transformed  to  <  1  ppm-P  in  6  days. 

From  these  results,  it  appears  that  phosphite  can  be  biotransformed 
in  natural  water  but  not  in  TBSM.  To  further  investigate  this  effect, 
phosphite-biotransformlng  organisms  wars  Inoculated  into  Coyote  Creek 
water  containing  10  ppm-P  of  phosphite.  Up  to  90X  biotransformation 
occurred  In  5  to  13  days.  When  basal  salts  wore  added  to  the  Coyote 
Creek  water,  the  blotransformatlon  was  slower.  Therefore,  It  appears 
that  instead  of  a  lack  of  some  specific  natural  water  nutrient  in  the 
TBSM,  soma  components  of  the  TBSM  may  Inhibit  the  phosphite 
blotransformatlon  at  the  above  concentration. 

The  blooxldatlon  of  phosphite  to  phosphate  by  microorganisms  has 
been  reported  (Casida,  1959).  However,  the  microorganisms  required 
glucose  as  carbon  source  and  oxidation  ceased  when  the  glucose  was 
exhausted'.  The  finding  that  two  water-source  microorganisms  can 
blotranaform  phosphite  without  addition  of  glucose  and  that  the  addition 
of  glucose  and  yeast  extract  extends  the  acclimation  period  is  an 
interesting  phenomenon  in  the  metabolic  disposition  of  phosphorus 
compounds  by  microorganisms. 

Phosphite  added  to  SRI  soil  at  200  ppm-P  level  did  not  oxidize 
during  42  days  of  incubation.  Alexander  (1977)  reported  that  biological 
oxidation  of  reduced  phosphorouo  compounds  was  evident  when  phosphite 
was  added  to  soil.  Accompanying  the  disappearance  of  the  phosphite  was 
an  increase  in  the  concentration  of  phosphate;  the  reaction  was 
eliminated  when  toluene  wee  added  to  inhibit  e  biological  function. 

2.  Detailed  Studies 

We  conducted  a  detailed  kinetic  study  of  the  blotransformatlon  of 
phosphite.  Microorganisms  from  Coyote  Creek  water  were  grown  in  Coyote 


Crtek  water  containing  10  ppm-P  of  phosphite.  The  organisms  were 
centrifuged ,  washed,  and  resuspended  in  one-fourth  volume  of  Coyote 
Creek  water  containing  10  ppm-P  of  phosphite.  Microbial  mass  was 
determined  by  plating  properly  diluted  cell  suspensions  on  Coyote  Creek 
water-10  ppm-P  of  phosphite  agar  using  the  spread  method.  Grown 
colonies  were  counted  and  expressed  as  colony-forming  units  (CFU). 

At  10  ppm-P  of  phosphite,  the  transformation  rate  appeared  to  be 
earo-order.  The  rate  study  was  repeated  using  9.0,  4.6,  1.5,  and  0.7 
ppm  of  phosphite,  which  yielding  transformation  rates  of  0.117,  0.095, 
0.048,  and  0.037  ppm-P  hr-1,  respectively.  Linaweaver-Burk  analysis  of 
the  data  gives  Equation  6i 

Cn 

+  <l0> 

where  C0  -  the  initial  phosphite  concentration 

v  -  the  transformation  rate 

Vm  -  the  maximum  transformation  rate 

Ka  *  the  substrate  concentration  for  half-maximum 
transformation  rate 

Hence,  the  calculated  Kg  Is  2.7  ppm.  Thus,  first-order  behavior 
will  be  observed  only  balow  2.7  ppm.  From  the  studies  performed  at  1.5 
and  0.7  ppm-P  of  phosphite,  we  obtained  pseudo-first-order  rata 
constants  of  0.108  hr"1  and  0.073  hr"l,  respectively,  with  an  average 
colony  count  of  5.5  x  10®  CFU/ml.  From  the  average  pseudo-first-order 
rate  constant  (0.091  hr"*),  the  calculated  second-order  rata  constant 
was  1.65  x  10"®  ml  org"*  hr"*. 

For  a  natural  water  body  such  an  Coyote  Creek,  in  which  the 
microbial  population  is  5  x  10^  CFU  ml"*,  the  calculated  pseudo-first- 
order  rate  constant  is  8.3  x  10"  hr"*  and  the  expected  half-life  for 
phosphite  is  83  hr. 
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Biosorption 


Biosorption  of  WP  was  conducted  using  four  species  of  Gram-positive 
and  Gram-negative  cultures:  Aeotobacter  beijerlncki  ATCC  19366, 

Bacillus  cereus  ATCC  11778,  Escherichia  coll  ATCC  9637  and  Serratia 
aarcescens  ATCC  13880.  The  organisms  were  grown  in  Tryptic  Soy  broth  at 
25*C  overnight;  at  the  late  logarithmic  growth  phase  they  were 
harvested,  washed,  and  resuspended  in  water.  Equal  optical  densities  of 
the  organisms  were  mixed  together  and  diluted  to  a  desired 
concentration.  The  cells  were  dispensed  in  Teflon-lined,  screen-capped 
test  tubes  and  placed  in  an  anaerobic  chamber.  The  cell  concentrations 
used  were  4.0  and  2.0  mg/ul.  WP  in  water  was  added  from  a  stock 
solution  to  give  a  final  concentration  of  2.0  ppm.  After  a  1-hr 
exposure,  the  cells  were  centrifuged,  and  the  supernatant  and  the  cells 
were  separated,  extracted  with  hexane,  and  analysed  for  P^  by  GC.  From 
the  analysis  of  triplicate  tubes,  the  sorption  coefficient  for  the  cell 
suspension  was  calculated  to  be  136  dk  11  at  4.0  mg/ul  and  183  ±  39  at 
2.0  mg/ml.  The  ovarall  average  sorption  coefficient  was  176  ±  39;  thiB 
value  is  considered  to  be  low  and  not  environmentally  significant. 


H.  Biotransformation  Summary 


The  results  of  the  biotransformation  studies  showed  that  WP  is  not 
readily  utilised  by  anaerobic  aquatic  bacteria.  Polyphosphates, 
including  the  cyclic  metaphosphates,  can  be  hydrolyzed  by  both  aquatic 
and  soil  microorganisms  under  aerobic  and  anaerobic  conditions.  This 
rate  of  hydrolysis  is  related  to  the  amounts  of  organic  nutrients 
present. 


Phoophite  is  also  readily  transformed  to  orthophosphate  under 
aerobic  conditions;  however,  the  transformation  is  inversely  related  to 
the  amounts  of  extra  organic  nutrients.  Half-lives  of  less  than  100  hr 
are  expected  in  most  aquatic  environments.  This  material  is  stable  to 
chemical  transformation;  hence,  microbes  will  be  the  primary  route  of 
transformation. 


Biosorption  of  WP  was  found  not  to  be  significant  for 
microorganisms.  This  result  parallels  the  finding  for  the  sorption  of 
WP  to  soil;  however,  it  would  not  be  predicted  from  the  Kow  valve 
(1200).  Although  WP  was  found  not  to  be  toxic  to  aquatic  bacteria  at 
near  aqueous  saturation  levels,  it  appears  that  this  results  from  the 
material's  poor  partitioning  to  the  cell  wall. 

Thus,  all  of  the  smoke-transformation  products  appear  to  be 
susceptible  to  microbial  transformation  in  sail  and  water.  (Phosphine 
was  not  studied  because  of  its  rapid  partitioning  to  the  atmosphere.) 
These  studies  indicate  that  microorganisms  will  play  a  major  role  in  the 
transformation  of  smoka  products  to  orthophosphoric  acid. 
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VII  PARTICLE  CHARACTERIZATION  AND 
DEPOSITION  VELOCITY  MEASUREMENTS  FROM 
SCALE-MODEL  DEPLOYMENT  OF  WP/F  AND  RP/BR  MUNITIONS 

Scale-model  smoke-deployment  studies  were  performed  to  characterize 
the  smoke  particles  generated  from  RP/BR  and  WP/F  and  to  determine  the 
deposition  velocity  of  particles  from  the  aerosol  phase.  These  data 
then  become  Important  in  the  overall  environmental  assessment  of 
phosphorus  smokes. 

Smoke  munitions  were  ignited  in  a  10'  x  12'  x  12' (H)  chamber 
(Figure  3$).  With  tho  RP/BR  munitions,  one  30-g  batch  of  pellets  was 
ignited  every  five  minutes  for  a  30-min  period.  The  WP/F  munitions 
ignited  on  exposure  to  air  and  therefore  were  not  as  easy  to  control; 
one  67-g  wedge  was  ignited  at  the  start  of  the  run  and  another  67-g 
wedge  was  ignited  10  min  later.  These  methods  of  deployment  led  to  a 
time  variation  in  smoke  concentration.  This  variation  is  less  than 
ideal  in  terms  of  a  laboratory  study;  uncteadlness,  however,  is  probably 
characteristic  of  an  actual  field  deployment. 

Ambient  air  was  drawn  through  the  combustion  chamber  at  a  rate 
between  0.42  m-^/sec  and  0.74  nrfysec.  These  flow  rates  resulted  in 
average  smoke  residence  times  in  the  chamber  between  60  and  90  sec. 

Smoke  generated  in  the  chamber  flowed  out  through  a  sampling  duct  (1'  x 
2'  cross-section;  Figure  36).  Average  flow  velocities  in  the  duct 
ranged  from  5.1  mph  to  9.0  mph.  The  duct  contained  a  pool  of  water  on 
which  phosphorus  smoke  could  deposit.  The  pool  of  water  was  in  a  trough 
(19”  x  3';  area  -  4413  cm^);  during  the  run,  the  trough  was  filled  so 
that  the  water  was  level  with  the  duct  floor. 

The  size  distributions  of  the  smoke  particles  were  measured  with  a 
piezoelectric  (PZ)  cascade  impactor.  The  total  mass  concentration  of 
the  smoke  was  measured  using  the  PZ  impactor  and  a  total  filter 


FIGURE  36  SCHEMATIC  OF  THE  SRI  SMOKE  CHAMBER 

Smoke  drawn  out  of  the  smoka  chamber  by  a  blower  passes  through  a  sampling  duct. 
This  duct  is  equipped  with  laser  smoke  detectors,  several  sampling  ports,  and  a  water 
pool  for  particle  deposition  measurements.  Smoke  is  exhausted  through  a  scrubber. 
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sample.  In  addition,  the  mass  of  munitions  used  to  make  the  smoke  and 
the  residue  remaining  after  the  run  were  measured.  Hence,  the  time- 
averaged  mass  concentration  of  phosphorus  in  the  smoke  could  be 
calculated. 

The  flux  of  phosphorus  smoke  to  the  pool  of  water  was  calculated  by 
analyzing  the  water  for  phosphorus  (ion  chromatography).  The  flux  of 
phosphorus  (grams  per  unit  area  per  unit  time)  could  be  calculated  from 
the  laboratory  concentration  measurement,  the  known  run  time,  and  the 
known  area  of  the  water  pool.  The  time-averaged  mass  deposition 
velocity  could  then  be  calculated  by  dividing  this  flux  by  the 
concentration  of  phosphorus  in  the  smoke. 

Bach  munition  was  burned  with  a  high  and  low  flow  velocity  in  the 
duct  (approximately  10  and  5  mph)  and  with  a  low,  medium,  and  high 
relative  humidity  (approximately  20,  50,  and  80Z).  The  mass  median 
particle  size  varied  from  0.95  to  1.23  pm  for  the  RP/BR  smoke  and  from 
0.83  to  1.96  pm  for  the  WP/F  smoke  (Table  25).  There  was  a  slight  trend 
toward  larger  particles  at  the  higher  relative  humidities,  consistent 
with  the  ability  of  phosphorus  smoke  to  absorb  water. 

Typical  size  distributions  of  these  smokes  are  shown  in  Figures  ST- 
42  ,  These  figures  are  in  order  from  low  to  high  humidity,  with  RP/BR 
first.  The  distributions  tend  to  be  blmodal,  with  the  peak  near  1  pm 
to  2  pm  and  one  peak  near  5  pm. 

The  mass  concentration  of  smoke  was  highly  variable  (Table  26). 

This  feature  resulted  from  the  fact  that  the  PZ  lmpactor  takes  point 
samples  (approximately  1  sec  sample  time).  However,  even  the  filter 
samples  were  not  very  constant,  Indicating  an  actual  unsteadiness  in  the 
smoke  burn  rate  (filter  sample  time  was  2  min).  The  phosphorus  mass 
concentrations  calculated  from  the  difference  between  the  final  and 
initial  weight  of  munition  were  generally  lower  than  the  concentrations 
measured  by  the  PZ  lmpactor  and  the  filter  samples.  This  feature  is 
partially  explained  by  the  fact  that  the  smoke  is  not  all  phosphorus 
(smoke  had  a  phosphorus  weight  fraction  of  34Z). 


Table  25 

DEPLOYMENT  CONDITIONS  AND  MASS  MEDIAN  SIZE 
OF  THE  PHOSPHORUS  SMOKE  PARTICLES 


Run 

No. 

Fuel 

Run 
Time 
*  (sec) 

Duct 

Velocity 

(cm/sec) 

Duct 

Reynolds 

Number 

Ambient 

Relative 

Humidity 

(*) 

Mass  Median 
Smoke  Partial* 
Size 
(|Am) 

10 

R 

2280 

250 

6.78*104 

22 

0.984 

11 

R 

2280 

376 

1.02* 105 

22 

1.13  ±  0.12 

8 

R 

2400 

273 

7.40»104 

44 

0.955  ±  0.05 

4 

R 

1860 

380 

1.03*105 

54 

1.04 

6 

R 

2280 

228 

6*18*104 

86 

1.49 

7 

R 

2280 

228 

6.18*104 

78 

1.17  db  0.07 

5 

R 

1800 

380 

1 . 03*10'’ 

79 

1.23  ±  0.06 

14 

W 

1200 

268 

7.26»104 

26 

1.33  ±  0.15 

12 

W 

1620 

402 

1.09‘105 

25 

0.833 

9 

w 

1620 

246 

6.67»104 

39 

0.955 

15 

w 

960 

402 

1.09*105 

52 

1.32 

16 

w 

1200 

291 

7 .89»104 

75 

1.30  ±  0.11 

13 

w 

1260 

393 

l.07‘105 

72 

1.96 

*Fuel  Is 

designated 

as  R  for 

RP/BR  and  W 

for  WP/F. 
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Table  26 


MEASURED  TOTAL  SMOKE  CONCENTRATION 
AND 

CALCULATED  PHOSPHORUS  CONCENTRATION® 


Run 

No. 

Measured  Smoke 
Concentration 
(PZ  Impactor) 
(«/m5) 

Measured  Smoke 
concentration 

. (W 

Calculated,  Time-Averaged 
Phosphorus 
Concentration 

_ (a/«3) 

10 

0.97 

0.079 

11 

0.164  ±  0.066 

mm 

0.052 

8 

3.52  A  1.08 

HI 

0.056 

4 

0.702 

- 

0.078 

6 

0.428 

- 

0.081 

7 

0.108  A  0.017 

- 

0.076 

3 

0.404  A  0.137 

- 

0.061 

14 

2.03  ±  0.87 

0.792 

0.119 

12 

2.37 

0.391 

0.063 

9 

0.904 

- 

0.092 

15 

1.42 

0.920 

0.102 

16 

0.154  ±  0.01 

0.310 

0.096 

13 

0.974 

0.484 

0.080 

*Runs  10, 
runs  14, 


11,  8,  4,  6,  7,  and  S  ware  performed  with  RP/BR  and 

12,  9,  15,  16,  and  13  were  performed  with  WF/7. 
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Given  the  large  fluctuations  in  instantaneous  smoke  concentrations, 
the  deposition  velocity  should  be  calculated  using  the  time-averaged 
phosphorus  concentration.  Table  27  shows  the  measured  phosphorus  fluxes 
to  the  water  pool,  the  calculated  mass-average  settling  velocity,  and, 
for  comparison,  the  settling  velocity  of  a  particle  with  a  size  equal  to 
the  mass  median  diameter  for  each  run.  The  deposition  velocities  ranged 
from  0.015  to  0.05  cm/sec  for  RP/BR  and  from  0.07  to  0.019  cm/sec  for 
WP/Fj  there  is  a  trend  toward  higher  deposition  velocity  with  higher 
humidity,  consistent  with  the  trend  toward  larger  particles  with  higher 
humidities. 

The  settling  velocities  of  the  mass  median  particle  size  are  of  the 
same  order  as  the  measured  deposition  velocities.  This  feature 
indicates  that  gravitational  settling  is  dominating  the  deposition  rata 
(as  opposed  to  Browian  or  turbulent  motion*).  We  expect,  however,  that 
the  deposition  velocity  of  phosphorus  measured  in  our  duct  will  exceed 
the  settling  velocity  of  the  mass  median  particle  size  because  our  size 
distributions  showed  a  substantial  mass  peak  at  a  size  (approximately 
3  nm)  significantly  larger  than  the  median  size  and  because  settling 
depends  on  the  square  of  the  particle  size. 


Using  a  standard  equation  for  turbulent  deposition  (Friadlander,  1977, 
p.  116),  we  find  a  deposition  velocity  on  the  order  of  10“*  cm/sec. 
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Table  27 


MEASURED  PHOSPHORUS  FLUXES 
AND 

CALCULATED  DEPOSITION  VELOCITIES 


Phoephorua  Flux 
to  Water  Pool 


6.79*10 


8. 36*10“ 
1.15*10" 
1.31*10" 
1.34*10“ 
1.13*10" 
3.03*10" 


1.02*10“* 

4.34*10~10 

6.54*10“10 

1.21*10“9 

1.84*10"9 

1.14*10“9 


Depoaitlon 

Velocity® 

(cm/sec) 


8.64*10“a 

1.66*10“2 

2.04*10-2 

1.68*10”2 

1.65*10’2 

1.49*10“2 

5.0*10“2 


8.58*10 

6.94*10 

7.10*10 

1.19*10 

1.90*10 

1.42*10" 


Settling  Velocity 
of  Meet  Median 
Particle  Size 
(ca/aac) _ 

3.26*10“3 

4.63*10“3 

3.31*10“3 

3.93*10“3 

8.06*10~3 

4.97*10“3 

5.49*10“3 


6.42*10 

2.52*10 

3.31*10 

6.32*10 

6.13*10 

1.39*10“ 


*Runa  10,  11,  8,  4,  6,  7,  and  5  were  performed  with  RP/BR  and  rune 
14,  12,  9,  15,  16  and  13  were  performed  with  WP/F. 

“Flux  divided  by  time-averaged  phoaphorua  concentration  (Table  26). 
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VIII  ENVIRONMENTAL  ASSESSMENT 


A.  Aquatic  and  Terrestrial  Modeling 

Tha  laboratory  studies  showed  that  volatilisation,  oxidation,  and 
hydrolysis  processes  would  dominate  the  loss  and  movement  of  W?  in  the 
environment.  These  processes  ware  modeled  for  some  simplified 
conditions  to  determine  whether  more,  complex  modeling  was  necessary. 
Reaction  and  transport  phenomena  associated  with  solid  and  dissolved  P^ 
ware  modeled  for  three  simplified  scenarlost  solid  P4  pieced  in  a 
turbulantly  flowing  stream,  solid  P4  placed  in  a  stagnant  water  body, 
and  solid  P4  burled  in  soil.  A  glossary  of  terms  appears  in  Appendix  1 
for  Identification  of  the  mathematical  symbols  used  In  this  section. 

When  solid  P4  la  placed  in  a  turbulantly  flowing  stream,  it  will 
dissolve  at  a  rate  between  1  and  10  pg/cm^~hr.  The  dissolved  P4  may 
then  oxidise  and  hydrolyse  in  tha  water  with  a  half-life  of  42  hr 
(estimated  shortest  half-life)  or  may  volatilise  out  of  the  water.  The 
reported  Henry's  constant  of  1600  torr-L/mole  indicates  that 
volatilization  from  a  turbulent  stream  will  be  rapid.  The  equations, 
boundary  conditions,  and  solutions  to  this  problem  are  shown  in  Table 
28.  A  calculation  performed  assuming  a  stream  4  feet  deep  and  an  eddy 
diffusivity  of  1.0  cm2/sec  resulted  in  a  volatilization  half  life,  t^/2* 
of  48.5  min.  We  believe  that  volatilization  will  therefore  be  the 
dominant  removal  mechanism  of  soluble  P4  from  turbulent  streams. 

If  solid  P^  l§  buried  in  soil,  the  rate  of  oxygen  diffusion  through 
soil  will  determine  the  rate  of  P4  transformation.  The  equations, 
boundary  conditions,  and  solutions  for  this  problem  are  ohown  in 
Table  29.  Assuming  a  soil  porosity  of  0.1  and  a  water  fraction  of  0.3, 
a  simple  calculation  shows  that  10  years  would  be  required  to  completely 
oxidize  a  1-cm-diameter  piece  of  P4  buried  12  cm  beneath  the  surface. 
This  calculation  does  not  account  for  the  possibility  of  an  oxide  layer 
accumulating  on  the  surface  of  the  P4  chunk  and  restricting  O2  flux  to 


Table  28 
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the  unoxldlzed  P^.  Therefore,  the  calculation  represents  the  minimum 
lifetime  of  a  buried  mass  of  P^. 

If  a  solid  layer  forms,  the  problem  can  be  described  as  shown  in 
Figure  43*  A  rigorous  analysis  would  require  numerical  computations  to 
obtain  the  oxidation  race  of  phosphorus.  However,  a  rough  estimate  may 
be  obtained  by  simplifying  the  problem.  If  we  assume  that  oxygen 
diffuses  to  the  phosphorus  only  through  the  vertical  soil  column 
directly  above  the  phosphorus,  the  equation  describing  oxygen  transport 
through  the  soil  is 


,  i2c°' 

•""dp2" 


<  y  <  d 


(12) 


If  the  curvature  of  the  phosphorus  piece  is  neglected,  the  equation 
describing  oxygen  transport  through  the  solid  layer  to  the  unreacted 
phosphorus  is 


d2C, 


ox 


dy" 


r  <  y  <  R 


(13) 


We  further  assume  that  although  r  is  a  function  of  time,  the  rate  of 
change  of  r  may  be  neglected  compared  to  the  rate  of  change  of  the 
oxygen  concentration.  This  type  of  quani-steady-state  assumption  is 
explored  further  for  the  case  of  white  phosphorus  in  water.  The 
solution  to  Equations  12  and  13,  with  appropriate  boundary  conditions, 
is 


**  +  -  M  1  -  0,75  (  3  ^1/3fM  4/3  -  M 

*D  +  ~)(MC  V  D -  {5np)  (MC  MC  ) 

s  ox  c  ox  o 


0.21  C,lr.  »R2  •  §  •  t 


(14) 


Soil  Surface 


Solid  Layer  Formed  by 
Oxidation  Reaction  Products 

LA-4937-44 


FIGURE  43  PHOSPHORUS  OXIDATION  IN  SOIL,  WITH  FORMATION 
OP  A  SOLID  LAVER 
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Figure  44  shows  the  time  required  to  reach  various  stages  of  consumption 
of  the  piece  of  white  phosphorus,  with  and  without  the  formation  of  a 
solid  layer.  The  y  axis  of  Figure  44  gives  the  fraction  of  mass 
unreacted,  and  time  Is  shown  on  the  x  axis.  The  curves  in  Figure  44 
were  generated  from  Equations  11  (Table  29)  and  14.  Table  30  shows  the 
parameters  used  to  generate  the  curves. 

Table  30 

PARAMETERS  USED  TO  GENERATE  CURVES  IN  FIGURE  44 

Parameters  Value 

d  15  cm 

R  0.5  cm  .  . 

Dg  3.8  .10”;  cmZ/Aec 

D._  1.10“®  cm2 /sec 

p  l.l  i/cm*  , 

Ctlr  1.3  *  10  J  g/cnr 

These  results  show  that  the  extremely  long  oxidation  time  of  the 

piece  of  phosphorus  with  solid-layer  formation  is  due  to  the  low  value 
of  DQX.  Although  there  are  no  data  available  on  oxygen  diffusion 
through  solid  layers  of  phospht nous  oxidation  products,  the  available 
data  on  diffusion  through  solids  suggest  that  Dox  should  be  between  10”® 
and  10"12  cm2/eec.  Therefore,  if  a  solid  layer  forms  on  a  piece  of 
white  phosphorus  burled  in  the  land,  the  lifetime  will  be  extended  from 
about  10  years  to  about  10,000  years. 

In  many  cases,  red  phosphorus  may  be  buried  in  the  form  of  porous 
material,  the  same  type  used  in  experiments  reported  earlier.  That 
material  has  a  large  internal  surface  area  caused  by  the  porous 
structure  such  that  each  gram  of  material  has  1,500  cm2  of  surface. 


1  10  102  103  104 
TIME  (yean) 

LA-4837 -48 

FIGURE  44  SIMULATED  WHITE  PHOSPHORUS  OXIDATION 
IN  SOIL 


The  proportionality  factors  a  equals  1,500  cm^/g  for  the  material  that 
SRI  used  for  the  experiments*  A  mass  balance  on  a  porous  piece  of  red 
phosphorus  gives  the  following: 

dM 

ir  "  ~kaCo2  ‘  Mc  (l6) 

Solving  for  Mc  as  a  function  of  time  gives 

-k*c0  .  t 

-  "eo  *  <l7> 

The  factor  CQ  is  the  concentration  of  oxygen  at  the  phosphorus 

surface  and,  for  b&ried  material,  may  be  less  than  the  air 

concentration.  C  is  found  by  equating  oxygen  flux  through  the  soil  to 
°2 

the  rate  of  oxygen  consumption  by  the  phosphorus.  An  approximate 
relationship  is 


°-21  c«ir  *  ps 
kadR  ' 


T 


+  D 


(18) 


The  experiments  performed  at  SRI  for  red  phosphorus  oxidation  ware 
conducted  over  a  relatively  short  period  of, time  and  essentially 
measured  the  Initial  degradation  rata.  From  section  IV, 


lim 
t  ♦  0 


-1.7952  •  10“8  g/sec 


(19) 


Since  the  red  phosphorus  was  not  burled  in  the  experiments,  d  » 

C  ■  0.21  C.j..  The  reaction  rate  constant  can  be  calculated 
o_ 

the  initial  mass  M_  was  2  g. 


0,  and 
because 
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k  »  2.397  .  l<f8  cm/sec 


(20) 


Th«  time  required  for  95%  of  «  1-em-diameter  piece  of  red  phosphorus  to 
oxidize  if  it  is  buried  15  cm  below  the  soil  surface  is  approximately  11 
years,  assuming  D§  -  5.8  .  10”*  cm2/s.  As  with  white  phosphorus,  if  an 
oxide  layer  forms  on  the  red  phosphorus  surface  due  to  an  accumulation 
of  reaction  products,  the  lifetime  could  be  extended  to  thousands  of 
years. 

If  nonporous  red  phosphorus  is  burled  in  soil,  the  lifetime  will 

also  extend  into  the  thousands  of  years.  A  solid,  spherical  piece  of 

2 

phosphorus  has  a  surface  area  of  4xR  and  will  decay,  without  formation 
of  a  solid  layer,  according  to  equation  21: 
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A  solid  piece  of  red  phosphorus  initially  1  cm  in  diameter 
(Mg  «  0.524  g)  will  require  2,600  years  for  complete  oxidation  without 
formation  of  a  solid  layer.  The  depth  of  burial  does  not  appear  in 
Equation  19  because  the  rate  of  oxygen  consumption  Is  very  low  for  solid 
red  phosphorus  and  the  reaction  rate  controls  the  degradation  process. 

An  analysis  performed  for  the  presence  of  a  solid  layer  predicted  much 
longer  oxidation  times,  as  expected.  That  analysis  was  more 
complicated,  requiring  numerical  solution,  and  it  is  not  presented  here. 

If  a  piece  of  phosphorus  is  placed  in  a  pool  of  still  water,  a 
solid  layer  may  or  may  not  form.  The  scenario  of  white  phosphorus  with 
no  solid-layer  formation  results  in  a  relatively  short  lifetime.  In 
this  case,  a  process  of  dissolution  followed  by  simultaneous  diffusion 
and  reaction  may  proceed,  leading  to  eventual  consumption  of  the  piece 
of  white  phosphorus.  A  quasi-steady-state  analysis  was  performed  for 
this  case,  in  which  the  concentration  profile  of  dissolved  phosphorus  is 
assumed  to  attain  its  steady-state  value  at  all  times  as  the  piece  of 


phosphorus  dissolves •  Hess  balances  on  dissolved  and  solid  phosphorus 
yield  two  differential  equations: 
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where  C  ■  P^  concentration  In  water,  and 
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where  R  ■  radius  of  piece  of  phosphorus  in  water 
and 

C  -  C*  at  r  -  R 
C  **  0  as  r  ■*  “ 

R  "  Rq  at  t  ■  0. 

Equation  22  may  be  solved  for  the  steady-state  concentration  profile, 
thus  obtaining  the  concentration  gradient  at  r  ■  R  for  use  in  Equation 
23.  The  result  gives  the  radius  R  of  the  place  of  solid  phosphorus 
implicitly  as  a  function  of  time; 
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The  steady-state  concentration  profile  obtained  from  Equation  24  is 


(25) 
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A  aphere  of  influence  around  the  piece  of  phosphorus  may  be  defined 
as  the  distance  at  which  the  phosphorus  concentration  in  the  water  drops 
to  a  predetermined  fraction  of  the  solubility  value,  C*.  That  distance 
will  shrink  as  the  size  of  the  piece  of  phosphorus  decreases  over 
time.  Figure  45  shows  the  sphere  of  influence  as  a  function  of  time. 

In  this  figure  the  sphere  of  influence  is  arbitrarily  defined  as  the 
position  where  C/C*  -  0.01.  Equation  24  was  used  to  find  the  radius  R 
of  a  piece  of  phosphorus  at  various  times.  At  each  of  those  times,  R 

it 

was  used  in  Equation  25  to  find  the  distance  r  at  which  C/C  ■  0.01. 
Figure  45  shows  that  the  radius  of  the  sphere  of  Influence  does  not 
significantly  decrease  during  the  first  half  of  the  lifetime  of  a  piece 
of  white  phosphorus.  In  fact,  the  sphere  of  Influence  is  about  75%  of 
its  original  value  after  90%  of  the  lifetime  of  phosphorus  has 
elapsed.  The  value  of  tha  reaction  rate  constant  and  the  diffuslvity  of 
phosphorus  in  water  affect  the  size  of  the  sphere  of  Influence  through 
the  dimensionless  parameter  X. 

The  use  of  dimensionless  variables  often  simplifies  the  solution 
procedure  for  differential  equations.  The  dimensionless  variables  uued 
in  the  analysis  of  the  sphere  of  Influence  around  a  piece  of  white 
phosphorus  allow  a  very  general  solution  to  be  represented  in  a  single 
plot,  Figure  45.  Figure  45  represents  the  solution  for  all  values  of 
the  initial  particle  size,  RQ,  the  diffuslvity  of  phosphorus  in  water, 

D,  the  concentration  of  soluble  at  the  particle  surface,  C*,  and  the 
particle  density,  p.  If  a  plot  of  the  aphere  of  Influence  vs.  time  had 
been  presented  in  dimensional  units  (simply  r  vs.  t),  values  of  RQ,  D, 

A 

C  ,  and  p  would  have  been  specified.  Since  R0  may  change  from  case 
study  to  case  study,  and  D  and  C*  are  not  known  precisely,  specifying 
those  parameters  to  present  a  plot  of  the  solution  would  reduce  the 
utility  of  the  plot. 

Any  case  study  can  be  conveniently  evaluated  with  Figure  45  using 
the  moat  up  to  date  values  of  D,  C*,  and  p.  For  example,  to  analyze  the 
case  of  particles  initially  0.5  cm  in  radius,  tha  quantities  K,  D,  C*, 
and  p  must  be  obtained.  Data  presented  earlier  indicate  the  half-life 
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of  dissolved  In  witer  Is  approximately  42  hrs,  giving  a  value  of  K  ■ 
4.17  .  10“6  a-1 .  The  solubility  of  in  water  io  approximately  4  ppm 

giving  C*  -  4.10"®  g/cm^.  The  density  is  about  2  g/cnr*.  The 
dlffualvity  of  solubile  in  water  is  not  precisely  known  but  assume  D 
■  10”4  cm^/s  for  the  purposes  of  this  example.  The  dimensionless 
parameter  K  can  be  immediately  computed, 

kR  2 

K  -  -  1.04 


A  reasonable  question  to  ask  is,  how  large  Is  the  sphere  of 
influence  after  24  hrs.  First,  calculate  the  value  on  the  abscissa  of 
Figure  45  at  T  ■  24  hrs  (convert  24  hrs  to  seconds). 
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From  Figure  45,  the  line  labeled  K  *  1  approaches  a  value  of  about  4  on 

c* 

the  ordinate  for  small  values  of  T~,  therefore,  at  T  ■  24  hrs,  the 

o 

size  of  the  sphere  of  Influence  is  4  times  the  original  particle  radius, 

a 

or  2  cm.  After  2500  days,  or  60,000  hra,  T  •  C  /p  ■  0.18  and 
r/RQ  «  3,  leading  to  a  sphere  of  influence  of  1.5  cm. 

For  this  specific  set  of  parameters,  K0,  D,  C*,  p,  and  K,  a 
dimensional  graph  could  have  been  constructed.  The  utility  of  the 
dimensionless  graph.  Figure  45,  is  illustrated  with  a  second  example. 

To  analyse  the  case  of  particles  initially  0.35  cm  in  radius  the  same 
procedure  is  used.  In  addition,  postulate  that  dlffualvity  measurements 
made  in  the  interval  between  the  first  analysis  and  this  example  reveal 
that  D  ■  5.10"^  cm^/a  rather  10“^  co^/e  as  originally  assumed.  Now, 
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At  T  "  24  hrs, 


7.2  •  10 


Reading  Figure  45  according  to  the  line  labeled  K  ■  0.1  gives  r/RQ  «  9 
eo  that  the  ephare  of  Influence  is  3.5  cm.  At  T  •  2500  days, 

*£  - 1.« 

P 

In  thla  caee,  r/R  ■  o  because  the  line  for  K  ■  0.1  decreased  to  taro 

ft 

at  TC  /p  ~  0.4. 

Thesa  two  examples  are  meant  to  show  that  A  vide  variety  of  case 
atudlea  can  be  conveniently  analysed  using  Figure  45.  As  additional 
data  become  available  on  the  behavior  of  in  water,  more  precise 
values  of  parameters  such  as  D  and  C*  can  ba  Incorporated  easily. 

A  possible  problem  with  the  above  analysis  Is  the  quasi-steady- 
state  assumption  that  the  concentration  profile  maintains  Its  steady- 
state  shape  as  the  piece  of  phosphorus  dissolves.  For  this  assumption 
to  be  valid,  the  time  required  for  the  concentration  profile  to  reach 
steady  state  must  be  very  small  relative  to  the  lifetime  of  the  pleca  of 
phosphorus.  Figure  45  indicates  a  phosphorus  lifetime  of  about  0.3 

K  ju 

units  on  the  x  axis,  which  corresponds  to  T  ■  1.5  •  10  ,  assuming  C  - 
4  ppm  and  p  ■  2  g/cc.  The  lifetime  of  a  1-om-radius  piece  of 
phosphorus  would  then  be  about  48  years,  assuming  D  ■  10“^  cm^/sec. 

An  analysis  of  the  transient  behavior  of  the  concentration  profile 
was  performed,  and  we  determined  that  the  time  required  to  reach  steady 
state  is  much  less  than  the  phosphorus  lifetime  (i.e.,  much  less  than 
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1.5  x  10*  dimensionless  time  units).  This  analysis  was  performed  by 
solving  the  transient  diffusion  and  reaction  aquation, 
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where 

C  •  C*  at  r  «  R 
C  +  0  as  r  ♦  «* 

C  •  0  at  t  -  0.  <27) 

Initial  condition  27  represents  a  system,  initially  phosphorus-free, 
into  which  a  piece  of  phosphorus  is  placed.  The  solution  to  Equation  26 

is 

Jj  «rfc  (^-)  e‘Kids  +  i  erfc  (§^)  e“KT  (26) 


dimensionless  distance  (29) 

dimensionless  time  (30) 

dimensionless  reaction  rate  constant.  (31) 

Figure  46  shows  concentration  profiles  plotted  at  various  dimensionless 

times  for  K  ■  0.1.  The  concentration  profile  for  T  ■  100  is  identical 

to  the  steady-state  profile.  Therefore,  because  the  lifetime  of  a  piece 

e 

of  phosphorus  is  on  the  order  of  T  *  10  and  steady  state  is  achieved 
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FIGURE  48  TRANSIENT  CONCENTRATION  PROFILES  FOR  K  •  0.1 
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by  T  «  100,  the  quasi-steady-state  analysis  should  be  sufficient  for 

this  problem. 

An  analysis  for  porous  red  phosphorus  submerged  in  water  was 
performed  in  a  similar  manner  to  the  analysis  for  porous  red  phosphorus 
buried  in  soil.  Equations  17,  18,  and  19  are  used  in  this  case  but  C 
is  now  the  oxygen  concentration  in  water,  approximately  8  .  10“®  g/cnr 
at  the  solubility  limit.  The  first-order  reaction  coefficient  for  red 
phosphorus  oxidation  in  water  can  be  calculated  using  the  initial  rate 
data  in  Section  IV  and  Equation  19. 

if  »n  dM  _a  a 

t“  o  <ar>  ■  -5-6335  •  10  ik  ■  -^CV«C  <19) 
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sec 

The  time  required  for  95%  of  a  piece  of  porous  red  phosphorus  to  oxidize 
in  well-aerated  water  is  approximately  20  years.  As  with  red  phosphorus 
in  soil,  a  solid  piece  of  red  phosphorus  will  require  several  thousand 
years  to  completely  oxidize  in  water.  If  an  oxide  layer  forms,  the 
lifetime  of  red  phosphorus  will  be  extended  about  three  orders  of 
magnitude,  whether  it  is  porous  or  solid. 

B.  Aerosol  Modeling 

The  atmospheric  dispersion  of  phosphorus  generated  from  burning 
WP/F  and  RP/BF  to  produce  obscurant  smokes  depends  ou  several 
interrelated  phenomena  that  are  critical  in  evaluating  the  fate  of 
phosphorus  In  the  environment.  A  smoke  plume  is  formed  from  a  munitions 
burn  as  the  smoke  is  carried  by  the  prevailing  winds,  rising  aloft  due 
to  the  bouyant  forces  generated  from  the  combustion  reactions.  Ambient 
wind  speed  and  atmospheric  conditions  that  Influence  plume  rise  and 
expansion  are  important  factors  in  determining  the  downwind  phosphorus 
concentration  and  deposition  of  the  particles  produced  from  a  munitions 
burn  on  the  ground.  The  type  of  terrain  will  also  have  a  profound 
effect  on  the  particle  deposition  of  phosphorous  smokes. 
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In  addition,  aerosol  rate  processes  that  occur  within  the  plume 
affect  the  dispersion  of  phosphorus  in  the  environment.  Nucleation  just 
beyond  the  combustion  zone  provides  the  source  of  very  fine  phosphate 
particles.  These  primary  particles  grow  in  size  by  coagulation  and  the 
condensation  of  water.  Because  deposition  is  size-dependent,  the 
material  deposited  will  change  with  the  evolving  particle  size 
distribution.  As  the  smoke  is  carried  from  the  burn  site,  the 
composition  of  the  particles  is  transformed  by  chemical  reaction; 
phosphate  polymers,  formed  during  the  combustion  reactions,  are 
hydrolyzed  to  phosphate. 

To  model  these  processes  that  affect  the  atmospheric  dispersion  of 
phosphorus,  we  have  Incorporated  relationships  for  the  expansion  of  a 
smoke  plume  and  a  correlation  for  particle  deposition  onto  rough 
surfaces  into  the  MAEROS  computer  code  (Gelbard,  1982),  a  simulation  of 
the  dynamic  aerosol  processes  it  or  multicomponent  systems.  This  overall 
model  was  used  to  predict  the  plume  concentration,  the  fraction  of 
materiel  deposited,  the  change  in  the  particle  size  distribution,  and 
the  change  in  particle  composition  due  to  chemical  reaction  .  s  a 
function  of  distance  downwind  of  the  burn  site.  Simulations  for  the 
phosphorus  dispersion  were  performed  for  several  values  of  each  of  the 
key  parameters  that  Influence  particle  deposition  to  evaluate  the 
sensitivity  of  the  model  to  these  parameters. 

1 .  Theory 

The  fundamental  aerosol  dynamic  phenomena  of  nucleation, 
condensation,  and  coagulation  that  govern  f ine-particle  systems  are 
represented  by  the  general  dynamic  equation  (Friedlander,  1977).  To 
solve  this  equation  under  general  conditions,  Gelbard  and  Seinfeld 
(1980)  Introduced  a  sectionallzed  representation  of  the  particle  size 
distribution;  the  continuous  distribution  of  particla-mass  concentration 
Is  divided  into  discrete  sections  where  the  mass  concentration  In  each 
section  is  constant.  The  rosult  is  a  histogram  of  mass  concentration  as 
a  function  of  particle  size.  With  this  approximation,  the  general 
dynamic  equation  is  reduced  to  a  system  of  ordinary  differential 
equations  that  can  be  solved  by  conventional  numerical  methods.  The 


MAEROS  computer  program  (Gelbard,  1982)  is  based  on  the  sectionalization 
of  the  particle  else  distribution  In  the  general  dynamic  equation.  It 
provides  a  simulation  with  respect  to  particle  size  for  the  evolution  of 
chemical  species  in  an  aerosol  system  undergoing  coagulation, 
condensation,  and  intraparticle  chemical  reaction. 

To  simulate  the  atmospheric  dispersion  of  phosphorus,  the  munitions 
burn  time  is  assumed  to  be  negligible  compared  to  the  evolution  of  the 
plume,  so  that  the  munitions  bum  is  considered  to  be  an  instantaneous 
source  of  smoke.  In  addition,  tho  particle  concentration  is  considered 
to  be  uniform  in  the  horizontal,  x,  and  vartical,  y,  directions  of  the 
evolving  plume  (downwind  distance  it  in  the  Z-direction.)  The 
convective  transport  lr.  the  plume  is  incorporated  into  the  eectlonalized 
dynamics  equations  of  MAEROS  in  aquation  32  (see  Gelbard  and  Seinfeld, 
1980,  Equation  50). 

:*ar "  “  i  re  QAk  ~  ir  ^  ^k +  *QAk  (32) 

+  1/2 A  5-1  r^JtjiQjkQi  “  ^Jti^ikV  “  °Ak  §-1*^1  jQj 
+  °Ak^A  “  5-11  ®AiQAiQAk  ”  °A-l,iQA-l,k3  +  3°A-l,kQA-l 

dQxk 

where  represents  the  rate  of  change  of  the  particle  maas 

concentration  for  section  A  end  component  k.  For  a  system  consisting  of 
a  components  and  m  total  sections,  m  x  s  ordinary  differential  equations 
result. 

The  first  three  terms  in  Equation  32  represent  the  dilution, 

deposition,  and  chemical  reaction,  respectively.  These  are  specific  to 

the  model  for  the  dispersion  of  phosphorous  smokes  and  each  will  be 

discussed  in  detail  later  in  the  report-.  The  4th  and  5th  terraa  are  the 

sectionalized  reprsseutation  of  coagulation  and  the  6th,  7th,  and  8th 
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terms  are  the  sectionalized  representation  of  condensation*  The 
simulation  of  coagulation  and  condensation  for  phosphorous  dispersion 
requires  no  modification  to  MAEROS  and  a  thorough  description  can  be 
obtained  in  the  references  cited. 

The  dilution  of  the  smoke  plume  as  it  expands  is  represented  in  the 
first  term  of  Equation  32.  Relationships  for  the  plume  area.  A,  as  a 
function  of  downwind  distance,  Z,  were  obtained  from  a  least-squares  fit 
of  a  power  law  relationship  to  Paaquill's  empirical  correlations  for 
plume  dispersion  coefficients  reported  in  Seinfeld  (1975). 

The  plume  dispersion  coefficients  are  given  by: 

°i "  w)Yi  (33) 

where  V0  is  the  wind  velocity.  The  fitting  coefficients  A^  end  y^, 
which  have  been  established  empirically  for  three  types  of  atmospheric 
conditions  (moderately  unstable,  neutral,  and  moderately  stable),  are 
presented  in  Table  31. 


Table  31 

PITTING  COEFFICIENTS  FOR  THREE  ATMOSPHERIC  STABILITY  CLASSES 


a  Vertical 

ox  Horizontal 

Dispersion  Coefficient 

Dispersion  Coefficient 

-V  Yy 

Moderately 

stable 

0.134 

0.551 

0.057 

0.928 

Neutral 

0.402 

0.625 

0.164 

0.883 

Moderately 

unstable 

0.0526 

1.14 

0.328 

0.890 

Given  that 

product 

plume,  this 

the  cross-sectional  area  of  the 

the  expression  becomes 

expression  for  plume  expansion 

plume  is  proportional  to  the 

(y  +  Yv)/t  •  For  a  neutral 
x  y 

produces  a  32-fold 
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concentration  dilution  going  from  100  m  to  1000  m  downwind  from  the 
source* 


Particle  loss  by  deposition  on  the  ground  Is  represented  In  the 
second  term  of  Equation  32.  This  term  is  the  product  of  the  ratio  of 
the  horizontal  plume  width,  1^,  to  the  plume  area,  A,  times  the  particle 
mass  concentration,  Q , ^ ,  and  the  deposition  velocity,  • 

The  deposition  velocity  is  defined  as  the  deposltio&  flux,  to  the 
surface  divided  by  the  bulk  plume  concentration  and  represents  an 
overall  mass  transfer  coefficient  that  includes  the  particle  deposition 
mechanisms  of  diffusion,  interception,  impaction,  and  gravitational 
settling.  The  deposition  velocity  depends  on  particle  size.  It  is  high 
for  small  particles  where  diffusional  transport  predominates,  is  minimum 
for  particles  between  0.1  and  0.3  |im  In  diameter,  and  increases  rapidly 
for  larger  particles  as  impaction  and  gravitational  settling 
predominate.  This  general  functionality  is  exhibited  for  all  rough 
surfaces,  but  the  magnitude  of  the  deposition  velocity  depends  on  the 
type  of  surface  and  can  differ  by  several  orders  of  magnitude. 


The  general  correlation  developed  by  Schack,  Pratsinls  and 
Frledlander  (1983,  in  press)  was  used  to  compute  the  deposition  velocity 
in  the  model  for  the  dispersion  of  phosphorus.  The  correlation  was 
obtained  from  a  fit  of  the  available  experimental  data  to  various  rough 
surfaces.  Profiles  of  the  deposition  velocity  as  a  function  of  particle 
size  are  presented  in  Figure  47  for  the  three  natural  surfaces 
considered,  rye  grass  and  crushed  gravel  (B),  and  water  (C).  These 
profiles  were  generated  for  an  ambient  wind  velocity  of  3.0  m/sec  and 
Include  gravitational  settling. 


The  sectional  representation  of  the  deposition  velocity,  derived 
following  the  procedure  of  Gelbard  and  Seinfeld  (1980),  becomes: 
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dv 


(34) 
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The  correlation  for  the  deposition  velocity  is  integrated  over  a 

particle  mass,  v,  mass  of  section,  1,  so  that  the  sectional  deposition 

velocity,  V.  ,  represents  an  average  deposition  velocity  for  section  X. 
aX 


We  now  discuss  the  chemical  reaction  term  (third  term)  in  Equation 
32.  The  change  in  particle  compos tion  due  to  the  hydrolysis  of 
phosphate  polymers  is  modeled  as  three  consecutive  first-  order 
reactions  with  a  phosphate  monomer  formed  at  each  step: 


P3  +  P1 


(35) 


The  change  in  chemical  composition  within  the  particle  due  to 
lntrapartlcle  chemical  raact^n  can  be  related  to  the  change  in  the 
composition  of  the  aerosol,  through  the  sectlonallzed  rate 

coefficient,  which  is  in  general  a  function  of  particle  size. 

However,  when  it  is  assumed  that  the  lntrapartlcle  reaction  rate  is  a 
linear  combination  of  first-order  terms  and  that  the  density  of  each 
particle  remains  constant,  the  sectlonallzed  rate  coefficient  equals  the 
lntrapartlcle  rate  coefficient,  k^,  and  is  independent  of  particle  size. 

The  system  of  first-order  ordinary  differential  equations  that 
result  from  Equation  32  provides  a  general  model  for  the  atmospheric 
dispersion  of  phosphorous  rmokes.  The  model  incorporates  particle 
dilution  in  an  expanding  plume,  deposition  to  rough  surfaces,  and 
lntrapartlcle  chemical  reaction.  It  also  Includes  the  dynamic  particle 
growth  processes  of  coagulation  and  condensation.  The  growth  processes 
influence  the  dispersion  of  phosphorus  because  particle  deposition  is 
such  a  strong  function  of  particle  size. 


The  simulations  of  phosphorous  smoke  dispersion  were  initiated 
1  min  downwind  from  the  source,  the  time  that  corresponds  to  the  age  of 
the  munitions  smoke  sampled  in  the  scale-model  deployment  studies.  The 


Initial  particle  size  distribution  used  In  the  simulation  was  obtained 
from  these  experimental  measurements  of  munition  burns .  Measurements 
were  obtained  for  burns  conducted  at  three  humidities.  For  the 
simulation,  we  assumed  that  water  condensation  was  completed  In  the 
60-sec  period  prior  to  sampling;  consequently,  It  was  not  Included  In 
the  simulation.  The  effect  of  relative  humidity  on  the  phosphorus 
dispersion  Is  therefore  reflected  only  in  the  Initial  particle  size 
distribution. 

For  the  base-case  simulation,  we  chose  Initial  values  of  the  key 
parameters  that  were  typical  of  conditions  that  would  exist  In  a  field 
munitions  burn.  The  initial  particle  distribution  was  taken  to  be  a  log 
normal  distribution  having  a  mass  median  diameter  of  1.3  tin  and  a 
standard  geometric  deviation  of  1.8.  These  values  were  consistently 
measured  for  the  particle  distribution  from  the  WP/F  burns.  The  initial 
concentration  of  0.3  x  10"^  kg/m^  was  used  in  the  simulation.  This 
concentration  is  reported  to  be  required  to  produce  an  effective 
obscurant  (Berkowltz  et  al.,  1981).  Neutral  plume  expansion 
coefficients  (Table  31)  were  used  for  the  base-case  simulation.  The 
other  metarologlcal  variables  that  affect  plume  dispersion  wind  spaed 
and  plume  rise  were  chosen  to  be  3*0  m/sec  and  10  m,  respectively.  The 
munltions-bum  rate,  determined  from  the  initial  concentration,  wind 
speed,  and  plume  expansion  by  JQ  «  C0VQA0,  was  computed  to  be 
6.05  x  10“ 1  kg/sec.  Finally,  the  deposition  velocity  was  computed  for  a 
gravel  surface. 

Profiles  of  the  particle  concentration  In  the  plume  as  a  function 
of  downwind  distance  for  the  moderately  stable,  the  neutral,  and  the 
moderately  unstable  plume-expansion  conditions  are  presented  In  Figure 
48.  The  decrease  in  mass  concentration  is  a  result  of  both  dilution  and 
particle  loss  by  deposition.  To  obtsln  the  fraction  of  material 
remaining  In  the  aerosol,  the  mesa  flow  rats  vslocity  downwind  from  the 
source  is  normalized  to  the  initial  source  rate.  The  fraction  deposited 
is  computed  directly  from  the  fraction  remaining.  The  fractional 
representation  is  Independent  of  concentration  dilution  and  of  the 
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source-dependent  quantities  of  munitions  burn  rate  and  burn  period.  For 
the  base-case  simulation  of  a  neutral  plume,  the  fraction  of  material  in 
the  aeroaol  and  the  fraction  deposited  are  given  in  Figure  49.  The 
particle  deposition  is  rapid  early  in  the  simulation  but  decays  with 
distance  as  the  particle  concentration  is  diluted  in  the  expanding 
plume.  Under  these  conditions,  682!  of  the  phosphorus  burned  remains  in 
the  aeroaol  10,000  m  from  the  burn  site. 

The  sectional  representation  of  the  initial  and  final  particle  size 
for  the  simulation  are  presented  in  Figures  50  and  51.  Section  number, 
plotted  on  the  x  axis,  is  a  function  of  the  log  of  the  particle  diameter 
and  covers  the  diameter  range  of  0.01  to  40.0  fim.  Mass  concentration  is 
plotted  along  the  y  axle.  The  initial  particle  size  distribution, 

Figure  50,  is  generated  from  a  log  normal  particle  size  distribution 
with  mass  median  diameter  of  1.3  |im,  a  standard  geometric  deviation  of 

l. 8,  and  an  initial  concentration  of  3.0  x  10“^  kg/rn^.  The  final 
particle  size  distribution  was  obtained  for  the  simulation  at  10,000 

m.  Particle  concentration  in  Sections  13,  14,  15,  and  16  are  depleted 
in  the  final  distribution  relative  to  the  initial  distribution  because 
of  the  preferential  deposition  of  the  larger  particles.  Coagulation 
appears  to  have  no  significance  in  this  case  because  there  is  observable 
difference  in  the  relative  quantities  of  the  lower  sections,  8,  9,  10, 
between  the  initial  and  final  distributions. 

Experimental  measurements  of  the  composition  of  phosphorous  smokes 
yielded  polymers  that  contained  up  to  16  phosphorous  molecules.  To 
simulate  the  hydrolysis  of  the  polyphosphates  in  the  phosphorous  smokos, 
the  particle  composition  was  rapresentad  by  four  components— -the 
phosphate  monomer,  dimer,  trlmer,  and  polymer,  where  polymers  greater 
than  3  were  grouped  in  the  polymer  component .  The  experimentally 
measured  composition,  used  as  the  initial  values  for  the  simulations, 
were  as  follows:  ■  24.7721;  P2  ■  24.822!;  P3  ■  11.442!;  and 

Pn  ■  38.972!.  The  reaction  rate  coefficients  were  determined  from 
measurements  of  the  hydrolysis  in  aqueous  solution.  Values  for  the  rate 
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FIGURE  50  INITIAL  PHOSPHORUS  SMOKE  AEROSOL/SIZE  DISTRIBUTION 
FOR  C,  =0.3x10'3 


FIGURE  51  FINAL  PHOSPHORUS  SMOKE  AEROSOL  PARTICLE  SIZE  DISTRIBUTION 
AT  10,000  m 


coefficients  used  in  the  simulations  are  as  follows:  k^  -  4.9  x 
10^  sec"1. 

The  evolution  of  the  particle  composition  as  the  munitions  aerosol 
is  carried  downwind  is  presented  in  Figure  52.  The  composition  values 
are  reported  in  mass  fractions.  The  reaction  kinetics  for  the  full  four 
days  is  presented  in  Figure  53.  If  it  is  assumed  that  the  reaction 
kinetics  are  the  same  for  the  particles  deposited  on  land  as  for  the 
suspended  particles,  then  this  96-hour  simulation  of  the  hydrolysis  of 
phosphate  polymers  accurately  models  the  fate  of  these  compounds  in  the 
environment  regardless  of  whether  they  are  in  the  aerosol  or  the 
deposited  form. 

Further  analysis  of  the  phosphorus  dispersion  model  was  conducted 
to  determine  the  sensitivity  of  particle  deposition  to  several  of  the 
key  parameters  that  Influence  phosphorus  dispersion.  The  evaluations 
ware  performed  varying  the  selected  parameter  while  the  base  conditions 
were  maintained  for  the  others. 

The  results  for  the  effect  of  initial  concentration  on  particle 
deposition  is  presented  in  Figure  34.  The  fraction  of  material 
deposited  is  Independent  of  concentration  until  the  initial 
concentration  la  high  enough  for  coagulation  to  be  significant. 
Deposition  is  enhanced  as  the  particles  coagulate  and  grow  in  size.  The 
effect  of  coagulation  on  the  final  particle  distributions  at  10,000  m 
downwind  is  shown  in  Figure  55  for  an  initial  concentration  of  1.0  x 
10"3  kg/m3  and  in  Figure  56  for  an  initial  concentration  of  1.0  x  10"1 
kg/m3.  In  Figure  55,  material  is  preferentially  depleted  in  the  lower 
sections  8,  9,  10,  and  11,  whereas  Figure  56  shows  s  pronounced  shift  in 
the  distribution  toward  the  right  to  the  larger  sized  particles.  These 
concentrations  at  which  coagulation  is  significant  are  higher  than  what 
is  preferred  in  the  field  since  concentrations  of  1 .0  x  10"3  kg/m3  have 
been  reported  to  have  undesirable  health  effects  (Berkowitz  et  al. 

1981).  Therefore,  we  conclude  that  coagulation  is  an  unimportant 
mechanism  in  the  dispersion  of  phosphorous  smokes  except  perhaps  in  the 
immediate  vicinity  of  the  burn  site. 
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FIGURE  53  SIMULATED  PHOSPHATE  HYDROLYSIS  IN  SMOKE  AEROSOL 
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FIGURE  54  PHOSPHORUS  SMOKE  AEROSOL  PARTICLE  DEPOSITION 
AS  A  FUNCTION  OF  INITIAL  CONCENTRATION 

(a)  C0  -  1.0  *  10‘3  kg/rn3;  (b)  C0  -  0.3  x  10'3  kg/m3; 

(c)  C0  -  1.0  x  10*3  kg/m3 ;  (d)  C0  -  1.0  x  10'2  kg/m3; 

(a)  C0  -  1.0  x  10*1  kg/m3. 
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FIGURE  55  FINAL  PHOSPHORUS  SMOKE  AEROSOL  PARTICLE  SIZE  DISTRIBUTION 
FOR  C  =  1.0  x  10-3  Icg/m3 
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FIGURE  56  FINAL  PHOSPHORUS  SMOKE  AEROSOL  PARTICLE  SIZE  DISTRIBUTION 
FOR  C  =  1.0  x  10'1  kg/m3 


Simulations  were  performed  varying  the  coefficients  for  plume 
expansion,  the  wind  speed,  and  the  plume  rise  to  evaluate  the  effect  of 
these  atmospheric  conditions  on  the  phosphorus  dispersion  In  the 
environment.  The  results  of  the  simulation  for  three  conditions  of 
plume  dispersion  are  given  in  Figure  57.  The  moderately  unstable  plume 
is  dispersed  more  rapidly  than  the  base-case  neutral  plume,  and  more 
material  remains  In  the  air.  The  moderately  stable  plume,  however,  is 
dispersed  more  slowly,  and  nearly  502  of  the  material  Is  deposited  by 
10,000  m  downwind.  The  initiation  of  the  deposition  in  the  stable  plume 
Is  delayed  until  the  plume  has  expanded  to  intercept  the  ground. 

Results  of  particle  deposition  as  a  function  of  wind  velocity  are 
presented  in  Figure  58.  Curves  were  generated  for  wind  speeds  of  3.0, 
6.0,  9.0  m/sec.  Although  high  wind  speed  dilutes  the  initial 
concentration  of  the  plume,  particle  deposition  increased  ns  wind 
velocity  Increased  because  of  the  enhancement  of  the  deposition 
velocity.  For  the  high  wind  speed,  442  of  the  burned  material  is 
deposited  by  10,000  u. 

The  effect  of  plume  rise  on  the  aerosol  behavior  is  presented  in 
Figure  59.  Simulations  ware  performed  for  plumes  that  wars  carried 
aloft  10,  20,  and  30  m.  Particle  deposition  is  delayed  until  the  plume 
has  expanded  to  intercept  the  ground;  then  roughly  similar  deposition 
profiles  are  produced.  Plume  rise  has  little  effect  on  the 
overall  particle  deposition;  the  difference  in  the  deposition  fractions 
between  a  10-m  and  a  30-m  plume  rise  is  approximately  82. 

The  type  of  surface  and  its  roughness  characterlstlis  has  a 
profound  effect  on  the  deposition  of  particles.  Expressions  for 
deposition  velocity  can  be  derived  from  basic  principles  only  for  simple 
roughness  elements.  These  results  can  be  applied  to  represent  an 
actual-surface  through  a  surface-specific  density  factor.  When  this 
approach  is  taken,  it  is  found  that  deposition  velocity  profiles  differ 
by  as  much  as  an  order  of  magnitude  for  different  types  of  grasses 
(Davidson  et  al.,  1982).  The  general  correlation  for  particle 
deposition  developed  by  Schack  et  al.  (in  press)  from  an  analysis  of 
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AEROSOL  MASS  FRACTION 


FIGURE  57  PHOSPHORUS  SMOKE  AEROSOL  PARTICLE  DEPOSITION 
AS  A  FUNCTION  OF  PLUME  EXPANSION 

(a)  Modarataly  unatabla;  (b)  Neutral;  (c>  Moderately  ttable. 


t 


AEROSOL  MASS  FRACTION 


FIGURE  59  PHOSPHORUS  SMOKE  AEROSOL  PARTICLE  DEPOSITION 
AS  A  FUNCTION  OF  PLUME  RISE 

(•)  30-m;  (b)  20-m;  (c)  10-m. 
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experimental  literature  produced  deposition  velocity  profiles  for  grass, 
water,  and  gravel  (see  Figure  47)  that  are  similar  in  magnitude.  Yet 
these  relatively  small  differences  in  deposition  velocity  produced 
fairly  significant  differences  in  the  calculated  amount  of  smoke 
depoclted  (Figure  60).  The  differences  in  the  simulation  obtained  for 
grass  (11%  deposition  et  10,000  m)  and  that  obtained  for  gravel  (34%) 
and  water  (29%)  may  be  a  result  of  differences  in  experimental  design. 
The  results  for  gravel  and  water  were  obtains4  from  experiments 
conducted  in  the  same  laboratory. 

The  effect  of  particle  else  on  the  deposition  of  phosphorous  smokes 
was  the  final  parameter  studied.  The  results,  presented  in  Figure  61, 
are  for  simulations  wherein  the  initial  log-normal  distribution  had  mass 
median  diameters  of  0.5  1.3,  and  5.0  (Am.  Particle  size  exhibits  the 
most  significant  effect  on  the  distribution  of  the  material  between  the 
aerosol  and  surface  phases  of  any  of  the  parameters  studied.  These 
results  illustrate  the  importance  of  including  the  particle  growth 
mechanisms  in  the  plume  dispersion  model. 

2.  Summary  and  Conclusions 

The  model  for  the  dispersion  of  phosphorus  munitions  smokes 
incorporates  correlations  of  plume  expansion  and  particle  deposition 
into  a  program  for  multicomponent  aerosol  dynamics  that  includes 
particle  growth  by  coagulation  and  condensation.  The  Intraparticle 
reaction  for  the  hydrolysis  of  phosphate  polymers  was  also  Included  in 
the  model.  Using  this  simulation,  we  predicted  the  downwind  aerosol 
concentration,  the  fraction  of  material  in  the  aerosol  plume,  the 
fraction  of  aerosol  deposited,  the  effect  of  particle  growth  processes 
on  the  particle  size  distribution,  and  the  change  in  particle 
composition  due  to  chemical  reaction. 

The  series  of  successive  first-order  reactiona  used  to  model  the 
phosphate  hydrolysis  reduced  to  a  representation  that  was  independent  of 
particle  size  and  plume  conditions.  Rate  constants  used  in  the  model 
are  based  on  experimental  findings  indicating  that  the  phosphate 
hydrolysis  is  complete  within  four  days. 
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FIGURE  80  PHOSPHORUS  SMOKE  AEROSOL  PARTICLE  DEPOSITION 
TO  VARIOUS  SURFACES 

(a)  Rya  gran;  |b)  Cruihad  gravel;  (c)  Watar. 
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FIGURE  81  PHOSPHORUS  SMOKE  AEROSOL  PARTICLE  DEPOSITION 
AS  A  FUNCTION  OF  MASS  MEDIAN  PARTICLE  DIAMETER 

(«)  dp"  0.6  Jim;  <b|  dp»1.3^m;  (c>  dp«6.0^m. 
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Because  the  atmospheric  conditions  and  the  type  of  terrain  for  a 
munitions  burn  are  variable,  the  effects  of  these  parameters  on 
phosphorus  dispersion  were  studied.  Simulations  varying  plume 
expansion,  wind  speed,  plume  rise,  and  type  of  surface  showed  that  under 
a  wide  range  of  atmospheric  and  terrain  conditions,  at  least  50X  of  the 
aerosol  remained  in  the  atmosphere  at  10,000  nm  downwind. 

The  simulation  was  the  most  sensitive  to  particle  diameter  and  its 
influence  on  deposition  velocity.  The  mass  median  particle  diameter 
produced  from  a  munitions  burn  was  consistently  measured  to  be  between 
1.0  and  1.3  pm  in  diameter.  Therefore,  unless  otherwise  specified,  the 
simulations  were  conducted  using  an  initial  mass  median  diameter  of 
1.3  pm. 

The  sensitivity  of  particle  deposition  to  particle  slice  emphasizes 
the  importance  of  including  the  particle  growth  processes  into  the 
phosphorus  dispersion  model.  Coagulation,  however,  proved  to  be 
significant  only  at  concentrations  higher  than  those  that  are  likely  to 
be  encountered  in  the  field.  Condensation  was  not  included  in  the  model 
presented  in  this  study.  Condensation  occurs  relatively  rapidly  and  it 
was  assumed  to  have  occurred  before  the  initiation  of  the  simulation,  60 
sec  downwind  from  the  burn  site.  The  incorporation  of  condensation  into 
the  model  also  would  have  presented  some  numerical  difficulties. 

Because  it  is  the  most  rapid  rate  process,  it  controls  the  progress  of 
the  numerical  integration  routine.  The  result  is  that  the  simulation 
becomes  inefficient  and  costly.  Further  analysis  of  the  effects  of 
condensation  on  particle  deposition  can  be  investigated  by 
experlmentallly  determining  whether  the  particle  size  distribution  in 
the  plume  stabilizes;  if  it  continues  to  grow,  condensation  should  be 
Incorporated  into  the  model. 

The  correlation  for  the  deposition  velocity  is  perhaps  the  weakest 
component  of  this  model,  as  evidenced  by  the  fact  that  theoretically 
derived  expressions  and  the  experimental  correlations  can  give 
deposition  velocities  that  vary  by  an  order  of  magnitude.  Because  the 
deposition  computed  by  the  simulation  is  sensitive  to  the  deposition 
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velocity,  it  should  be  known  fairly  accurately.  The  correlation 
developed  by  Schack  et  al.  (in  press)  is  the  most  suitable  available 
representation  of  deposition  velocity;  however,  there  are  difficulties 
with  its  application.  Deposition  experiments  are  known  to  scale-up 
rather  poorly,  and  the  results  can  be  dependent  on  the  experimental 
apparatus.  The  deposition  data  used  in  this  model  were  limited  to 
controlled  rough  surfaces  of  a  homogeneous  composition.  A  natural 
terrain  is  heterogeneous,  consisting  of  a  variety  of  surface  types.  SRI 
believes  that  control], ed  laboratory  deposition  studies  on  heterogeneous 
surfaces  need  to  be  performed  so  that  a  more  confident  prediction  of 
deposition  velocity  can  be  obtained. 


IX  DISCUSSION  \ND  CONCLUSIONS 


The  use  of  smoke  In  military  training  exercises  will  deploy  many 
components  into  terrestrial,  aquatic,  and  atmospheric  environments.  The 
major  components  on  a  mass  basis  are  the  condensed  linear 
polyphosphates.  By  applying  solution-phase  chemistry  to  the  aerosol 
particle,  we  projected  that  the  linear  polyphosphates  would  undergo 
acid-catalyzed  hydrolysis  to  orthophosphorlc  acid,  with  a  half-life  of  2 
days.  That  fraction  of  aerosol  that  deposits  on  terrestrial  or  aquatic 
surfaces  where  neutralization  can  occur  will  be  subjected  to  both 
aerobic  and  anaerobic  microbial  transformations,  which  are  nutrient- 
dependent.  In  soils  snd  water,  half-livea  of  up  to  15  daya  are 
projected  in  system  a  containing  at  least  1  x  105  organlsms/ml  (g).  The 
cyclic  metaphosphates,  which  ar«  apparently  formed  by  catalytic  events 
on  linear  polyphosphates  in  natural  waters,  are  much  more  stable  to  acid 
hydrolysis  and  microbial  transformations  than  their  linear  analogs. 
However,  aerobic  and  anaerobic  bacteria  were  found  to  hydrolyze  these 
compounds,  with  half-livss  of  15  to  20  days.  This  transformation  wca 
also  observed  to  bs  nutrient-dependent. 

Phosphorous  acid  (H3PO3)  was  found  to  be  stable  in  aerated  water  but 
could  be  oxidized  to  orthophosphorlc  acid  by  microorganisms.  Half-lives 
of  25  to  35  days,  including  acclimation  periods  of  15  to  25  days,  were 
observed  in  two  natural  wetsr  systems.  This  biotransformation  was  found 
to  be  inversely  related  to  the  amounts  of  extra  organic  nutrients  present 
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in  the  system,  which  tended  Co  slow  the  transformation.  However,  the 
transformation  still  occurred  and  suggests  that  the  persistence  of 
phosphorous  acid  will  be  controlled  by  microbial  populations. 

Elemental  phosphorus  (P/4)  in  the  vapor  phase  will  react  extremely 
rapidly  with  atmospheric  oxygen  to  produce  the  phosphorus  oxides  mentioned 
above.  Thus,  the  overall  fate  of  P^  in  the  vapor  phase  will  be  controlled 
by  oxidation  and  acid  hydrolysis  to  produce  orthophosphorlc  acid.  This 


reactivity  to  oxygen  minimizes  the  possibility  that  any  deposition  in 
terrestrial  or  aquatic  environments  will  occur  from  the  vapor  phase. 

Fhosphine  (PH-j)  was  observed  to  be  produced  from  the  incomplete 
combustion  of  WP/F  and  from  the  solid-phase  oxidation  of  RP/BR  as  a 
function  of  relative  humidity.  Based  on  its  high  vapor  pressure  and  low 
water  solubility!  PH3  will  rapidly  volatilize  into  the  atmosphere.  In 
the  atmosphere,  PII3  will  be  subjected  to  oxidative  transformations 
resulting  from  ozone  and  hydroxyl  radical,  yielding  half-lives  of  8  and 
5  hr,  respectively.  In  the  presence  of  sunlight,  transformation  rates 
are  enhanced  (t^  “2.3  hr)  due  to  interactions  with  tropospheric 
agents  that  are  net  completely  understood.  We  project,  however,  that 
PH3  will  not  be  a  persistant  chemical  in  the  environment. 

The  modeling  of  the  smoke  aerosol  projected  that  at  least  30%  of 
the  aerosol  remained  in  the  atmosphere  at  10,000  a  downwind.  Because 
acid  hydrolysis  predominates  In  the  aerosol  phase,  the  major  component 
of  aerosol  deposition  sfter  20  hr  will  be  orthophosphorlc  acid;  la  the 
immediate  vicinity  of  ths  burn,  the  heavier  particles— including  the 
condensed  polyphosphates— could  contaminate  surrounding  environments. 
Overall,  however,  the  products  of  WP/F  and  RP/RR  combustion  are  readily 
controlled  by  environmental  mechanisms  that  convert  these  materials  to 
orthophosphorlc  acid. 

Incompletely  combusted  WP/F  and  RP/BR  or  disposed  smoke  devices 
from  manufacturing  practices  were  found  to  persist  in  terrestrial  and 
aquatic  environments  for  an  exceptionally  long  time.  The  fate  of  WP/F 
or  RP/BR  burled  in  soil  will  be  controlled  by  the  diffusion  of  oxygen 
through  the  soil,  the  diffusion  of  oxygen  through  surface-oxide  layers 
that  build  up  on  phosphorus,  the  surface  area  of  the  phosphorus,  and  the 
depth  at  which  the  phosphorus  is  buried.  For  a  l-crn-*  piece  of  white 
phosphorus  buried  12  cm  in  soil,  the  lifetime  is  expected  to  range  from 
10  years  In  the  absence  of  a  surface-oxide  layer  to  10,000  years  In  the 
presence  of  a  surface-oxide  layer.  RP/BR  is  expected  to  persist  longer 
due  to  its  slower  reactivity  with  oxygen. 


The  persistence  of  WP/F  in  water  is  controlled  by  the  readiness  of 
the  phosphorus  to  dissolve  into  the  aqueous  phase .  Due  to  the  low  water 
solubility  of  dissolution  is  slow  and  a  lifetime  of  48  years  is 

projected  for  a  1-cm-diameter  piece  of  WP.  On  the  other  hand,  RP/BR  has 
essentially  no  water  solubility  and  oxidation  at  the  water-solid 
interface  will  control  the  lifetime  of  RP.  On  the  basis  of  the  slower 
reactivity  of  RP  with  oxygen  relative  to  WP,  lifetimes  are  expected  to 
be  in  excess  of  48  years. 

The  lifetime  of  soluble  P^  in  the  aquatic  environment  will  be 
controlled  by  volatilization  in  shallow  or  highly  turbulent  water 
bodies,  a  factor  overlooked  by  previous  Investigators,  leading  to  a  wide 
range  of  reported  rate  constants  for  the  loss  of  P4  in  water.  In  a 
stream  1  a  deep,  the  half-life  is  estimated  to  be  50  min.  As  depth 
increases  and  turbulence  decreases,  the  persistence  of  soluble  P4  is 
controlled  by  hydrolytic  and  oxidative  processes.  These  processes  occur 
at  approximately  equivalent  rates,  yielding  a  half-life  of  42  hr.  These 
processes,  however,  can  be  mitigated  by  naturally  occurring  organics 
that  may  provide  stability  to  P4  solutions  in  water.  The  exact  nature 
of  these  stability  mechanisms  is  unknown  and  caution  should  be  used  in 
estimating  persistence  of  soluble  P4  in  deep  water  bodies  with  low 
turbulence. 
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Appendix  A 


NOMENCLATURE  FOR  TERRESTRIAL  AND  AQUATIC  MODELING 


Area  of  phosphorus  surface 

Concentration  of  dissolved  phosphorus  in  water 


Concentration  of  air 


Concentration  of  oxygen  in  soil 


Solubility  limit  of  phosphorus  in  water 
Diffuslvity  of  phosphorus  in  water 
Diffuslvity  of  oxygen  in  soil 
Diffuslvity  of  oxygen  in  solid  layer 
Depth  below  surface  phosphorus  piece  is  burled 
Complementary  error  function  .  , 


Dimensionless  reaction  rate  constant 


Reaction  rata  constant  for  red  phosphorus  oxidation  or 
white  phosphorus  destruction  in  water 


Mass  of  phosphorus 
Initial  mass  of  phosphorus 
Radius  of  piece  of  phosphorus 
Initial  radius  of  piece  of  phosphorus 
Radial  distance 


Radius  of  unreacted  phosphorus  in  a  piece  where 
a  solid  layer  forma 


Dummy  variable  of  Integration  in  Equation  17 


Distance 


Dimensionless  radial  coordinate 


Density 


Dimensionless  time 


Appendix  B 


DEPOSITION  VELOCITY  MODEL 


Deposition  efficiency,  t}  is  defined  es  follows: 

particlea/sec  depositing  on  single  cylindrical  element 
11  ”  particles/ sec  passing  through  projected  cylinder  area  dfh 

Deposition  efficiency  by  diffusion  (Parnas  1984} 

,  QO  «  -1/2  _  -2/3 
•  1.88  Re  Sc 

U  d 

Re  -  “2-r  Sc  -  v/D 


'M.  ; 


where 


velocity 

kinematic  viscosity 
fiber  diameter 
dlffualvity. 


i;- 


Deposition  efficiency  by  Interception  (Parnas  1984) 


0.80  R*  Re 


2  n^l/2 


wha  re 


d  / d, 
P  t 


particle  diameter 
fiber  diameter 


ff. 

Sl’< 

|Vn[ 


The  overall  efficiency  becomes: 


rj  ■  1  -  [(1-rjj)  (l-r^)  - 


> 

•i.: 


i*'N 


i'A 


In  terms  of  the  dimensionless  mass  transfer  number  (Sherwood  No.)  for  a 
cylinder 


kd 

_ e 

D 


J£S£ 


Pe  - 


d.U 

f  o 


The  mass  transfer  coefficients  to  an  individual  cylinder 


[m/sec] 


Particle  deposition  rate  to  singla  roughness  element 

J  -  kA^C  [kg/ sec  -  element] 

where 

Af  *  surface  area  of  element 

Flux  to  roughneso  surface 

Jg  ■  pfkA  C  [kg/m2  sec] 

where 

Pf  -  density  of  roughness  elements  [elements/m2] 

Filtration  deposition  velocity 

V£  -  Jf/C  -  pfkAf  [m/sec] 

Overall  deposition  velocity 


Where  Vg  is  gravitational  settling  velocity 


where 

Pp  -  particle  density 
g  ■  gravitational  constant 
|i  ■  viscosity 

and  the  slip  correction  factor  is  given  by: 

C  -  1  +  ~  <1.37  +  0.4  exp 
P 

where 


X  ■  naan  free  path 


